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ABSTRACT. The three-dimensional development of the nucleon component of extensive 
air showers is investigated. A vectorial diffusion equation is derived for the phase-space 
distribution of the nucleons, and solved as far as required to determine the mean square 
angular deviation, and the mean square distance from the shower axis. Account is taken of 
the variation of density in the atmosphere, and numerical results are obtained for all altitudes. 
These are exhibited graphically. The lateral spread of the nucleon cascade in lead is also 
determined, and the results tabulated. A full discussion is given of the significance of the 
numerical values obtained. 


§1. INTRODUCTION 

NTIL the last few years it was commonly supposed that the primary 
| component of cosmic radiation consisted mainly of electrons. On this 
hypothesis, the development of extensive air showers was explained in 
terms of a cascade of high energy electrons and gamma-rays. The theory of the 
longitudinal and lateral development of such cascades was given by Heisenberg 
(1946) and Moliére (1942) among others. Most subsequent work (Nordheim 
and Roberg 1949, Nordheim 1941) has also been based on the assumption that 
extensive air showers were initiated by electrons. Only recently it has been 
shown experimentally (Critchfield et al. 1950) that the proportion of electrons, 
if any, in the primary cosmic radiation is less than 3%, and that actually protons 
and alpha-particles predominate in number. ‘These nucleons, through their 
collisions with the nuclei of the air, are instrumental in the production of a 
nucleon cascade (Messel and Ritson 1950, Janossy and Messel 1951, Messel 
1951 a,b). Besides protons and neutrons, charged and neutral mesons 
(7-mesons and mesons of greater mass) result from the nucleon collisions, 
among which the neutral z-mesons decay to produce gamma-rays (Carlson 
et al. 1950). 
It is now generally accepted that these gamma-rays, a mere by-product of 
the nucleon cascade, are mainly responsible for the soft component of extensive 
air showers. ‘This circumstance removes the physical basis of previous theories 
of the development of the soft component, a correct treatment of which would 

take into account its continuous evolution from the nucleon cascade. 
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The inadequacy of the earlier theories is particularly evident if one considers 
the lateral spread of extensive air showers. ‘The spread of the soft component 
will in fact be determined primarily by the spread of the nucleon and meson 
components. Strangely enough the first experiments (Cocconi et al. 1949,” 
Greisen et al. 1950) appeared to support the theory of Moliére, and indicated 
that the spread of the penetrating component was practically the same as for the 
soft component. Succeeding experiments from another source (Eidus, Alymova: 
et al. 1950, Vernov et al. 1950, Eidus, Blinova et ai. 1950), however, showed that 
the extent of air showers was much greater than predicted by Moliére, and 
could not be explained by any mechanism hitherto proposed. 

The purpose of this paper is to show that the above discrepancy can be 
resolved by the theory of the lateral spread of the nucleon component of cosmic 
radiation. As a preliminary to this work it was necessary to know the mean 
square angle of scatter of the nucleons involved in nucleon-nucleon and 
nucleon—nucleus collisions. We have already found (Green and Messel 1951 a, b) 
that the mean square angle of scatter of particles resulting from the nucleon— 
nucleon collision should be (U+U’)/UU’ where U and U are the energies of 
the incident and scattered particles measured in proton masses. We found also 
(Messel and Green 1952) that the corresponding mean square angle for a 
nucleon-nucleus collision was proportional to 1/U. ‘These results are rather 
insensitive to the particular form of cross section for nuclear collisions adopted, 
and are a natural consequence of the properties of relativistic transformations. 

Our considerations relate only to particles with energies greater than several 
103 Mey. At these energies, the angular deflection of a particle resulting from a 
nuclear collision is small; we therefore systematically neglect the average value 
of (1 —cos @)?, which varies as the mean fourth power of the deflection 6. This 
is the only approximation employed. In the atmosphere, the density of the air 
varies with height in a manner which depends on the temperature distribution. 
In evaluating our results we have considered an isothermal atmosphere as the 
most faithful simple representation of the physical reality. The result differs 
widely from what would be obtained under the assumption of constant density. 
Previous work on the lateral spread of the soft component in the atmosphere 
assumed a constant density, and is open to criticism on this account alone. 

Results are obtained for absorbers of both constant and variable density, 
applicable to lead slabs, for example, and the atmosphere respectively. The 
mean square angular deviation from the shower axis, and the mean square 
distance from the shower axis are evaluated from various initial conditions. 
Information is thus obtained bearing on the energy spectrum of the primary 
component of the cosmic radiation. ‘This and other experimental implications 
are discussed in the final section. 


§2. THREE-DIMENSIONAL DIFFUSION EQUATION 
In this section a diffusion equation will be derived for the phase-space 
distribution function /, which measures the probability* f(p, r, z) dp dr of finding 
a nucleon with momentum in the range p, dp, at height z, and with a horizontal 
displacement r, dr, from the axis of the shower. The height z, and the two 
components of the vector r are measured in centimetres; the three components 
of the momentum p in units of Mc, where M is the proton mass and ¢ the velocity 

* The differential element dp stands for dp/27p?, i.e. dpydpsdp3/27p?. 
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of light. At the ultra-relativistic energies which we consider, the magnitude p of 
the momentum is indistinguishable from the energy U, measured in units of Mc?. 

Let v(p’,p)dp be the differential probability of finding a nucleon with 
momentum p, dp among the nucleons which result from the collision with a 
nucleus of a nucleon of momentum p’. This will not be required explicitly for 
our present purpose, though it has been determined in our previous work 
(Messel and Green 1952). 

Let 5(z) be the density of the medium in which the cascade develops; then the 
probability that a nucleon will collide with a nucleus of the medium in traversing 
a layer of thickness dz is kd(z)dz/C where C is the cosine of the angle between 
its direction of motion and the vertical, and k is a constant whose value is 
1/75 cm?/g in air, or 1/160 cm?/g in lead. 

Consider a nucleon with momentum p arriving at height z. The probability 
that it has traversed the layer of depth €—sz immediately above, without 
collision, is 


exp[ — {0(z)—6(D)}/C], where 6(z)=h | F Six odcine Inne ates (1) 


The probability that it resulted from a collision between depths ¢ and €+d€ is 
—6'(2)dt/C. 

The probability that its immediate ancestor had momentum p’,dp’ is 
Ft (p’, e, £) dp’ dp, where pe, de was its displacement from the shower axis. The 
probability that it generated a nucleon with momentum p, dp is 1(p’, p)dp- 
Hence 


fip.r,2)=—|_ | AP’, (Pp) exp [-{0(2) -A(}/C] 


6’(f) dt/C dp’ +f(p, p.., ©) exp[—A(2)/C] ...... (2) 
where the components of e are connected with those of r by 
Pi=%—Pi(S—2%)/Ps,  P2=Ta—Po(S—2)/Pp sees (3) 


and g,, is the value of e for €= «0. 
By differentiating (2) with respect to z, 7; and 7,, one obtains the differential 


equation* 


C_ (afte) _ ts 20) _ 22 YO) 4 yp) | pprvie’ map ; 
Fey oe — BO) Pe BP fp) = | eR PAB. (4) 
Since C, the cosine of the angle made by p with the vertical, may be written 
p/p, one obtains, on making the substitution z= —rs, the vector equation* 
1 p o(p) es 
PS 8 = ADI Sd E55 AL 
Tam) par thP)= | APP", p) dp (5) 


It is customary to measure the vertical distance in g/cm”, which may be effected 
by taking 6(z) instead of z as the variable. ‘Then the eqn. (4) becomes 


%_1 (no , mH) a 
cB-a(' ines +f= |/@) p)dp’ ss (6) 
where @’ continues to represent the function 6’(z), expressed as a function of @. 
In an isothermal atmosphere, 
A(z) =(Rpo/8) exp (—89o2/Po), 8 = 88 9/Po ss se (7) 
* The variables on which f(p, r, s) depend are dropped wherever no confusion is likely to arise. 
45-2 


692 H. S. Green and H. Messel 


where p, and 8, represent the pressure and density at sea level, and g is the 
acceleration due to gravity. For a medium of constant density 5, on the other 


hand, hig PSielaed Guy Gaerne Sree (8) 


which is constant; the unit of length may then be chosen so that kb =1. 


§3. EQUATIONS FOR THE ANGULAR AND LATERAL SPREAD 
To obtain the equations satisfied by the mean values of the quantities 71”"72", 
it is convenient to follow Borsellino (1950), by considering the two-dimensional 
Fourier transforms of f with respect to 7, and 79: 


g(p, A, 0) = I | fpr, 9) exp(—iden)drjdry secon (9) 
Then one has £= 2g. (PCAs Ee eae (10) 
where Yi | | (—ir,)™(—irs)"f(p, 0,0) dr, dry. cease (11) 
The equation satisfied by g is obtained from (6): 
og tA. P Te / t —7, 
Oe ap g+e= | a(P )o(P’,P)ap’. sawn (12) 


By differentiating (12) m times with respect to A,, and m times with respect 
to Aj, then setting A; =A, =0, one obtains a set of equations for the g,,,. For 
example, if 


So0=*% —Wio=th, —Wor=Us, —(SootSoe2)=% «+++: (13) 
one has 
0 ~ 
C aa +«(p) = | x(p')(p', p) dps 5 tists ees (14) 
0 “#, 
cAMP) _ Pe ap) +uy(p)= | u(P')(',) dB 
ead 15) 
Ou — 7 ( 
Cc ae ¢ 7 ¥() +u:(Pp) = | us(P’)»(P’, P) dp 
Oz 2. ae 
CP Fe MP) +H0)= | PP PP. ee. (16) 
Writing u, =)p,y/p and u,=p.4/p, the last equation becomes 
Oz 292 aes 
mC ae — Gr H(P) + 2(P) = | P(e, p) dp’j.m =e (17) 
where Si=(p240 Vii a. (18) 


Multiplying the two equations (15) by p,/p and p,/p respectively, and adding, 
one has 
Oy(P) 1 (Pibr' + P2P2’) or fp! 9) ga! 
{CO _ zp) +y(p)} = | PPL EPPe) vpru(p', pap’... (19) 

The equations (14), (17) and (19) will now be written in polar coordinates; 
the following notations will be used. The cosine of the angle between the vectors 
p and p’ will be denoted by ¢, and its sine by s; the angle between the plane con- 
taining both p and p’, and that containing both p and the vertical, will be denoted 


by ¢. Then i dp’ = [- dp' | * de ibs agion 
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and the cosine of the angle between p’ and the vertical is C’ =Cc+ Ss cos ¢; 
also (p1P1' + PeP2')/PP' =(P - P’ — Pps’) /pp’ =c — CC’ = S(Sc— Cs cos ¢). Hence 


fal o rl 20 
Bene) ar Oy x(U", Ce + Ss cos $)(U’, U, c) dU’ de dd /2n 
a0 
U¥ -1/0 


ae (20) 
s{cX9 — Sal gi o} = in {_ [stse-es cos ¢) 
x y(U’, Ce+ Ss cos 6)(U', U,c) dU’ dedd/27_...... (21) 
ee NCU 
= | i | . { Ces cosh LLU: CV AU dcds (on ae (22) 
U0 /-1/0 


where p and p’ have been replaced by U and U’. 

The physical quantities of primary interest are the mean square angular 
deviation and the mean square distance of the particles from the shower axis. 
They are given by 

1 1 1 
wo U, 0) = { | x(U, C, 0).S? ac} {| aU; C0) ac} oes (29)) 
=1 1 


and 12(U, 8) = { { “a U, C, 8) ac} { { ig x(U, C, 8) ach”. eae (24) 


In order to determine these, it is not necessary to solve the equations (20) to (22) 
completely; one requires only certain coefficients in the expansions of x(U, C) 
and 2(U,C) in Legendre polynomials. One therefore writes 


x(U,C)= (K+ })vq(U)Pa(C, y(U,C)= ¥ (K+ dy a(U)Pa(O) 


Ba a (25) 

pe Nae ariel Eat) 
then a(t j=xaixn(U)\—x(U)ia(U) ss.» (26) 
and Lex (a ee) Se Ae ee A Ss (27) 


The coefficients required are determined in the following section. 


§4. DETERMINATION OF THE LEGENDRE COEFFICIENTS 
When the expansions (25) are substituted in the integrals of (20)-(22) the 
latter simplify. Use is made of the addition theorem for Legendre polynomials: 


| hP-(Oe4 Sveoe Basen = P,LO)P 0) 


also ik P,(Cc+ Ss cos ¢) cos ¢ d/27 = Px’(C)Px'(c)/K(K +1) 
‘ 2 IGE 
~ K(K+1) dC dc 


1 
If ra, = { MO Ure\Pale)de aaa (28) 


1 200 
one has | [ x(U’, Cc+Ss cos 6)(U’, U, c) de dp/2r 
0 


by 
1 Mes 
T=3) 


(Kee etl", UNPHey . 2 ..(29) 
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therefore, multiplying (20) by P,(C), and integrating with respect to C from 
—1 to 1, one obtains 


a(K+1 K 

Frey telU)+ pcg tx a(U)) +9(U) 

_ ip: eA eee (30) 
U 


Again, one has by a straightforward but tedious calculation 


rl 1 2 
| dc| | S(Sc—Cs cos ¢)y(U’, Cc+Ss cos ¢)(U’, U, ¢) de d/27 
ae 10 


=Hy(U')—y(U (UU); nn eee (31) 
hence, writing 
1 

wU)=| 9, C)S*dC=8fy(U)—y(W)} veer (32) 

one obtains by integrating (21) with respect to C 

aa (C=1)8%y(U, Cac + SL a 5 Ss Ty al UD 
50 30° a 
x [- Uk Oa ee (33) 
U 


Now, this equation is exact, but the integral of the first term contains a factor 
(1—C)?, the average value of which is negligible, since it varies as the mean 
fourth power of. the deflection from the shower axis. Hence, using (26), (33) 
becomes 
0b(U) 
06 
Similarly one obtains, on integrating (22) with respect to C from —1 to 1, 
dz,(U - : : 
2 220) 30) = (Ne x,(Udvd U2 Uae a eee (35) 


Now it may be seen from as that 2(U, C) is proportional to S?; hence 


~ x9(U)w*(U)/6’ + b(U) = | : AU eA 0 ab ee (34) 


=(U)—a(U)=_a(U, C\C~1)de 


may be neglected; i.e. one may substitute z)(U) for z,(U) in (35). 

We shall first determine the mean square angular deviation, from (30). 
Following the approximation of neglecting the mean value of (1—C)?, we shall 
neglect 1 

( (Uy C(I — CP dC = y+ $9 — 2, 


by setting X =4Nx_+ FX. Tere es) 


From (30) one then obtains by substituting k=0 and k=1, and subtracting the 
results, 


= {%o( U) — x2(U)} + { {%o(U) — x2(U)} =| {X( U") — xa(U")}v9(U', U) dU’ 
: fs to( U'ivg( UW, 0) v0) (37) 


The coefficient 
n(U)= [a U,C)dc=[ | | far dryac 
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is simply the average number of nucleons at depth @ in the medium, and is 
known to be 


(U, 8) = 5 wT, fits (U/Uyy-@+? exp {-H(s)O} ds... (38) 
where " Ms)=1- is (a) U'v(U’, v)d(z) a (39) 


if the shower is supposed to be generated by a primary nucleon of energy Up. 
The equation (37) can then be solved by means of a single Mellin transform; 
the result is 


HU) —m(U}= ra 73[ (Ge) slo) lexp (M0199 
aie AS\e aoe My, Gr Pelehant sone (40) 
where 
{h(s)—h(s— 1}i(s) =3 | ; (7) U'y,(U', U) —v,(U", U)}d () . (41) 


We have here used the result of our previous paper (Messel and Green 1952), 
that U’*{v,(U’, U)—v,(U’, U)} is a function of U/U’ only. We state here, for 
convenience the explicit values of vo(U’, U) and »,(U’, U)—v,(U’, U): 
: 1 1 setio /Y\-6+ 
vo U' ; U) — Daa a ina (z) {1 — f[D 4«(s)]} ds ee eeee (42) 
and 
${vo(U', U)—v,(U’, U)} 


=O (oa ea} PD eOl-Daals— 1) 


~ Qni U® U" a(s)—a(s 
fh il (43) 
where eee en | (44) 
a(s)=1-2{ |” eteo(eyq) de de; Ee (45) 


in which D, is the average number of collisions suffered by a nucleon 
in a diametrical passage through a nucleus of atomic number. A; and 
w(€,, €) de, de, is the probability that, in an inelastic nucleon-nucleon collision, 
the scattered nucleons have energies in the ratios ¢«,,¢«, to that of the incident 
nucleon. From (39) and (42) it follows that 
LOG | rere (46) 
also, from (41) and (43) one has 
J(s) ={2—a(s)—a(s—1)}{a(s)—a(s—1)}. (47) 
The mean square angular deviation from the shower axis, according to (26), 
is obtained immediately by taking the quotient of (40) and (38). The result, 
it may be observed, is independent of the variation in density of the medium, 
which is entirely taken account of by the variation of 6 with z; the denominator 
6’ does not appear in the result. 
The mean square distance from the shower axis is determined by 2,(U, @); 
however, eqn. (35) can only be solved when the solution of (34) is known. This 
equation is readily solved by Mellin transforms; one has 


b(U, 0) — = | a U-@+) exp {—h(s)0} i tig exp fh(s)A} dads ....0. (48) 
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where as, 0) = 5 an U*{x9(U, 8) —9(U, 8)} dU 


= U,* Y(s) [exp { — h(s— 1)6} —exp {- h(s)}] Sones (49) 
according to (40). 
The determination of z9(U, @) from (35) is now precisely similar : 


b(s, A) 


2(U,0)= 55 U-©+D 2 exp {—(s)9} f. Ty exP UHUS)AY Adds (50) 
where 
b(s, 0) =U," exp {—h(s)9} fe on exp {A(s\A—A(s—1)}—1] aA 2a-.s (51) 


is the Mellin transform of 5(U, 6). 


SS. bX Ph Ci (RS Wildes 

We now proceed to summarize results for the various situations of physical 
interest: in an isothermal atmosphere, and a medium of constant density, for a 
single incident primary particle of energy U, and a primary power law spectrum 
with exponent y, and for differential and integral spectra. Small letters will ~ 
relate to differential spectra, and the corresponding capitals to integral spectra. 
The suffix i will characterize the value for a single incident primary nucleon; the 
suffix p, that for the incident power law spectrum. Primed symbols will indicate 
values correct only for a medium of constant density, e.g. lead. 

The average number of particles per unit energy range at depth 6, measured 
in units of the interaction mean free path, has been denoted so far by x)(U, @); 
we shall now revert to the customary notation by replacing this quantity by 
n(U,@). Thus one has the boundary conditions 


nT, 0 =0) =U 20a). ON, (060) = eee 
UU, US. 
m(U,0=0)= 4 0, U<U, 
U,"/U", U>U, rf 
N,(U,0=0)= 4 ie Des loky ata aati (52) 


where U, is the geomagnetic cut-off energy, measured in proton mass units. From 
previous work (Messel 1951 a, b), it is known that 


1 | So +400 U\-(G+) 
n(U,9)= 5 a | (x) expt ace 


0% 8)—t00 


N,(U, 6) = at | si (7) exp { — h(s)6} ds/s 


8otic~o / U7 \—(s-+1) a 
n,(.U, 0) = mar | (7) | exp {-— BE Y 


Cw 8 —1i0 


89-400 —s 
N,(U, 8) = ml. (z) exp (2H) oat. fine (53) 
where A(s) is given by (46), (44) and (45). 
For U>U,, the last two integrals may be evaluated immediately 
ny(U, 8) =y(U,/U)’ exp {—h(y)0}U, 
N,(U,0)=(U,/U) expt hye}, Saaae = ae (54) 
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From (40) one has for the mean square angle 


1 1 1 sett / U7 \-6+1 , 
w;*(U, 0) = n(U, 0) 2ni val (x) J(s) [exp {—h(s — 1)6} 


8-10 


—exp {—/A(s)6}] ds 
F 1 13 Vl) pert Tp \ S29 
Q?(U, 6) = NU, 0) 2ni lite (z) j(s) [exp { —h(s — 1)6} 
—exp {—h(s)0}1S 


“ 1 1 1 tio / U7 \—-s+)) . 
(U8) = am a Ga], (ae)  selexp (-Ms— 9} 


| —exp {= A(6)8})] PS 
7 1 1 1 pete 7), 
Q,7(U, @) = REPUB) ad U. ilee @ J(s) [exp { —h(s — 1)6} 
d 
—exp {A }] ay mes (55) 


where sy<y+1, and j(s) is given by (47) and (45). 
For U>U,, the last two integrals become 


w,°(U, 6) = = vu, it 1)[1—exp {h(y)6 —h(y + 1)0}] (radians)? 


©,,7(U, 0) = Boat) é D Hy +1) —exp {h(y)@ —h(y + 1)0}] (radians)?. .... (56) 


To evaluate these results, we have taken y =1-1, and 
W(€1, €2) = 120€,€o(1 —€, —€.) 

a(s) =1—240{(s+2)(s+3)(s+4)(s+5)} *  ...... (57) 
which were found (Messel 1951 a,b) to give best agreement with experimental 
data for the average numbers. With these values, eqns. (56) simplify to 

w2(U50)= 2:55 (UU) Ale 2) 
Oey s+ ee) er (58) 
The results for the mean square distance from the shower axis, derived 


from (50), are as follows: 
For media of constant ee 


/2(U, 0) = Bs, i 3 tet ico Fae 2 atta 
Be N, Gave 9) 2m re a (zy oy TOES Bae. | (59) 
where J’(s, #) = esp (ED -exp Ko 


6 exp {—h(s)6} 
{h(s) —h(s — 1)} 
For the power-law spectrum 
BU= ep ra tal (a) se" Ardsley) 
FE BU O\ Oni Ue TN gL NEN 
lhe od ei ete (Al 
Ng(U, 0) 277i Ue J «,-ic0 \Ue 


—1exp{—h(s\O}.0 naa es (60) 


L,'*(U, 6) = ie J(s)J’(s, O)y ds/s(y—s+1). ...... (61) 
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For U>U,, (61) reduces to 
ly'®(U, 8) =2( UU) Yyly + 1)I"(y +1, 8) 
L,'(U, @) = 2 a+ AiG) Vee ETS See (62) 
In lead, Ds =8-9, and @ is measured in units of 160 g/cm?; then 
1, 2(U, @) =1-02(UU,)“ eo18 fe? _ 1 — 0-16 — $(0-10)?} x 10° cm?. 
L540, 6) =5:3U> € oe 1-010 SOO) x 10S cnr ee ee (63) 
The following results are for an isothermal atmosphere with surface density 55 
and surface pressure pp: 


1 1 2 ptetio (04) | yds 
itu, 9) = (Bol?o8) G7, 8) an ae (x) I(s)J(s, 8) 


89—1.0 y—sti1 

So +40 USES. d 
L53(U, ®)=(Pol908) OH Bei TA, (tT) HOM ery | 
where J(s;0) =exp {—h(s)0} eas - ae : +... 7 65) 


8(s) =h(s) —h(s— 1). 
For U>U,, and y=1-1 as before, one has 


.99)2 .99)3 
1.2(U, 6) =309 (UU, c 02 {0-204 Cae so bm? 


Ze tim, Seas! 
Bes a 0-20)" 2 (0=26)" 

WP ACUr0\ = 162 Uae {0-204 — + — - ashy Kai". 2 eee (66) 
for an isothermal atmosphere of 0° c. 

The root mean square angular deviation in the atmosphere, for an integral 
spectrum, obtained from (58), is exhibited in fig. 1. The root mean square 
lateral spread from the shower axis for an isothermal atmosphere of 0°c is 
shown in fig. 2. 

In the table below, we give the root mean square lateral spread in centimetres 
for an integral spectrum with power law 1-1 in lead. 


Depth (cm) Energy in Proton Mass Units * 
Wes U>10 U>30 
1-41 - 0-04 0-03 0-02 
2-82 0-12 0-08 0-05 
7:05 0-46 0-32 0-19 
14-1 1-26 0-89 0-52 
28-2 3-42 2:42 
42:3 6:18 4-37 
56:4 9-16 6:48 
84-6 15-63 


* Proton mass unit & 10° mev. 


§6. DISCUSSION 
The mean square angular deviation from the vertical, shown in fig. 1, is 
inversely proportional to the minimum energy of the particles considered. 
Assuming the primaries are incident in the vertical direction, the deviation 
increases rapidly down to the first 250 g/cm? of air encountered, but thereafter 
tends to an almost constant value. This is in accordance with observation 
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{Lord and Schien 1950). Actually the primaries are inclined to the vertical, 
their distribution about the vertical having been determined (Van Allen and 
Gangnes 1950). If, however, this initial distribution is inserted as a boundary 
condition to the equation which determines the mean square angular deviation, 
the result for large depths is not affected, as the initial value appears multiplied 
by a factor exp (— 0-26). 

In determining the mean square distance from the shower axis, assumed to 
be vertical, we have not taken account of the angular deviation from the vertical 
of the primary particles. If this were done, one would obtain misleadingly 
large values owing to the large lateral distances travelled by particles which 
are initially almost horizontal. On the other hand, if one wishes to know the 
lateral spread of showers initiated by particles inclined to the vertical, this is 
easily obtained by taking a shower axis inclined to the vertical at the same angle. 

The mean square distance from the shower axis, exhibited in fig. 2, is also 
inversely proportional to the minimum energy of the particles considered. 
Leaving aside for the moment the magnitude of the results, it will be seen that 
there is a rapid increase down to an atmospheric depth of about 400 g/cm?. ‘This 
is solely a result of the fact that the density is so small, which enables a particle, 
once deflected at a high altitude, to travel large distances in a lateral direction 
‘without collision. 

Top of Atmosphere 
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Fig. 1. Root mean square angular deviation of nucleons Fig. 2. Root mean square lateral spread 
with energies greater than U plotted against the of nucleons with energies greater 
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The curves are for an incident integral primary depth. ‘The middle line corresponds 
power law spectrum with exponent 1.1. to the axis of the shower. The 


curves are for an incident integral 
primary power law spectrum with 
exponent 1.1. 

Following the first collision, a particle may easily travel a distance of several 
kilometres from the shower axis before making the next impact. ‘Thereafter 
the density increases and the collisions become correspondingly more frequent. 
It can be seen from this that there will be large fluctuations in the distance 
travelled by a particle from the shower axis, particularly in the first few hundred 
grams of material traversed. ‘This will need to be taken into account in 
interpreting the experimental data. 
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Below 400 g/cm? the mean square distance from the shower axis decreases, 
owing to the loss of particles from the energy range considered. The mean 
lateral extent of showers is therefore greatest at about 8000 metres. For an 
equivalent atmosphere of constant density, no such maximum extent would occur. 
From the table it will be seen that, in lead also, no maximum spread is attained 
in the first few metres; a maximum exists, but it occurs much deeper than in a 
medium of variable density. It is therefore most important to take the variation 
of density with height into account in considering the development of a cascade 
in the atmosphere. This, however, was not done in previous work on the lateral 
spread of the soft component of the cosmic radiation. 

The actual values for the mean square distance from the shower axis are 
unexpectedly large. This does not necessarily imply that individual showers 
are diffused over a wide area. The relation between the mean square distance 
from the shower axis and what one would normally regard as the radius of the 
shower depends critically on the behaviour of the radial distribution at large 
distances. In the present instance, it is necessary to interpret the large value 
found for the mean square distance as due to a long ‘tail’. In a qualitative way 
it is easy to see how such a ‘tail’ may arise from the fluctuations, noted above, 
in which a particle can travel many kilometres away from the shower axis. It has 
already been pointed out by Messel (1951 c) that the number distribution has 
also a long ‘tail’, for a precisely similar reason: the ability of particles to traverse 
long distances in the upper atmosphere, with large fluctuations in the collision 
mean free path. ‘To discuss this phenomenon from the quantitative point 
of view, it will be necessary to determine the radial distribution function 
relating to particles in a shower; this we hope to do shortly.* 

The present results for the atmosphere have been evaluated using an integral 
power law spectrum for the primary particles with the exponent 1-1, which is 
now most generally accepted. We have also worked out results for the 
exponent 1-7, and find values at sea level which are still much larger than those 
given here. 

Experimentally, the lateral spread of extensive air showers has been 
investigated recently by Cocconi et al. (1949), and Greisen et al. (1950), up to 
a distance of 100 metres from the core; and by the Russian school (Eidus, 
Alymova et al. 1950, Vernov et al. 1950, Eidus, Blinova et al. 1950) up to 
distances of the order of akilometre. Cocconi etal. and Greisen et al. found densely 
populated cores to the showers, with the soft component decreasing more rapidly 
than the nucleon component towards the periphery. The Russian school found 
a core, consisting of both the soft and hard components, of mean radius 
100 metres ; outside this core is a region, of the order of 1 km in radius, in which 
the nucleon component predominates. In this region the particle density 
decreases as r-*® with the distance from the core, a behaviour which, if continued, 
would give an infinite mean spread. Since we find a finite value, there is 
presumably a ‘cut-off’ at a distance of many kilometres. We wish to stress at 
this point that our results do not indicate high density showers (detectable by 
ordinary counter arrangements) diffused over wide areas. Assuming our 
interpretation of a long tail then what would be observed, experimentally, are 
showers of high density diffused over distances of several hundreds of metres 
from the shower axis, and surrounded by a very wide area in which the particle 


* This has now been done. The results will appear in the Physical Review. 


Angular and Lateral Spread of the Nucleon Component Jol 


density is comparatively small and difficult to detect experimentally. It is, 
however, desirable to carry out experiments over the wide areas of low density 
which would provide satisfactory data for a quantitative comparison with the 
theory. 

The reason why Moliére’s theory, based on a physical picture which is now 
generally admitted to be erroneous, gave apparent agreement with the 
experiments of Cocconi et al. (1949) was the fact that it is actually correct for the 
central region to which these experiments were limited. 

We infer that the centre of a shower consists of a core, composed of 
descendants of the particles of very high energy produced early in the nucleon 
cascade. Among these particles there will exist neutral mesons, which through 
their decay into y-rays, generate the soft component, which predominates in 
the core. ‘The particles in the outer region are the descendants of particles with 
somewhat lower energies, among which the z-mesons are not sufficiently 
energetic to generate electrons and photons in the same proportion as found in 
the core. 

These conclusions, which were already drawn from the result of our 
calculation of the mean square angular deflection in nucleon-nucleon collisions 
(Green and Messel 1951b), have been confirmed experimentally by McCusker 
and Millar (1951) who found that the proportion of the penetrating to the soft 
component increased with distance from the core. It is reasonable to suppose 
that a satisfactory quantitative discussion of the soft component can be based 
on the foregoing considerations. 
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The Altitude Variation of Penetrating Showers 


By A. L. HODSON 
The Physical Laboratories, University of Manchester 


Communicated by P. M.S. Blackett; MS. received 6th November 1951, and in amended form 
21st April 1952 


ABSTRACT. The altitude variations of local penetrating showers and of extensive 
penetrating showers have been studied between sea level and 33500 ft. 

The attenuation length in air of the radiation producing local penetrating showers is 
found to be (129+ 2) g/em®. This value has been corrected for the effect of deviations of 
the radiation from vertical incidence on the apparatus. Local penetrating showers observed 
under a transition layer of lead 10 cm thick are about 750 times more frequent at 33 500 ft. 
and 500 times more frequent at 30000 ft. than at sea level. 


§1. INTRODUCTION 
ROADBENT and Janossy (1947) showed that showers at sea level containing 
penetrating particles are of two types: ‘local penetrating showers’ 
produced in the materials of the apparatus by single incident particles 
and ‘extensive air showers’ containing penetrating particles or radiation capable 
of producing them. Most of the local penetrating showers are produced by 
fast nucleons, and the altitude variation of these showers is a measure of the 
attenuation of fast nucleons in air. The events which discharge a penetrating 
shower detector in extensive air showers are not necessarily of this type, and it 
is necessary to separate the two types of showers if we are to measure the 
attenuation of nucleons. 'Tinlot (1948 b) and George and Jason (1950) attempted 
this separation at mountain altitudes, but in the existing measurements at higher 
altitudes (Wataghin 1947, 'Tinlot 1948 a, Walsh and Piccioni 1950) no separation 


was made. 


§2. EXPERIMENTAL ARRANGEMENTS 

The penetrating shower detector used in our second series* of measurements 
is shown in fig. 1. ‘The events recorded were sevenfold coincidences between 
the trays of counters T (split into three interlaced counter groups), M (split 
into two) and B (split into two). Each shower recorded therefore contained 
at least three particles, of which at least two could penetrate 22 cm of lead. 
This detector is similar to that used by Broadbent and Janossy but smaller. 

An unshielded counter tray Eg of total sensitive area 1190 cm?, was placed 
78 cm from the centre of tray 'T. Sevenfold coincidences (denoted by P) and 
sevenfold coincidences accompanied by a discharge in Ey (denoted by P, Eg) 
were recorded simultaneously. ‘The events in which a sevenfold P coincidence 
was not accompanied by a discharge in Eg are denoted by P—Eg. 

Broadbent and Janossy (1947) found that a 10 cm layer of lead placed over 
the P set produces a large increase in the rate of local penetrating showers, 


_* In 1947 two flights were made in which sixfold coincidences between T, M, B were recorded 
with and without discharges in En. The results were similar to those obtained with the sevenfold 
set, and the mean apparent attenuation length deduced was 117 g/cm®. These results were 
communicated by Professor P. M. S. Blackett to the Cracow Conference, 1947. 
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but little change in that of extensive penetrating showers. In our experiments 
two blocks of lead were mounted on horizontal rails so that they could be moved 
from fairly remote positions to positions directly above T. During flights the 
blocks were moved from one position to the other every four minutes. The 
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Fig. 1. The penetrating shower detector. 


difference due to the lead, in the rate of P events unaccompanied by discharges 
in Ep, is expected to be due almost entirely to local penetrating showers. The 
apparatus was placed inside a duralumin case (a bomb case) and carried in the 
bomb bay of an R.A.F. Mosquito aeroplane. 


§3. EXPERIMENTAL RESULTS 


Three flights were made in February 1948: the aircraft climbed as rapidly 
as possible and flew at constant altitude on an E—W line, the average geomagnetic 
latitude being 56° N. 

The sea-level rates used for comparison were corrected for casual coincidences 
and for triple-knock-on events (Janossy 1942). It was found that the altitude 
increase of (P, E,) events was about the same as that of all sevenfold coincidences, 
and since it seemed unlikely that local and extensive penetrating showers should 
show the same altitude variation, a cause was sought. It seemed possible that 
penetrating showers produced in the aircraft structure, and especially in the 
fuselage petrol tanks above the bomb bay, might discharge the tray Ey as well 
as the sevenfold set. Four more flights were therefore made in November 1949 
with additional unshielded counter trays in the rear compartment of the 
aircraft, remote from the P set and separated from the outer air only by a thin 
layer of plywood. 

Two of these trays E;, Ey each contained three counters each of active area 
119 cm? and one half-size counter of 60 cm?. ‘Tray Ey; contained seven large and 
three half-size counters. In order to obtain information about extensive pene- 
trating showers the discharges of all these counters were recorded separately in 
a neon lamp hodoscope, whenever they were coincident with a master pulse 
from the sevenfold P set. 

In the flights of this third series there was no petrol in the tanks above the 
bomb bay. The results were in good agreement with those of the second series. 
The large altitude increase of (P, E;,) events was confirmed, and a similar increase 
was found for the showers which discharged the distant trays E;, Ey,, Ey; as well 
as the P set. The results of the two series are combined in table 1, and the data 
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on local penetrating showers are discussed below. The results on extensive 
showers will be treated elsewhere. 

In fig. 2 the hourly rates of events, P, P—E and the transition difference 
(P—E with 10 cm lead) —(P-—E with no lead) are plotted logarithmically against 
the atmospheric depth. In each case the variation is well represented by an 
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Fig. 2. 
Top curve: P events 
Middle curve: P—E eee 10 cm lead, 


Bottom curve: ‘Transition difference. 


exponential law. ‘The best-fit attenuation lengths for the various events, 
calculated by least squares, are shown in table 2. 

Since the events (P, E) are a small and nearly constant fraction of all P events, 
the fraction of the P events due to extensive penetrating showers which fail to 
discharge the unshielded trays may also be taken as constant over the altitude 
tange with which we are concerned. ‘The altitude variation of events P—E 
may thus be taken as representing closely that of single fast nucleons. 


Table 2. Apparent Attenuation Lengths (in g/cm? of air) 


x=0 2=10 cm Pb ‘Transition 
Altitude range pao P—E 1p P—E difference 


33 500 ft.-sea level 112:04+1°5 110-4413 1185-04-11 TSS te 152 1 2055 ae 


A correction should be applied to the result because we do not measure 
the attenuation of a parallel beam of nucleons. We assume cosmic-ray 
particles are incident isotropically on the top of the atmosphere and that 
fast secondary nucleons preserve the direction of the primary. ‘The intensity 
of fast nucleons at depth x in the atmosphere at an angle @ to the vertical is given 
by I(x, 0) =I(« sec 0, 0) = A exp (—px/cos @) assuming exponential absorption. If 
the detector has an angular variation of sensitivity P(?) for penetrating showers 
formed by an incident particle making an angle @ with the vertical, the rate of 
local penetrating showers recorded at a depth x is 


m2 — px ; 
R(x) = | exp (=5) P(6) sin 6 0. 
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Since it is not possible to calculate P(@) for a given detector it was determined 
experimentally in a sea-level experiment. The bomb case containing the P set F) 
was rotated about its horizontal axis and the counting rate was determined for 
30° and 50° rotation. Since any fast nucleons at sea level must be incident — 
nearly vertically, P(@) is determined directly in this way. The angular 
sensitivity of the P set is approximately represented by P(#)=cos@ for the set 
without the movable lead blocks and by P(#)=cos? 6 for the set with the blocks. 


0 


| Curve | e-? 
\ » 2 Isotropic Detector 
NN : » 3 cos 6 Detector 
» 4 cos? 6 Detector 


» 5 cos*é Detector 


Log Intensity 


Fig. 3. Altitude response curves for various detectors. 


The angular variation of the difference (lead in—lead out) is well represented 
by P(@)=cos*@. It is reasonable on geometrical grounds that the angular 
dependence should be greater for the showers produced in the 10 cm lead blocks 
than for other showers. We now evaluate R(0) for the three cases. 


(1) No lead blocks, P(@) =cos 6. 


Writing = =2 


=$e-“(1 —a)—4a°®Ei(—a), 
where Ei(a) is the exponential integral. 
(ii) 2 =10 cm lead (blocks in place), P(#) =cos? 6. 
R(a) = ge “(2—a+a?) + 4aE(—a). 
(iit) ‘The transition difference P(@) =cos* 0. 
afl 1 eis a* 5 
R(a) =e (5 — 354+ apt im) + 5 El —a). 
The values of R(a) are plotted logarithmically against a in fig. 3, together with 
R(a) for a detector of uniform angular sensitivity. ‘The P set thus follows 
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curves 3—5 instead of curve 1, the exponential fora parallel beam of particles. (Over 
the range of depths (a=2-5 to 8) with which we are here concerned curves 3-5 
are nearly exponential.) The ‘apparent’ attenuation lengths deduced from 
the experimental curves are clearly shorter than the ‘true’ attenuation length 
for a parallel beam. Over the range considered the corrections to be applied 
are +16% for an isotropic detector, +14% for a ‘cos6’ detector, +12% for 
a ‘cos?@’ detector and +10%% for a ‘cos*@’ detector. 

On making the Sona rete corrections we find the following values for the 
true attenuation length: 


A=126 +2 g/cm? from P—E with no lead blocks. 
A=130+ 1-3 g/cm? from P—E with 10 cm lead. 
A=132+43-5 g/cm? form the transition difference. 


It will be observed that with the correction for inclined primaries the three 
values of A are in much better agreement than the uncorrected values. 

The most probable value for the attenuation length of fast nucleons in air 
is 129+2 g/cm? where the estimated error allows for statistical fluctuations and 
the uncertainty of P(@). 


§4. COMPARISON WITH THE RESULTS OF OTHER WORKERS 

While preparations were being made for the third series of flights, Tinlot 
(1948 a, b) reported measurements of penetrating showers at aeroplane altitudes, 
mountain altitudes and sea level. At first sight it might appear that our final 
result for the attenuation length of fast nucleons in air (129+2 g/cm?) is 
significantly higher than the frequently quoted value obtained by Tinlot 
(118+2 g/cm?). In deriving the latter value no distinction was made between 
extensive penetrating showers and local penetrating showers, and no correction 
was made for inclined primaries. Although it is necessary in principle to separate 
the two types of penetrating shower, it turns out that the rate of extensive 
penetrating showers is so small, and varies with altitude in such a way, that the 
lack of separation is not very serious (see table 2). ‘The difference between the 
two results is mainly due to the correction for inclined primaries. Since Tinlot’s 
detector was geometrically very similar to the one used here with X=10 cm of 
lead, his result may be compared with the apparent attenuation length applicable 
to coincidences P for & = 10 cm Pb (table 2), i.e. 118-0+ 1-1 g/cm?. The agreement 
here is very satisfactory. 

From measurements with a large sevenfold coincidence set at sea level and 
at the Jungfraujoch (3457 m) George and Jason (1950) deduce an apparent 
attenuation length for the primaries of local penetrating showers A=114 +410 
g/cm”. 

Using a much less selective type of detector, Walsh and Piccioni (1950) found 
an apparent attenuation length of 112+2g/cm*. By assuming that their detector 
was not direction sensitive they concluded that the true attenuation length was 
A=140 g/cm?. The correction which they make for the effect of inclined primaries 
seems rather high. However, although not stated explicitly, it may be that 
Walsh and Piccioni derived their value only from measurements made above 
25000 ft., since at low altitudes many knock-on events produced by single mesons 
were recorded. A higher correction would then be applicable. 
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The Calculation of Scattering Amplitudes * 


By G. J. KYNCH 
Department of Mathematical Physics, University of Birmingham t 


MS. received 26th May 1952 


ABSTRACT. In two earlier papers a differential equation was given for the asymptotic 
phases required in the two-body scattering problem. ‘The method implied a resolution of 
the wave function into eigenfunctions of the orbital or total angular momentum as the 
potential was central or non-central. This feature can be avoided and in this paper a 
integro-differential equation is given for the total scattering amplitude corresponding to 
a given incident wave. It can be applied to any scattering or ionization problem, single 
or multiple, where no particles are bound after scattering which were not bound before. 
‘The method is applicable to the Schrédinger or Dirac equations. 

If used with time-dependent equations the usual time-dependent perturbation theory 
results, owing to the apparent necessity of adhering to a definite time distinction between 


initial and final states. The same is true for the quantized field equations in the Heisenberg 
representation. 


§1. INTRODUCTION 

HE importance of scattering problems in atomic and nuclear physics has 
led to a wide selection of theoretical methods of handling them. The 
general problem is to find the probability or cross section that a system of 
particles which starts in a state A should end in a state B after a reaction or 
collision process. The state A is completely described by the states of the 
individual parts and their directions of motion and energies. 'The same 

information is needed to describe the state B. 
The wave equation gives us rather more information than this; it leads to 
the scattering amplitude f(A: B) for the process, and while this gives the cross 
section, which is proportional to |f|?, the amplitude f also contains the relative 


ot This paper forms a continuation of two papers entitled ‘ The Two-Body Scattering Problem 
with Non-Central Forces ”, published in Proc. Phys. Soc. A, 1952, pp. 83, 94. 
t Now at University College of Wales, Aberystwyth. 
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phases of the waves. This is as it should be, since this extra information is. 
needed to obtain predictions on more complicated experiments on polarization, 
double scattering, and so on. In fact, the scattering amplitude f(A: B), given 
for all A and B, can be in principle determined by scattering experiments and 
is the most that can be learned from them. This was first clearly stated by 
Heisenberg, who introduced his now-famous S-matrix theory, this matrix. 
being simply related to f in the matrix formulation of quantum mechanics. 

One very useful method of calculating the scattering amplitude is a 
perturbation theory based on an expansion in powers of an interaction constant, 
especially when the first-order term by itself is adequate. However, this is not 
true for the collisions between nuclei or nuclear particles, where the first-order 
term or Born approximation, as it is usually called, is not accurate enough even 
at high energies. The inclusion of several terms of higher order is then 
necessary to improve the accuracy sufficiently. In polarization experiments. 
for example, Born approximation often leads to the answer that there is no effect: 
at all. In problems of nuclear solutions it is more usual to make an expansion 
in terms of angular momenta and obtain an exact solution. This method is quite 
efficient if only a few eigenvalues of the angular momentum are needed. It is 
probably seen at its best in the theory of nuclear reaction, where, coupled with 
the idea of the compound nucleus, an immense range of results have been 
successfully interpreted. In the two-body problem it is cumbersome when 
tensor forces are introduced, and the three-body problem is even worse. 

In these circumstances it seemed worth while to attempt a generalization of 
a method recently published (Kynch 1952a, b, to be referred to as I and II) 
for the two-body problem.:A very brief account of this method illustrates some 
features of the more general method now being given and the reasons why it 
was chosen. We expressed the wave function for a given total angular momentum 
as a matrix 4=Yua/r where u=u,+u,S; Y is an angular function, u, and uz 
are free wave functions, a(r) is an amplitude and S(r) is the phase-matrix whose 
limit as r— 0 is effectively the f(A: B) of the problem. A suitable definition 
of a(r) leads to an equation for S: 


S(R) = if u(r) V(r)u(r) ar. 


This is an exact equation and S occurs in the integrand. Now this equation 
bears a strong resemblance to Born approximation and reduces to it if S is 
neglected on the right-hand side; also it does not contain the amplitude 
factor a(r), and can be derived without introducing that factor at all. The 
important feature of the equation, however, for its generalization turns out to 
be the fact that S(R) and u(R) have a physical meaning. Since S(0o) is the 
phase-matrix for the actual potential, then S(R) is the phase-matrix for the 
same potential cut-off outside a sphere of radius R, and u(R) is with constant S 
the free wave function outside this sphere. By an examination of Dirac’s method 
(1932) of deriving Born approximation as a perturbation theory we were led 
to a suitable equation which can be applied not only to individual angular 
momentum but to a scattering problem as a whole. 

Using the stationary state method we use an extension of Born approximation 
to calculate the change in the scattering amplitude f(A: B) due to a small change 
dV in an existing potential V. Except for constants this is 


df(A: B) = [ (—B)dVf(A) dr, 
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where (A) is the total wave function of that state of the system with interaction V 
where its constituents are incidents in state A, and can emerge after scattering 
in various ways, including B. Let #b(—B) be the total wave function describing — 
the system incident in the state (—B), i.e. the same as B but with the direction 
of motion of the various parts reversed, and scattered in various ways, one of 
them being (— A). Notice that, except for the incident waves, ¥(A) and 4(— B) 
contain only outgoing waves. Thus ¥(—B) and %*(B) are not the same, as 
*(B) contains the incident wave plus incoming waves (cf. § 2). 

This equation is not by itself sufficient to give a practical method of 
calculation. It is combined with a method of assembling the potential, piece 
by piece, so that at every stage (A) and y( — B) are calculated from the asymptotic 
forms using the wave equation either with no potential or with a simple potential. 
This means that, apart from the incident wave and these simple solutions, we 
derive equations which only contain the scattering amplitudes f(A: B). 

In view of the generality of the problem it is not to be expected that the 
resulting equations are simple, or that their solution is easy, but we can show 
by suitable examples that our equation does lead to considerable simplifications. 

At this stage in fact we try to establish three advantages of the method. First 
we are able to derive new results, e.g. equations for the total scattering amplitude 
in the two-body problem. Second, when we discuss such problems as the 
scattering of electrons by atoms, we see that we have a simple extension of 
methods used previously (Mott and Massey 1949), which helps to make clear 
the approximations made. ‘The reason for this lies in the third advantage of 
our procedure, that the functions %(A) and ¥(—B) are determined directly from 
the asymptotic forms and the scattering amplitudes. If any process is neglected 
on physical grounds, for instance that it makes little difference to the cross section, 
then a term of a known form is neglected, so that the error can be estimated 
directly. 

At the moment the field of application does not include all bound state 
problems. Where the grouping of the incoming particles into atoms or nuclei — 
is different from that of the outgoing particles, it can be used only to discuss the - 
effect of any additional interactions which come into play. / 

Equation (1) is very similar to equations given by Schwinger (1951). In 
fact, if we apply our method to time-dependent equations we obtain just his 
equations (2.14) and (2.129). For this reason we have not discussed the 
quantized fields equations at all in this paper, although most of the results for 
example on the S-matrix can be derived very quickly in the Heisenberg 


representation. In §5, non-relativistic time-dependent perturbation theory is 
mentioned briefly. 


§2. THE SCATTERING OF A PARTICLE 
(1) Generalization of Born Approximation 

The general method to be adopted is illustrated by the scattering of a particle 

by a centre of force. The wave equation is written 
(V2+ Rk?) = Vif, Lee GU (24) 

where k? is the energy and V the potential, in suitable units. In the 
method of Stationary states, to obtain the Born approximation we write 
(Ry) =exp (tk, .r) + b(A;) where exp (7k, .7) represents an incoming plane wave, 
and solve the resulting equation for %, with the condition that it only contains 


The Calculation of Scattering Amplitudes 711 


outgoing waves, i.e. for large values of 7 it can be written 
exp (ik 
th ~f(Ry 7) sable TRE i rand (2.2) 


It should be borne in mind that in this and later formulae the function f depends 
only on the direction k, of the incident wave and the direction of r, but not on 
their magnitudes. To first order in the potential we obtain the Born approximation 


(Ry: Re) = — = | exp (—2k,.r)V(r) exp (tk,.r) dv. ...... (2.3) 


Let us now construct an analogous formula when there is a small change 
in the potential. Let (k,) be the correct wave function for an incoming plane 
wave exp (7k, .r) with a potential V, and 4(—k,) be the correct solution with the 
same potential and an incoming wave exp (—7k,.7) 


W(h) ~ exp (ikyr) + f(b 7) SPE) 


( — Ry) ~exp (— ik, .7) + f(— Re :7) eR) 


First, since #(k,) and %(—k,) satisfy the same wave equation, we notice that, 
in a sphere of radius R 


0= [ [VM —ho)(hex) — W(— Aa)V9A(R,)] do 
= | (A ve) Wm) HEY) ean, 


We proceed in the same way as Dirac (1932), using the asymptotic form (2.4) 
for large value of R. In the limit R- o we find that 


fh Beare Me Vas | Ceriite: cs Key Jer Sal =, SN) Or Pam) sao (2.5) 


_In fact the non-zero terms arising from f(k,:7) are of the form 


: ; 7] exp (zkr) ‘ 0 (exp (tkr) " 
Him fb 20) 45 (exp (thy) “PED — exp (—it.r) 5 (SPE) i! R* dQ. 
= ares Thaker eS mee: (2.6) 
It is now easy to prove that an increase dV(r) in the potential results in an 
increase in the scattering amplitude 


df(k, : ky) = — = { Wah) Vik) dv.  cacans (2.7) 


Let ¢’ =4+d be the correct solution for the potential V’=V+dV. Then 
the wave equations give 


[ il’ — he) AV h,) do-= [ [V2 (— ha) (Ry) — i (— ha) VU) 1 do 
= | [V%dyp( — Aa yb(,) — di — hy) Vb(hy)] de. 
Since asymptotically dib( —k,) =df( —k,:1r) exp (tkr)/r we derive 
[ (Re) dV p(y) do = = 421 df( — hy : =) = — 47 dffhy hy). 
In the limit as dV +0 this goes into eqn. (2.7). 
(ii) The Method of Assembled Potentials 


We now establish a method which enables the calculation of f(k,:k,) for any 
potential V(r). This is done by choosing V and dV in our equations so that the 
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integral only extends over a very small region in space, and in that region (hj) 
and ¢(—k,) are determined by the simple equation (V? + k?)/ =0. 

Let o and o’ be two closed surfaces such that o’ encloses o and let us assume 
that V=Vi(r) inside o, V =0 outside a, 

V+dV=Vi(r) inside o’, V+dV =0 outside o’. 

Moreover let , be the solution for the potential V and , the solutions for 

the potential V+ dV such that y,(k,) has the form (2.4). Then eqn. (2.7) gives 


1 © 
— Ge les ths) fol ia)} = | Ye — hy) dV daha) do. 


In the limit as o’>o this gives a differential equation for f, where the right-hand 
side involves a surface integral, and the %, are only required on the surface of o. 
But on the surface of o and outside they are determined by the free-space wave 
equation and the asymptotic form. Thus the only unknown in the equation 
is f = . 

In practice it seems to be most convenient to choose as the surface o a set of 
spheres, since %, assumes a simple form. Thus if o is a sphere of radius R and 
o’ a sphere of radius R+dR we find that eqn. (2.7) becomes, in the limit dR 0, 


= he at) | fin — Valen 4D, «6 (2.8) 
the integral extending over all surface elements dS = R?dQ, dQ =sin 0 d0d¢ of 
the sphere. 

This can be transformed into an explicit equation for fp(k,:k,) using the 
solutions of the wave equation in spherical harmonics. 

A far neater form is obtained by using the integral representation of the solution 
of the wave equation 


: La nas R,: 

p(y) =exp (tk, .7) + sa] dp | d2,p° exp (ip. ue ee (2.9) 
the integral over p being taken along the p-axis from — 00 to + oo in the complex 
plane above the pole at p= —k and below the pole at p= +k. 

The resulting equation for fp is (dpt =p? dp dQ,) 
1 dfr(ky: Re) k? 1 fr(k,: p) dptV(p: ky) 
BeBe dice = OA) peat er 
BE | V (Ri: 2) dp'frl P : a) 
Fa Rp 


k2\2 Re eet "tk 
+ (3) ie) ap'V(p:p') dp tae, 


where V(p:p’) =| exp (—ip'. R)V(R) exp (ip. R) dQn.* 

‘T'wo points may be mentioned. First the V(p:p’) are not the matrix-elements 
of V in a momentum representation; the equation cannot be used therefore 
with non-local potentials. The second point is that although the last integral 
is convergent it leads to angular 8-functions if the integrations over p and p’ are 
carried out before the integrations over angles, as it is usually more convenient 
to do. 

This equation is immediately available for the two-body problem, and in 
fact any problem of the collision of two particles where only elastic scattering 
occurs. It can also be applied to multiple scattering problems. 


* See author’s correction added in proof on p. 764. 


| 
| 
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(it) Inclusion of Spin and Exchange 

With non-central forces it is necessary to include spin in our equations. 

This is quite trivial. We write 
: é exp (zkr) 
Lm(Ry) m™ exp (tk, : Nae ay Xfm’ RG : r) r Xm’ > 

where x,, form a complete set of spin-functions, e.g. for n—p scattering it is 
convenient to choose for y,, the triplet and singlet functions. The eqn. (2.10) 
is now replaced by a set of equations for the elements frm, z. ‘It is perhaps 
easier to arrange the %,, as a column matrix, so that f,,.,- is also a matrix. This 
fits in better with our general exposition. 

If the potential is of an exchange type, this character is taken into account 
in evaluating V(p:p’). It occurs nowhere else. 


(iv) A Lemma on Asymptotic Wave Functions 


Before we consider a more general system it is necessary to examine an 
essential step in the calculations which is not obvious for a general system although 
it may be so for a two-body system. ‘This is the statement that the asymptotic 
form of the wave function determines its values elsewhere and, in particular, 
at those points where dV is not zero. ‘This can be demonstrated with the aid 
aid of the following lemma: 

Lemma. If ¢ is such that in the region outside a sphere o (i) (V?+k?)¢ =0 
outside o, (ii) lim,_, ., (06/0r —tk¢) =0, (111) ¢ has assigned value ¢, on the surfaces 
g, not all zero, then, as r> 210, d—>g(0,¢)e""/r where g is not zero. The proof 
of this lemma is simple. Outside the surfaces o the solution ¢ can be expressed 
as a series rf =e'*"D_g,(0,¢)r-" where r is the polar distance from any point. 
Because of (i) the g,, satisfy the recurrence relation 


0 0 
= boas hale: 
20k sin? O(n + 1)2,4, = {sino 5, (sino, a) +5 - +n(n+ 1) ben 


Thus, if g)5=g =0, then d=0. 

It follows from this proof that there cannot be two solutions of the equation 
(V?2+k?)=0 satisfying (i) and (ii) with the same asymptotic form, but having 
different values on the surfaces o, because their difference would satisfy (i) 
and (ii) and have g=0, so that their difference would be zero. 

This result can be extended easily to more than three dimensions. 


§3. A GENERAL SCATTERING PROBLEM 

The calculation of the previous section can be extended at once to calculate 
the scattering amplitude for the process where n particles A, B, .... N with 
momentum k,, k,, .... k, come together and are scattered by a potential V 
and emerge with momenta k,’, k,’,....k,'. Since bound states are considered 
later we shall assume that none are possible. We neglect spin since its inclusion 
merely involves an extra suffix. ‘The potential V includes interactions between 
the particles and external forces if they occur, Le. 

{(Ve+ViA-+ ... EV )+ RA HVd, wee (3:1) 

where K?=k,?+...+h,,. 

The calculation is made in a 3n-dimensional space. If we assemble the 
potential in the same way as before using a (3n—1)-dimensional surface o, then 


our equation for fis __ df. (k:R’) = C| b,{ —K)Vibd) do, eer) 
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where f,(k:k')=fo(Ras Ry. ++ RniRay-+-Ry 3 do represents a volume element 
between o and o’, C is a constant depending on 7, K and the rest of the notation 
is obvious. Many-body potentials cause no difficulty in this formulation. 

The proof that %, is determined on the surface o by its asymptotic form 

b(k) ~exp {i(k,.7at-.-+Ry-Pa)}tf(k: Reet/RO OP, ....,. (3.3) 
where R is the distance from the origin in the 3”-dimensional space, follows at 
once from the lemma of the last section. 

This method could be used, for example, in collisions between three nuclear 
particles under their mutual interactions. The centre of gravity would first 
be eliminated, leaving only two coordinates 7, and 7, chosen so that the kinetic 
energy term in the Hamiltonian is (V,?+V,). The potential would be expressed 
in terms of R and suitable angular coordinates. 

Although this procedure shows the validity of the method, one objection at 
least can be raised which suggests that calculations should be made in other 
ways. Since potentials usually arise from two-body interactions, the integrations 
in (3.2) are bound to introduce a number of delta-functions corresponding to 
one, two or more particles not being scattered at all. ‘Thus we could better 
express the asymptotic form (3.3) in a form which recognizes this explicitly and 
uses the individual 7,, 7,, etc. 

To make this clear we choose an even simpler example of particles A and B 
scattered by an external potential V=V ,(r,) + V,(r,) with no interaction between 
them. It is clear that either A or B or both can be scattered, and that we write 


elk) ~ exp {ilheg ty + Rn-te)} Hf Rai 7a) ee Feo) exp (ik, 1, 


Va 
+ faler,) pee) ee 
b r re 
This is correct, of course, since the wave equation is separable and f,” =f, f,’.* 

The three-body problem is not considered further in this section for the 
following reasons. In practice only two bodies collide so that we must consider 
the problem of bound states. It is true that S-matrix theory suggests that an 
analytical continuation of the results to complex relative momenta should give 
the required answers; the only reply to this is that the behaviour of the S-matrix 
during this process is still a matter at best of controversy and at worst of complete 
ignorance. 

Three further generalizations are useful. The first concerns coulomb 
potentials. It is possible to include some potentials on the left-hand side of 
eqn. (3.1), such as a coulomb term. ‘This requires the use of the corresponding 
coulomb solutions instead of the plane waves exp(ik,.r) and outgoing waves 
exp (7kr)/r in the asymptotic form of % and the use of the wave equation with a 
coulomb term in the derivation of % on a hypersurface from the asymptotic form. 
It must be possible, and the analytical difficulties are not enormous. The second 
suggestion is that % can be obtained in an integral representation in terms of its 
asymptotic form, and although the completely general expression has not been 
found, many terms in the integral equation for the f, can certainly be expressed 
as in (2.10) and it seems likely that the rest can also. 

Finally the occurrence of identical particles causes no difficulty. If two 
particles are identical the calculation proceeds as above, and at any stage the 

* Incidentally, this can be proved using only eqn. (3.1). 


a 
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correct solution is obtained by making the solution antisymmetrical in the two 
particles, if they have spin }, or symmetrical otherwise. 


§4. BOUND STATES 

As a further illustration we consider such problems as the scattering of an 
electron by a hydrogen atom. When the nuclear motion is ignored, this is in 
effect the scattering of two interacting particles A and B by a centre of force. 
The potential energy is V=V,(r,)+V,(r,)+Va,(ra,) where r, and r, are 
distances from the centre or nucleus. The wave equation is, in suitable units 

Vb =(V24+V,2)¢=(V —E)y. aor Baan << (4.1) 
The two particles are assumed distinct, in agreement with the remarks made at 
the end of §3. . 

In the absence of the interaction V,, this equation is separable and has 
solutions of the form 4,,(A)x,,(B) where the solutions ¢,,(A) are a complete set 
of regular normalized solutions of (V?+E,,)¢m=V am and similarly y,,(B) are 
solutions of (V?+ E,,)x,=V,x,. The total energy is E=E,,(A)+£,(B).t+ Some 
of these states are assumed to be bound states. 

When there are bound states and sufficient energy one possible process is 
that B collides with A, which is assumed to be in a bound state «, and A is ejected 
leaving B in a bound state 8. In accordance with our general procedure we 
consider two solutions of eqn. (4.1) of the form 


P(1) =exp (thy .7y)b.(A) + ¥°1), 

y( —2) =exp(—2k,.7.)xB)+¥"(—2), . «2.00. (4.2) 
when the incident particle B comes in from the direction k, and A is emitted 
in the direction k,. ‘The total energy is E=hk,?+E,=k,?+ Ep. 

The scattering amplitude which interests us is contained in the asymptotic 
form 


* ex 1R aha 
| xe (B) dv, ~f(a, ky: Ra, ee pay (4.3) 
Since (1) and %( —2) satisfy the same wave equation, 


[£924 — 21) — H( — 2) V24(1)} do, de, =0. 
Let +/'(1) =4(1) + df(1) satisfy the wave equation with a potential V’=V+dV. 
‘Then, as in deriving eqn. (2.7), 
fu —2) dVip(1) du, dv, = fivy( — 2) dib(1) —( —2)V? dib(1)} dv, dv,,  ... (4.4) 
neglecting products of first-order corrections. 


Because of eqn. (4.3) we expand the wave functions in terms of the complete 
set of functions x,,(B) for particle BT 


a1) =XdF,,x,,(B), (—2) =G,,x,, (B), 
where dF, and G,, are functions of r,. Using the wave equation satisfied by the 
x,(B), the right-hand side of (4.4) becomes 


* Y * 7 AI 
bagd | (eee Va Z V, 79 En)GuXn ad Gok Aes as V ba» E) CHIR A dv, dv. 
mn 

+ To avoid difficulties in the subsequent discussion due to the asymptotic form of the wave 
functions when Vg and Vp are coulomb potentials -we assume that they are cut off outside a radius r. 
The rest of the potentials is added to the term Vgy. ‘This separation of the potentials into parts 
can be related to the screening effect of one particle A when bound to the nucleus. , 

t The expansion of ¥(—2) in terms of x,* corresponds to the change of direction of motion 
of B in ¥(—2), for an unbound state, compared with the emerging waves from 4(1), when x», is 
.complex. For bound states we can choose real xn. 
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The orthogonal properties of the x,, reduce this to a single sum 
¥ | [dF (VG) - (V2 dF ,)] 40 
n 


which is essentially the same as eqn. (2.6) of the previous section. Using Green’s 
theorem and the fact that all the dF,,, and G,, represent outgoing waves except 
G, which contains in addition the incoming wave exp(—1k,.r,), it reduces to 
a single term in dF;. Now 


i ot ae 
| dib*(1) yp"(B) do, ~ df(on, Ry: pees Aa bs (4.5) 
Consequently after some simplification we obtain 
4 df(a, ky: Ra, B) = — J o(—2)dV f(1)dv,dv,. ...... (4.6) 
By an expansion of dj(—2) and (1) in terms of the ¢,(A) we can show that 
Af(a; Re? Ris B) = GCak, Bee Bel eee (4.7) 


so that if the two f’s are equal when the interaction is first included then they 
remain equal subsequently. 


It is quite clear that this procedure applies whatever the transition to be 
discussed. ‘. 

These equations are now transformed so as to be suitable for calculations 
or approximations. As an illustration let us suppose that E is greater than the 
ground state energies E,(A) and £,(B) of A and B but less than the first excited 
states. Then the asymptotic form of (1) and %(—2) are 


ee +hup (Ryir yore) re) xa(B), 
Y( — 2) ~ x@(B) {exp (—tky.1a) +hpp( — Raia) ce a2 sete 


exp c a 


+ fd — Beat?) PO) 4.4), 


the abbreviated notation being obvious. There are three equations for f,,, fgg 
and fug(2y' ha) =foal — Ra? —hy)- 

The prescription for building up the potential rests on the fact that f,5=fg, =0 
when V=V,+V,; we therefore start at this point, which is quite convenient 
as the wave equation is separable, and f,,, and fgg can be calculated. The potential 
due to the interaction of A and B is now built up so that at every stage the values 
of #(1) and ¥(—2) can be determined on a certain surface from the asymptotic 
form and the known solutions for the separable equation with the potential 
VitV,. 

When the energy of the incident particle B becomes sufficiently large to make 
possible the use of the excited states of A (or B) then additional terms must be 
added to the wave functions, and additional equations appear for the new scattering 
amplitudes. 

These calculations are very similar to those previously made for the scattering 
of electrons by atoms, and there is no difficulty as long as V ,, isa repulsive potential, 
or, when it is attractive, not so large that particles A and B can form bound states. 
When we increase it up to the magnitude where a bound state first appears, a 
singularity appears in the direction r,=7, which cannot now be ignored. Our 
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present method breaks down and a discussion is now necessary of the nature 
of this singularity along the lines suggested in the previous paragraph. Such 
difficulties arise whenever complex particles A and B collide to give particles 
C, D, ... and the potential whose effect on the scattering is under consideration 
happens to be responsible for the binding of any of these particles. 


§5. THE DIRAC EQUATION AND TIME-DEPENDENT EQUATIONS 
Once the method has been established there is little point in giving details 
of all the applications. For example, the results given in I and II can all be 
derived directly by the present method after resolving into separate angular 
momenta, or alternatively by expanding the scattering amplitude f(k,:k,) of 
eqn. (2.10). 
In this section, therefore, we give a few results on other wave equations. 


(i) Dirac Equation 
In II we considered an electron in an external field, using equations resolved 
into angular momenta. We now establish a more general result. Given the 
Dirac equation (the notation used is that of II and §2) 


fe cA) Spe Ei View yao eee (5.1) 


and its conjugate complex for %* with a slightly different potential we can easily 
obtain 


| Wha)” aVU(hy) ao 
1 

~ 2 

where n is a unit radial vector in the solid angle dQ and the arrows over V indicate 


which way they operate. The second term tends to zero as Roo. If the 
‘wave functions have the asymptotic forms 


Why) =U, exp (ths) +f :7) PE), 


| a (Re) (pak + E) .[n( V _v) —io.n x (V + v)] dip(k,)r? dQ, 


Wa) * = ug exp (— thy. 7) +f*(—hy 7) SPD, 


where u and f are spinors, the first term reduces so that 
— Zarit," (py +E) df : Ba) =? [ y( eg)" dVib(Fy) do 
(= —2ndf"(—Ry: —Ry)(pse + E)uy), «ae (5.2) 
which reduces easily to eqn. (2.6) in the non-relativistic limit. ‘There is also an 
equation analogous to (2.10). By using the simplest forms for mw, and u,*, and 
p3 diagonal we obtain equations for the corresponding component of the 


amplitude f. Moreover, with the usual restrictions on the potentials it can 
easily be verified that the components of f are related, as suggested in II. 


(ii) Time-Dependent Equations 
Time-dependent perturbation theory considers the development of the 


‘wave function with time, due to a given perturbation. ‘This exactly corresponds 
to our procedure for building up potentials. In fact, starting from the Schrodinger 


equation iO 
(ve+ims) EN (5.3) 
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and proceeding as in §2, we arrive at the equation 


[w'(2) avy) dr =i f [vow ao| 


The integration is taken over a four-dimensional volume between times ei and 4 
and dV =V’—V say. If (2) =(1) and we expand in terms of the eigenfunctions 
of V then we obtain the usual time-dependent expansion in the form of an integral — 
equation. 
Alternatively, let (1) =, at t=, describe the incident state, and (2)* = $2" 
describe one of the states in which the particle can emerge. Then we can easily — 


verify that 
ae { | boul)=de | = | | FGA ao| ies (5.5) 


= 


t=ts 
t= 


1 


The function /*(2) develops in time reversed in the same way as (1) does in 
time, cf. (5.3). Finally from eqn. (5.4) 


im dao = J Way dVil \dodt,. eee eee (5.6) 
where Qj, is the probability amplitude of a transition from state ¢, to the state do. 
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ABSTRACT. Previous treatments of nuclear evaporation have dealt with average values,. 
neglecting fluctuations. In the present paper the fluctuations are considered in detail. The 
distributions of the number of neutrons and of the total energy associated with evaporation 
stars of a given size and the size distribution of stars of given total energy, are evaluated 
explicitly. 

Knowledge of the fluctuations is useful for interpretation of data concerning small stars, 
such as those produced by z-meson capture or by artificially accelerated particles. The 
fluctuations are not very important in the large stars and do not modify any of the results of 
the author’s previous mean value calculations. 

The analysis of the mean evaporation process has also been put into a simplified form, 
from which numerical results can easily be obtained. 'The method has been used to compare 
experimental data with various theoretical energy—temperature relationships ; the Fermi 
gas law previously used by the author seems to be the most satisfactory connection between 
the high-cnergy star data and the observed level density at low energies. 
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§1. INTRODUCTION 


ETAILED investigations of Harding, Lattimore and Perkins (1949), 
Fujimoto and Yamaguchi (1949, 1950) and Le Couteur (1950, referred 
to as I and II), have shown that a very large class of nuclear disintegration 

stars can be described by evaporation theory. The different treatments have 
recently been compared by Yamaguchi (1950) and reviewed by Rochester and 
Rosser (1951) and Camerini, Lock and Perkins (1952). Since the calculations 
of I were completed further independent experiments of Perkins, Bernardini 
and Rochat, described in the second review, have confirmed the accuracy of the 
expected ratios of proton, deuteron, triton and «-particle emission. Also 
observations of Hodgson (1951, 1952) of the frequency of emission of heavy 
splinters, hammer tracks and beryllium nuclei are in good agreement with 
predictions of Le Couteur from the formulae of I. 

The experimental material used in I referred mainly to stars in silver and 
bromine with an average size of about 10 prongs and the formulae of I must be 
a good representation of the nuclear temperature at the average excitation, about 
200 Mev, of such stars. 

It is rarely possible to measure the energies of all the charged particles emitted 
from a particular star and the neutrons are not observed, therefore the total 
energy released must be deduced from the star size by a combination of theoretical 
and statistical arguments. According to fig. 21 of Camerini ez al., the scale of 
average total excitation energy established by fig. 2 of I for silver and bromine 
is in agreement with experiments. With stars of a given size there must, however, 
be associated a whole distribution of excitation energies about this average 
value. 

In this paper the fluctuations of the evaporation process are considered in 
detail. The distribution of star sizes arising from evaporation of a given total 
energy and the distribution of total energy associated with evaporation stars of 
a given size are worked out explicitly. One might expect the fluctuations to 
increase the width of the energy spectrum of particles emitted from stars of a 
given size over that calculated by the mean value methods of I; the increase 
however turns out to be negligible. Some types of fluctuation were very briefly 
considered by Fujimoto and Yamaguchi (1950) but previous treatments of 
nuclear evaporation have dealt mainly with the mean emission process; this 
procedure is justified by the results of this paper which show that the fluctuations 
are quite small. 

In many experiments nuclear disintegrations are produced by artificially 
accelerated particles or by cosmic rays and only the size distribution is observed. 
To interpret these experiments one needs to infer the associated energy 
distribution which can usually be related to the primary collision process 
(e.g. Bernardini, Cortoni and Manfredini 1950, Barton, George and Jason 
1951). For the smaller stars such as are produced artificially or by 7-meson 
capture, knowledge of the fluctuations is very desirable. 

One would like to go further and analyse the experimental data to yield some 
information about the variation of nuclear temperature with energy. ‘This 
requires consideration of the. change of disintegration characteristics, such as 
energy spectra and relative particle frequencies, with the total energy or rather 
with the star size, and knowledge of the fluctuations is essential. 

Consideration of the fluctuations entailed the development (§ 4) of a simplified 
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treatment of the mean evaporation process which turns out to be quite accurate 


and easy to use. In the final sections of this paper, this method is used to compare — 
the experimental energy spectra from moderately large stars with the predictions 


of various temperature laws. 

To fix the law of temperature variation one must cover the largest possible 
range of energies, therefore in parallel with this work the disintegrations produced ~ 
by low excitation energies of about 20 Mev are being studied by Miss Lang, in 
this department. 

The combination of material may give a reasonable indication of the nuclear 
specific heat. 


§2. FORMULATION OF THE FLUCTUATION PROBLEM 

The fluctuation problem presents no formal difficulties. Let 
PA Zp, Op Ag Ze SEU, = OB Cd. ae (2.1) 
be the probability that a nucleus b of mass A, charge Z,, with excitation energy Uj, 
(often denoted by Bor b, (Us) disintegrates to leave a residual nucleus c with 
excitation energy in the range dU,, by emission of a fragment of type x (neutron, 
a-particle etc.) of mass A, = A,—A,. If U, is below the threshold «, for nucleon | 
emission from b the probability vanishes, if U,, is above «, the nucleus must 
disintegrate somehow and 


U), 
ay | pCR Cydthaie i) "eles eee (2.2) 
0 
In the notation of I, the total probability of emission of x from B is 
P,=P(B:c)=[p(B:C)dU, and %,P,=1.. ...... (2.3) 


The binding energy Q,,, of x in b is defined by I eqn. (50), and the kinetic energy 
of the emitted particle is 

Tp=U), = USO 6 Oa eee (2.4) 
The elementary probabilities, which may be derived from the theory of 
Weisskopf (1937), are given explicitly in I, §§3.2 and 3.3. The probability of 
emission of x from B is of the form I (56) 


Pl = PUB: 6) =y,po0 p= soe eee (2.5) 
Here y,, defined by I(51), is a statistical factor proportional to the mass and 
number of spin states of x, and p, and 7, are the level density and temperature 
of the residual nucleus c at excitation energy R, = U,—Q,,—V,’, where V,’ is 
the effective potential barrier for emission of x from b. The factor p is chosen 
to make XP,=1. The distribution of 7, and U, is 


Ja Hie) aT, = exp [— (fs a Ve) tele Fate 


1 
= p. PB nO) (oe a ee eT Se (2.6) 
With this distribution the mean value and variance* of 7’, are 
T= 2+ Ve and o(T)=To (Fatt ea nee (2.7) 


In the notation of 1 (66), (67) the mean excitation energy of the residual nuclei c 
formed from b is U;,— H,, which is determined by (2.4) and (2.7) as 


H, = "Ph, =, PLO aa) ee (2.8) 


* Mean square deviation, always denoted by o?. 
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Now let P(O:f)=P(A,,Z,, U,: Ay, Z;) be the probability that an excited 
nucleus O finally reaches stability as f with residual excitation energy below the 
threshold e- for emission from f. Conventionally P(O:f)=8(0,f) for U,<e,; 
for energies above the threshold the probabilities are determined by the integral 
equation 


U 

P(O:f)= | BOTA) CULPOlet) atk) Lf bee (2.9) 
av0O 

The iterative solution is 


P(O:f)= ["(O:F)du,+3 || pO: A) dU, pA; F) dU, 


res | || p(0:A) dU, p(A:B) dU,p(B:F)dU,+..., Gee (2.10) 
a,b 


with U,>U,>Ug....>U; and «->U,;20. The successive terms represent 
contributions from processes in which one, two, three emissions take place. 
The number of terms is limited by the initial energy U,, since each process 
requires a threshold energy. The normalization 


pl (Ok) te eS Ge i A (2514) 
follows immediately from (2.2). 
The probability distribution of numbers of each type of emitted fragment 
may be described by the moment generating function 


RAE 2 P(Oe 7, 5... .)EXP(TA SY ae..44)) 0 Geaee- (Zeal) 


where P(O:7,s....) is the probability that 7 particles of type x, s of type y etc. 
have been evaporated in the disintegration of the initial nucleus O and X, Y.... 
are parameters. WM is determined by the integral equation 


WOK) =>, [HO rA) dU, exp(X,,) MALX, Yn. ie... (2.13) 


where by X,, is understood the parameter associated with the particle (x, y....) 
emitted in the nuclear transition o—a. 

The method may be further generalized to yield the joint probability 
distribution of the kinetic energies of the emitted particles. Let X,Y,. 
represent arbitrary functions of the kinetic energies of the corresponding 
particles x, y, ...., so that’ in (2.13) X,, is a function of U,—U,—Q,,, then 
M becomes a characteristic functional (Bartlett and Kendall 1951) defining 
the energy distribution. ‘This is easily understood by consideration of the 
special case in which all parameters vanish except one which is constant in a 
certain kinetic energy interval and zero outside it; then M gives the distribution 
of the number of particles emitted in this energy range. The substitution in 
(2.13) of 

PEM eS ea FCC ema | geen. (2.14) 
leads to a moment generating function M(O:A) which defines the probability 
distribution of the sum of mth powers of the kinetic energies 7, of the particles 
of type x emitted in the disintegration of O. 


§3. A SIMPLIFIED PROBLEM WITH ENERGY FLUCTUATION ONLY 


To clarify the content of the formalism, let us consider a simplified problem 
in which only one type of particle can be emitted so that the nuclei may be 
specified by mass and energy only. In this case (2.9) or (2.13) are alternative 


PROC. PHYS. SOC. LXV, ¢—A 47 


] 


722 K. 7. Le Couteur 
ways of working out the total number of particles emitted. The generating 
function for the number of particles emitted into the kinetic energy range «, f is 
determined by 

M(O:X)= f[p(O: A)dU,e*7M(A:X), a aes (3:1) 
where X vanishes outside the range /=U,-Q-—«>U,>U,—Q-—f=m, 
determined by (2.4), in which it is constant. The mean 7(O) and mean square 
7>(O) of the number of particles emitted into this range by disintegration of O 
are determined from the coefficients of X and X? as 


7(O) = [.°9(0:a) aUs(A)+ | “WOTAldU | oe ae (3.2) 
72(O) = [.(0:) 0,71) + { p(O: A)dU,{1427(A)}. ...... (3.3) 


The integral equations may be solved approximately by expanding everything 
about the point U=w say. Define 


Uo Uo 
U,= | °PO:A)U,dUy a= J, p(O:A)(U,—w)" aU, 


yet ig Ap be ie es ae gee vate a 
v= | p(O:A)(U,—w)" aU, 
and then (3.2) becomes 
UU) Sie ae = = 
( Omm em + 31( et) — Hal ae TROD Ute oD ot ( C ) 
and (3.3) becomes 
= dr ie Mae 
(U,- Us) Ta +31(U,— wv)? — pe] Fs +....=%42 (Fn Ze +...) : 
ee ae (3.6) 


where 7, 7? and their derivatives are evaluated for U=w. 'To simplify matters, 
choose w to make the terms in d?7/dw* vanish. This requires 


(U,—w)2=p,=(U,—U,)?+0%U,) or w=4h(U,+U,)—}4o%(U,)/(U.— U,), 


i anand (3.7) 
where o7(U,) is the variance of U, given U,. The equations are then 
ath ae d— - Qy\id? a has 
ete eye aap ele Se Pe 
eee an {72 — (r)?} (14 a) + Fae ets (3.8) 


where, as in (2.8), H,=U,—U,. For the mean and variance of the total number 
r of particles emitted, these equations reduce to 

ar 1 dane ee dr o%(U, oy, 

At) aie, Lae ite comes 
dw ie dw ’ “do HS dw? H, 
For integration it suffices to take w= U,—4Q, which is precise at low excitation 
energies when the nuclear temperature is negligible so that U,=U,—Q, o?=0. 
Thus one has, for U>2Q 


3 a rU+42 dU, } U+3@ dU, 
CS ee a aga ( E., ‘ddpash alte deca (3.10) 
a a rU+4Q dra U; U+10dU, e ‘ 
cA) - Moma | ae elgg ee 
~ 5Q/2 0 = 5Q/2 oO 
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and (3.8) gives an equation like (3.10) for the energy spectrum of the particles. 
The constants of integration are settled by reference to the exact solution 
which is obviously r=0 for U<Q, r=1 for 20>U>Q. Therefore (3.10) was 


adjusted to give the mean value 7(Q)=4. Similarly the variance r? — (r)? vanishes 
if U<2Q and must increase rapidly to } at the value of U, slightly above 2Q, 
for which r=1 and r=2 are equally probable; therefore (3.11) was adjusted 
to give the mean value { at U=20. 

In fig. 1 these approximations are compared with the exact solutions derived 
by numerical integration of (3.2) and (3.3) for the Weisskopf spectrum (2.6), 
with 7,=[4(U,—Q)]}*2 and Q=10mev. In the range U>2Q, for which a 
solution is required, the approximation (3.10) is of remarkable accuracy. In 
(3.11) neglect of the negative term in d?7/dw? leads to a negligible overestimate 
of the variance of r at high energies; at energies below 3Q the solution is too 
simple to follow the bumps in the curve of variance. The accuracy of these 
approximations is also good in the less realistic case of constant temperature 
which can be solved in terms of incomplete [’-functions. 


4 


) 10 - 20 30 40 50 
U7 (MeV) 


Fig. 1. Exact and approximate (broken curves) values of the mean 7 and variance o*(r) of the 
number 7 of particles evaporated from a nucleus with initial excitation energy U. Binding 
energy Q=10 mey and nuclear temperature (U/6)1/* Mev are assumed. 


Equation (3.10) for the mean emission is in agreement with the formulae 
used in the detailed mean value calculations of Le Couteur (1). 

Equation (3.11) for the variance of r is easily understood. Emission of one 
particle leads to a spread +o(U,) in the residual excitation energy, which 
implies a spread +0(U,)/H, in the number of further particles emitted; so 
the variance of r is made up of contributions o?/H,” from each particle emitted. 

For the evaporation distribution (2.6), (2.7) the variance of r becomes 

yo 27 dr i 
is 6 wee uh Pasieoat <2, 
because the temperature is usually much less than the binding energy Q. ‘The 
variance of r is much less than that, 7, of a Poisson distribution with the same 
mean and for this reason the total number of particles emitted is a good indication 
of the total energy evaporated. The physical reason is that only the kinetic 
energy, but not the binding energy, removed by each particle is subject to 
fluctuation. 

In experiments usually only the emission of charged particles is observed, 

but, since (§6) the neutron emission is fixed within very narrow limits by that 
47-2 
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of the charged particles, the number of prongs remains a good indication of 
total energy. 

This analysis is easily extended to the case where different types x, y.... of 
particles can be evaporated, if the emission probabilities are considered to be — 
functions of excitation energy only. The mean and variance of the number n, 
of particles of type x are given by 

dn jdw=P iy 3h Oe ee (3.12) 

in agreement with 1(68) and 

do*(n,) = P,(1 — Px) dr + P,2drd?/H? +2P2dr(1—h,/H), ...... (3.13) 
where P,, and H are given by (2.5) and (2.8) and &? is the variance of the residual 
energy after emission of any type of particle. Equation (3.13) shows that the 
variance of m, is the sum of contributions from (i) a binomial fluctuation in the 
relative number of different types of particles emitted, (ii) a proportion of the 
fluctuation in the total number 7 of particles emitted, (iii) a small term depending 
on the difference between the mean energies removed by different particles. 

The differential energy spectrum is given by similar modification of (3.8). 
Alternatively one can use the method of (2.14) for the moments of the spectrum, 
in particular for the Weisskopf distribution (2.6) 


ZT, = [(27,+ V»’)PxdUy|Hp, 
DTP = [[(2r¢+ Ve’? +272] Py dUy/Hy, 


§4. EXPLICIT SOLUTION OF THE SIMPLIFIED EVAPORATION 
PROBLEM 
The formulae of §3 can be integrated analytically if eqns. (3.10), (3.11) are 
combined with a simple energy temperature relationship. Consider the power 
law 


Lota 2 Ne eee (4.1) 


The Fermi-gas nuclear model corresponding to 6=2 was used in I and is 
approximately correct; here 6=2 or 3 will be considered. 

If Q represents the sum of binding energy and potential barrier, (2.8) gives 
H,=Q+27(U,—Q) and (3.10) becomes 
nua [Nat a (TaN aA tl 
Q+27(U,—Q) 0 Q+2r Tm Jo On 2r° 
where 7,,=7(U—}Q) is the maximum temperature involved. This may be 
expressed in terms of dimensionless parameters, 


EO 


t=27/O-. (sD Oe nn ees ee (4.3) 
2 COR sat U,—3sQ2 
i \ aaa Saas oe 0 2 
as OE | Prine) i. Oe (4.4) 
5 5 
where X~ (1 sSeet gage) ere (4.5) 


may be expressed exactly in terms of logs and is graphed in fig. 2. The average 
nuclear energy removed by each particle is Q/X, and so the average temperature is 


7=40(X1-1)~ hoe for 6=2 


au —37,2/0 for 6=3 
which is also shown in fig. 2. 47m ~ 407m 
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Similarly (3.11) gives the variance of r as 
soni baie sil U-430 t? 
%(r) = parereprren lati ed ieegr ese Vl ws ete 
op (r)=$+ {, (O+27)3 st O. Ore 
ee ya Bf ttt BB 
= 542), 7240) 542543 See ee ae 
is drawn in fig. 2, Also, from (3.14) the expected value of the sum of squares 
of the emitted kinetic energies is given by 


2dU U-30O 60m (2 Mat r2 
See bomen. ho Be ea ce ee 
(T-V’) Se 0 ™ 54D TAS (4.9) 
6=2, Uae ct? 6=3, Vat? 


Z 


05 


0 04 1-0 1-4 2:0 0 0:4 1-0 1-4 20 
§=2Tn/Q S=2Tn/@Q 


Fig. 2. Particles, with binding energy O are evaporated from a nucleus with initial excitation energy 
U and temperature T,,=7(U—4Q). The number emitted has mean 7=X(U—4Q)/Q and. 
variance o°(r)=4+ Y(U—4Q)r7,,7/O* and the variance of their kinetic energy spectrum is: 
o(T)=Tp’” Z. The figure shows X, Y, Z and the average temperature 7 as functions of the 
dimensionless parameter s=2T,,/Q. 


The variance of the kinetic energy is therefore 


o(T)=672—(27)2=ty2Z, nee (4.10) 
where 
a Zou te, Ps 27 She # 
Z = (5 — 135 *) for 6=2 and Z=(5 — me) for 6=3 ..(4.11) 


and is drawn in fig. 2. 

The curves in fig. 2 are convenient for practical use. Strictly the constant 
of proportionality in (4.1) varies slowly with the nuclear mass and so changes 
during the evaporation. ‘This can be allowed for adequately by using the average 
mass to evaluate 7,,. 

The results can be applied to stars where different types of particle are emitted 
if the mean H is approximated as in (4.2) with an appropriate Q. ‘The results 
concerning energy spectra remain valid if in (3.14) the variation of P, during the 
evaporation can be neglected. If these approximations are invalid, no question 
of principle arises, but use of the graphs must be replaced by numerical integration 
of the equations of §2 and § 3, as was done in I. 
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The approximations give o°(T) =117?/4 for }=2 and 12r°/5 for 8=3. Thus 
the width of the integral energy spectrum exceeds that of the Weisskopf distribution 
(2.6) corresponding to the same mean energy and temperature by a factor (11/8)'” 
for 8 =2 and (6/5)? for 8=3: the broadening is greater for the square law than 
for the cube because the former leads to a more rapid cooling of the nucleus 
than does the latter. Therefore, as pointed out by Bagge (1941), experimental 
energy spectra should not be interpreted directly in terms of a Weisskopf 
distribution at constant temperature, though this has often been done. 

The integral spectrum can be represented by the gamma distribution 


P(T)dT =exp(—(T—V)/s*\(T-V)dT/T()r*, (4.12) 


which has mean T=V'+/r* and variance /7**, if the values 7* =117,,/12, 
1=16/11 for §=2 and 7* =97,,/10, /=5/3 for ’=3 are chosen to secure agreement 
with (4.6) and (4.10). The elementary Weisskopf spectrum (2.6) corresponds 
to [=Z, 

If the total energy U is varied, the two temperature laws considered lead to 
significantly different variations of the mean energy 7 and spectrum width o(T), 
the latter is approximately proportional to U'”? or r¥?._ This is considered, in 
relation to the available evidence, in § 11. 


4.1. Variable Potential Barrier 
There is strong evidence for a reduction of the potential barrier at high 
nuclear temperatures, as was assumed in I. This has an appreciable effect on 
the energy spectrum of the emitted particles as can be seen from the examples 
given in I, fig.5. For simplicity suppose that V falls off linearly with temperature, 
so that for protons say 
= Veh eee (4.13) 
and the parameters d, e of I are well represented by k=}. In this case the mean 
kinetic energy of the protons is 
LHVi~t (ZR oy (4.14) 
which can be evaluated easily from (4.6) or fig. 2. The change in V also slightly 
narrows the energy spectra. Analytically 


(T—V,)? =(2—k)?7? + 272 = (6-4 + 2)? 
by (3.14) which may be calculated as (4.10) to yield 
o%(T) =(1— $+ GR) {ty 2 + 4(7)?} — (2-7) 
=Ty LZ — ($k— th) 722 — 2a), ee (4.15) 
which is easily evaluated from fig. 2. Since k +4 the reduction of o?(7) is small, 
being approximately 


0:33" foro —=2, 


(3k— tk? oes eretete ere All 
MBS I crue ee, sae 


Comparison of this result with (4.14) shows that the width of the energy 
spectrum should be a better indication of nuclear temperature than is the mean 
energy, because V’ is not accurately known. 


§5. NEUTRON-PROTON FLUCTUATIONS 


In most experiments on nuclear evaporation the neutron emission is not 
observed and must be calculated as was done in §3 of I. In I, $3.3 it was explained 
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that the equations of emission I (65) contain controlling terms which tend to 
keep the neutron excess A= N-—Z of the evaporating nucleus close to the 
equilibrium value A,=N,—Z, of a stable nucleus of mass N+Z; thus the 
disintegrating nucleus tends to follow the Heisenberg valley. This effect does 
not depend on the particular energy-temperature relationship assumed in I. 
Consider the more general case 


I LG WE TNT al a ta a G31) 
where y defined for 8=2 by 1(38), is a small coefficient approximately equal 
to 1-5; then the entropy is 


AMS /77\e4 Aue § /U\e-ne 
<. ef i ae ved td ech 
A(1+ sai) =} (1+2 7) al) See ae, (5.2) 
aS 1 as 1 cay 0s tL 2yA Oo 
d i —— —— => CU — ape ee 
ee oll” (ae 3-1 Ar es), SE Re FG oe 


Now P,/P,=exp{S(A—1, A—1, R,)— S(A-—1, A+1, R,)} from (2.5) and the 
relative emission probabilities become, as in I, § 3.3, 


P, ee ee ple A 
P, = exp eae Vo (0—-1)> = | BO POs i tm aiange teres, (5.4) 
P 1 TOGA A \eaetae. >? 2yA\ R 
Beef ots SHAN rr -07(te 8) SD 
Rte (5.4) 
P 1 3-28(A —A, 4yA\ R 
p =exp{_| -38+ AEA) sy, ey a(S ty (2+) |}. 
eae (5.4¢) 
ee se ers : eck: (5.44) 
P, 0-8 ae A, 2yA\ R 
FE = exp{- =| 244 a dy Vy -V-@-1)2(34 2°) 5]h. 
we (5.4) 


In each equation 7 is the mean of the temperatures and R the mean of the maximum 
excitation energies of the two alternative residual nuclei. £, defined by I (44), is 
about 23mev. With 6=2 these equations are identical with I (65) in which 
6=A/A and 7=2,/R/A. As shown in I, fig. 1, the term in A— A, tend to make 
the final value of N—VJN, insensitive to details of the energy-temperature 
dependence and to the initial conditions of the disintegration. 

It will now be proved that fluctuations of the neutron excess from its 
calculated mean value are strongly damped by the controlling terms. Simple 
formulae will be derived for the final mean value (5.21) and variance (5.24) of 
(A-A,). 

The calculation is much simplified by neglect of the fluctuation of the energy 
carried away by the different particles (which was considered in §§3 and 4); 
all particles are assumed to remove precisely the mean energy H defined by (2.8). 
Then egn. (2.10) simplifies to 

PO) = S(O 4) EUAe be PCC T fd eye tes (O25) 
a, D.. 


with (=U,—H,, U,=U,—H,, ete. and U;<Q, 
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since the evaporation proceeds in steps from excitation energy U, to U, to — 


U,.... Let P(r, A) denote the probability of a neutron excess A at stage 7, then 
P(r+1, A) = P(r, A+1)P,+ P(r, A—1)P, 4+ P(r, A) Pat P(r, A+1)P, 


cP hd \ Pp VE bere tale tee saa (5.6) 
and thus the moments of A satisfy 


D,AP(r +1, A) =EP(r, A)[(A—1)(P, + P;) + (A+ 1)(Pp +P.) + A(Pat P,)], 


or A.5, =A, = SPO, PPee.— 2 Pon eee (5.7) 
and 
A432 =A,2—22AP(r, A)(P,+ Py— Py — P,) + P(r, A)(P,+ Py + Py + Pi) 
. (5.8) 


and therefore the development of the variance is given by 
{R2— (A)}, 4 = {A?— (A)},— (Arua —A,)®— 28(A — A) P(r, A)(P, + P,— P,— Pr) 


+ 2UP(7, A) P,+ Py + Py +P). cae cela (5.9) 
These equations involve 


Lhe a ee ie oC Pi=a OF Beste (5.10) 
say. ‘The approximation 
P,—P,+P,— P,=a tanh 26(A—A*)/7A 
where A*=A,—{V, —(8—1)4yAR/A3A/4B i... .. (5.11) 
which was used in I, §3.5 is exact at low temperatures and within 5% at hie 


temperatures. Then at the high temperatures which prevail until the last stage 
of the evaporation, the equations reduce to 


A,w—A,=—a tanh 26(A,— A*)7 A eee (5.12) 
Or412(A) —0,2(A) =a—(A,,, —A,)?—«, a sech* =] aE ae Nt ee (5.13) 


Equation (5.12) is equivalent to I (70) or I (73) for @=A/A which was solved 
in I. To test the validity of the approximations, the fundamental equations (5.6) 
were solved numerically in the special case 


H,=(10+2r,) Mev, +=(4U)!” or «/A=2 Mev and 48/A=1 Mev ...... (5.14) 


with initial nucleus {?Br at excitation energy 112 Mev. The small «-particle and 
deuteron emission was neglected. It is convenient to reckon A from the initial 
value A,=9 and the calculation yields the values in table 1. 


Table 1 

Stage 0 1 2 3 a 5 6 a 
U (Mey) oe 94 77 60 45 30 17 5 
TA/2B SO) fe? 6:3 05) 4°6 3°4 2:0 — 
eG —2 —=2:3 0 =X) = SoH] SoS one —4:1 
A— AG exact 0  -—0:250 —0-:522 —0-887 —1:18 —1:55 —2-01 —2-62 

— Ap approx (5.12) 0 —0-245 —0-522 —0-840 —1:23 -—1-64 -—2-13 —2-814 
o*(A) exact 0 0:94 1:63 1:98 2°30 2°38 2:23 1:78 
o*(A) approx (5.13) 0 0-94 1-60 2:05 DD 2°30 $597, 1-46 
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The numerical values show that A,,,—A, is always small so fair approxima- 
tions to (5.12) and (5.13) are 


hes Ba 
7 (A-A*)=-a ayik Ac ae 


by 


as) s (1- = o*(A)) sech® = (@- A*) 


with solutions 


(A—A*),=(A—A*), exp (- i, ar) s [Sa rexp (— ip =f adr) 


TA d 
ENRON (5.16) 
and 
PATS [" a sech® == ae 7 (A-A*) dr exp (- ie neg sech® =F ai 5 (h-AN) a r), 
ricoh (5.17) 
where o and f denote initial and final values. As in I (§3.5), one can write 
dr=—dUjH=—8A(rfe)  drf/H, aa (5.18) 
Te2B _ 20pa_ ([7\°"* _ 2BaS 
so that iia ne (:) aon ort EA (5.19) 


by (5.2). A mean value of the slowly varying quantity a/H is to be taken. 
Further integration yields 


eee a pee aay Aa 28 Bay 
[Far exp | - Ceara (=) |-Isa%> dy exp (Sj i) 


= 3B alae Pe oe (5 =, i) soa (5.20) 


for large U,. The values of H and a should be taken at + =e[(8 — 1)H;,/258]"0-P 
where the integrand is a maximum. For high initial energy, S>4A say, the 
first term in (5.16) is negligible and the final value of the mean neutron excess 


tends to 
dA*| «A | #(H/Ba) for 6=2 
dr |2Ba | 4(7H/3Ba)!!2 for 8 =3 
where A-* =(A,— AV’/48), and H is evaluated at 
7 =e(H,/48) for 8=2 or tr =e(H,/38)* ford=3. wae (5.22) 
This result is in agreement with I (76), though more simply expressed. 

To calculate the variance o,2(A), one may replace A—A* in (5.17) by its 
final value since the sech? differs from unity only at the end of the disintegration 
when the temperature is low. Below some temperature 7’ the sech? is negligibly 
small and the approximation 


 43Bar?2 Ala Gee aah a\* 1" 
2 Ae dr sech a, (A = A Ve = cant = 1 wii (:) ] 
SOM ety oN nce ols (3,23) 


where 7’ =48(A—A*),/A for 6=2 or 4-48(A —A*),/A for 6=3, may be derived 


(A—A*),= 


, 
if 
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by evaluation of |’? dy sech*1/y and leads to 


45 cA/3—1 H\"eD (8 \ Ar? 
ae Wess tas =) lala 45° Ba) at SS) 
eee HAs ee 


46 8Ba “He 
6B [7'\2 1A (7H Ar’ 5.24 
-[1+5(2) lag) - 7s for d= oe eteNousts ( 5 ) 
with H and a calculated at 
7 =e(H,/88) for 6=2 or T= e(H6B)" for 03. ~ 7 eee (5.25) 


The accuracy of these SHAE formulae is adequate. In the oe G14) 
|0A*/dr| =0-3, a=1, <A/8B? = 3, and so (5.21) gives (A—A*),=0°3 x 3 x 13-3 =1-33 
or A—A,=—4141:33=—2:77, in excellent agreement with “tbe accurate 
results of table 1. Similarly for o?(A) one finds 7’ =1-33 Mev from (5.23) and 
(5.24) gives o?=(141-6) x}x 4x 11-7—6 x 1-:77/10 =2-6 — 1-06 =1-5, which is 
also rather accurate. 


§6. THE NEUTRON EMISSION 
The total number of prongs (charged particles) emitted is 
No=Nek gt he ei”, eee (6.1) 
and the total number of particles emitted is 


r=N,+M,. OS eee (6.2) 
The final value of the neutron excess differs from its equilibrium value by 
pe 
My + Mye— My — Met aa < (My +My, + 2ng + 3n,+ 30, + 4n,) 
=(A—A,)-=(A—A,)e# o,(A), % etelelele (6.3) 


which was calculated in §5. If, for simplicity of writing, the numerical values 
OA /0A =OZ™ or GN 072 FS 

used in I are assumed, eqn. (6.3) fixes the neutron emission as 

n= 1:50, +, + 0-5ng— 0-50) + lng 1-254, = o( de eee (6.5) 


In many experiments the charged particles are not distinguished from each 
other, and then (6.5) must be expressed in terms of the total number of prongs N, 
and the corresponding expected numbers 7, of protons, etc. (which for silver 
bromide may be taken from I, fig. 2). The resulting expression for the neutron 


emission of an N.-prong star is 
My, = 1-5N,—0-5%,— ig — 27%, + 0-57, + 1-25(A,—A), 


+ {1:5602(A) + 4%, +7%+4m,+4}myr, (6.6) 
the standard deviation Rete obtained by adding to the variance 1-56 o,2(A) in (6.5) 
the contribution from the statistical fluctuations found in (3.13) of the relative 


numbers of «-particles, etc. amongst the N, prongs. For silver bromide this 
formula for the standard deviation gives approximately 


a(n, given N,) = 2-81? for N, <2, {2:8 +0-7(N,—2)}1? for N,>2, 


if the numerical results of §5 and of I, fig. 2 are used. ‘ 
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§7. THE PRONG DISTRIBUTION 
The prong distribution of an evaporation star with given initial energy U 
can be derived from §§4 and 6. If H,, and H, are the average energies removed 
from the nucleus by neutrons and charged particles respectively, the fluctuations 
of the number of particles emitted are subject to 


H.N,+Hot,j=U-4O=H(N,+%,)=HF (7.1) 

which, in conjunction with (6.6) gives 
Neale L25(At A) lic, 71 6 eave Sie (7.2) 
o°(N,) =[02(r) + (Hy /H)*(1-560,2(A) + 4m, +7ig +40, +42, 60... (7.3) 
where c=2-5—(0-5%,+M%y+2n,—0°5%,)/N, eevee (7.4) 


and of course H,/H=1—4V’/H, where V’ is the mean potential barrier for 
protons. 

Equation (7.2) for the mean number of prongs has the merit of simplicity 
and should be useful for the analysis of experiments, in which case the value of c 
can be taken directly from the experiments. Numerical integration as in I 
should yield more accurate values for the mean emission but is much more 
laborious. 

The main result of the calculation is the explicit formula (7.3) for the 
variance of the number of prongs in terms of the known mean emission and of 
the variance of 7, which was calculated in § 4 with sufficient accuracy for the present 
purpose. 

The main results of this paper are simple to use. As an example, table 2 
shows the full calculation for silver bromide stars, assuming the mean of 
the energy—temperature relationships d, e of I, §4. ‘Thus «e=12-4 Mev or 


4 =(12-4U/A)!?=(U/6-7)"" if the average value 83 for A is used. Also 
H=Q+2r is used with O=10 Mev. 
Tiable.2* 

(parameters d, e) 
U (mev) 50 100 200 400 600 900 
Tm (Mev) Dike eee STi 760. 9:40. 41-6 eqn. (4.3) 
Cie 0-488 0-732 1-066 iheey 1:88 232 
y 3°4 6°45 Tie, 20°5 28-0 36°8 fig. 2 
c. 25 2-3 222, Pye Pei 2 eqn. (7.4) 
Ne 0-4 1:8 4.3 8. MCT y983) igg-4 eqn. (7.2) 
o°(r) Go77 i4e 1-09 «24 3-8 6-1 fig. 2 
o? (mp, given N,) 2:8 2:8 4-4 7:5? 40 13-6 eqn. (6.6) 
o” (Nz, given U) 0-4 0-52 1:0, 2:0 2°8 4-7 eqn. (7.3) 
o?(T, given U) (Mev?) 7 16 33 66 100 151 fig. 2 


§8. SPECTRUM OF EXCITATION ENERGIES FOR GIVEN STAR SIZE N, 
_ The spread o(n,) of neutron emission implies a spread of amount 
H,o(n,, given N,)=Qo/X in the total excitation energy U; the approximation 
is a slight overestimate of H,. Also the thermal fluctuations cause a variance 


* For U< 400 mev these values of Ne agree perfectly with those plotted in I, fig. 2, which assumes 
the same parameters. At higher energies there is a small error in I, attributed to the use in this 
-energy range of an approximate path of integration in the (N, Z) plane. 
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ro(T) in the total kinetic energy of the emitted particles. Since different 
particles have different binding energies Q and barriers V’, the spectrum (2.6) 
implies that with N, fixed the statistical fluctuation in the relative numbers of 
protons, «-particles, etc. gives a further contribution 


(Op+ Vo’ —Q0)*%ip + (Qu+ Ve —O,)*H, +... =No(O,+ Ve) 


to the variance of U, where QO, is the mean value of Q+V’ for the charged 
particles. This contribution is much smaller than the others. Thus the 
variance of U given N, is 

o*(U, given N,) =(H,)?o(n,, given N,)+707(T)+N,c((Q,+V,).  ...+-- (8.1) 

It follows that 7,, is distributed about the mean 7,, =7(U —4Q) with standard 
deviation o(7,,) =o0U) dr/dU. 

If, as is usual, measurements from a large sample of stars with given N, are 
combined together to derive say a proton energy. spectrum, this will be wider 
than that calculated in §4 for the mean energy U and temperature 7, 
corresponding to N,. In the sample the average value of X(7'— V’)? is, according 


to (4.9), proportional to Tne =(7,,)?+07(7,,) and so the fluctuation of U increases 
the variance of the kinetic energy of the protons in the sample to 


o(T, given N,)=0%(T, U=U)+ sor given: Nol—2.4as8 (8.2) 


5 being the exponent introduced in (4.1) and (5.1). 
However, the values in table 3 show that the effect is very small. 


Table 3 
Nuclear constants are as in table 2, from which the first two lines are taken. 
N, (prongs) 0-4 1:8 4-3 8-8 123 17-4 
U (mev) 50 100 200 400 600 900 
(H,)?o?(my) (Mev?) 478 605 1220 2270 4505 8100 
ro (T) (Mev?) 24 101 382 1350 2800 5550 
N,o*(Qe+ Vo’) (Mev’) 1 4 30 135 240 430 
o? (U, given N,) (Mev?) 503 710 1622 3755 7545 14080 
a (U, given N,) (Mev) 22-4 26°6 40-4 61-3 86:9 118 
o (7m, given N,) (Mev) 0°55 0-49 0:53 0:58 0-68 0-76 


a(7) is so small that the difference (8.2) between o (7, given U) and! 
o(T, given N,) is always negligible; the former, which may be taken directly 
from fig. 2 or table 2 is easier to work with. So far experiments have always been. 
analysed by comparing experimental spectra of given N, with theoretical spectra 
of given U; the above numerical results show that this procedure is quite 
legitimate. 


The values of o(U, given N,) should be useful for the interpretation of 


experiments in which star size distributions are observed. Note that the 
uncertainty in neutron emission is the largest contribution to o(U). According 
to the results of §6 this contribution would be much smaller if the different: 
charged particles were identified and in this case for stars of more than eight 
prongs the temperature fluctuations would give the main contribution to o(U).. 
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8.1. Note on very large Cosmic-Ray Stars 


The previous calculations have assumed that the evaporation starts from a 
nucleus, normally stable, which has been highly excited so that initially N— Z=A 
has the equilibrium value A,. In large cosmic-ray stars this may not be the case. 

Suppose that , charged mesons are emitted from the stable target nucleus 
by the primary collision. Positive and negative mesons are equally likely, so the 
charge removed from the nucleus has average 0 and standard deviation 1/7,. 
Since N +Z is unchanged this implies a standard deviation 2,/n, of A. 

A number R of fast recoil nucleons with energy above the limit of the 
evaporation spectrum are also emitted. Of these N,=RdZ,/0A will be protons, 
observed as ‘grey’ tracks, and R— N,=RON,/0A will be neutrons. This process 
contributes 2, to the variance of A. ; 

Therefore after emission of mesons and fast nucleons in the primary collision 
processes, there is left a highly excited nucleus with neutron excess A close to 
the equilibrium value A, appropriate to its mass, but with standard deviation 
o(A) = V/ (4n, 1 2N,)- 

In calculating the neutron evaporation as in §6, the variance o,?(A) must be 
added to o,?(A), the variance of the final neutron excess. 

Experimental statistics suggest the following average numbers of grey and 
shower tracks for a given number JN, of evaporation prongs: 


N, (<5) 8 13 
N,(average) rat) 3 8 
N,(average) 0 1 3 
o,2(A) 0 10 28 


which inserted in table 3 would lead to the values in table 4. 


Table s 
N. 0-4 1°8 4-3 8:8 5 ss) 17:4 
o(U, given N,) (Mev) DAP 26-6. 240-4) | 90° —-- 167 «260 
o(Tm, given N,) (Mev) O55. 0:49: 054 | 0-85 1-31 | 1-68 


As in §8 these small values of o(7,,) do not give a significant contribution to 
o(T, given N,). 


§9. NOTE ON 7-MESON INDUCED DISINTEGRATIONS 

The disintegration of silver and bromine nuclei by capture of slow negative 
z-mesons has been studied experimentally by Menon, Muirhead and Rochat (1950) 
who find N,=1-1, o(N)=1-1 for the mean and standard deviation of the number 
of prongs. According to table 2, this corresponds to a mean excitation energy 
of 75Mev. This is probably an underestimate because of the peculiar method 
of nuclear excitation, which is discussed by Menon et al. ‘The 7-meson is probably 
absorbed in a primary process like 7» + P+P+>N+P ora +P+N-+N+WN and 
it is likely that one fast neutron is emitted before thermal equilibrium is reached 
and only the balance of energy remains.to excite the nucleus which reaches 
stability by evaporation. On this assumption, the excited nucleus starts with a 
neutron surplus of about 4, since its charge was reduced by one by capture of the 
meson and 0N,/0Z,=1:5. Since the total neutron evaporation is fixed by the 
final value of the neutron excess, the neutron emission must exceed that 
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calculated in §6 by this amount. ‘This requires a correspondingly greater 
excitation energy, so that the mean value may be estimated as approximately 
80 Mev. 

The variance 1-2 of the observed prong distribution is much greater than 
the value 0:5* given by table 2 for the mean excitation energy. A spread 
of excitation energy must be assumed with standard deviation o(U) such that 
{o(U)ON,/0U}2 +0-5 =1-2, which gives o(U)=30 Mev so that the distribution of ’ 
excitation energy is U=80+30mev. This wide spread of U is quite in accordance 
with what would be expected from the method of nuclear excitation, for the 
fast primary particles may sometimes escape immediately from the nucleus and 
sometimes traverse a nuclear diameter with large loss of energy which appears 
as nuclear excitation. 


§10. RESULTS FOR THE CUBE LAW U=47%/e 


Numerical results have also been calculated for 6=3, that is for 
Us Ate =U A = (10.1) 

which Miss Lang’s results suggest as a reasonable connection between nuclear 
temperatures at U=20 Mev and the temperatures found in large stars. Two 
numerical values were considered : 

f  ¢=8-2Mev, giving 7=(67U/A)#*, 

g « =6°8 Mev, giving 7 =(46U/A)". 
The values f were chosen to give the same temperature at U=200 Mev as the 
square law Uocr? which was used in I and in tables 2 and 3. As explained in 
the introduction, the choice ensures reasonable agreement with most of the 
available data from cosmic-ray stars. ‘The values g were considered for 
comparison. 


With these temperature laws and f, dA*/dr as in § 5, eqns. (5.21) and (5.24) give 
f (A—A*),=2-25, (A—A,);= —1-75, o,%(A) =3-9 
g  (A—A*),=1-8, (A—A,)p=—2:2, of(A)=3. 
The final variance of o(A) is higher than was calculated with the square law 
because of the higher nuclear temperatures during the last stages of evaporation. 


For comparison with the square law (tables 2, 3), the numerical values for 
law (g) with H=Q+27, O=10Mev are given in table 5. 


Table 5. (Parameters g) 


U (Mev) 50 100 200 400 600 900 
Tn (Mev) 2:9 3-7 4-7 6-0 6-9 7-94 
Z 3-1 6-2 14-5) 620:9) 2 29-7" aaa 
N. 0:3 1:7 4-2 9-0 13-000 “20 

a? (r) 0-27»... 0:51. bey 407) ea ot 3-8 5-7 
o? (Np, given N,) 4:7 4:7 6-2 9-6 12 Scemmelris 
o? (N,, given U) 0-7 0-9 1-4 2:5 3-5 5-4 
o? (T, given U) (Mev)? 11 18:3: 2p 47 62 82 

o (U, given N,) (Mev) 31-6 34:9) 45:0 "66-40 780-7 Mion 
o (tm, given N,) (Mev) 0-61 0:43)" 0*36» 0:33" 50:32 =no-at 


The values of c used were the same as in table 2. The differences between the 
square and cube laws arise because at low excitation energies the cube law gives 


- * The numerical integration of § 5 applies directly to this case so the estimate of variance should ~ 
e accurate. 


. 
| 
/ 
; 
| 
| 
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the higher temperature and at high excitation energies the square law gives the 
higher temperature. As in §8 the increase of o (T, given N,) over o (T, given U) 
is negligible. 


$11. COMPARISON WITH OBSERVED SPECTRA 

Experimental energy spectra of protons emitted from silver and bromine 
stars of several sizes are given by Harding et al. (1949) in their fig. 1*. The 
interpretation is not quite straightforward; one has to estimate where the 
evaporation spectrum ends and the tail of ‘knock-on’ or Heisenberg protons 
begins. ‘Two methods of analysis have been used: (a) The mean T and mean 
square deviation o7(7) were calculated for the proton tracks of kinetic energy 
T<V +37, and the number of evaporated particles was taken to be the number 
of prongs in this energy range, for large stars about 85%, of the total number 
recorded. This is thought to be an underestimate of T and o7(T) because the 
theoretical spectrum (4.12) predicts that 5°% of the particles are evaporated with 
a kinetic energy above this limit. The experimental points are plotted as — in 
fig. 3. (6) All proton tracks with T<V +47, were analysed as in (a). ‘These 
experimental points are plotted as + in fig. 3 and may be overestimates of T and 
o*°(T), because in the evaporation spectrum (4.12) 98%, of the particles have 
lower kinetic energy and so some ‘knock-on’ particles must be included amongst 
those analysed. For small stars of 1 or 2 prongs formed by 7 capture, the energy 
spectrum of Menon et al. (1950, fig. 3) is available. This has been analysed as 
above. 


o*(T) (Mev?) 


ee) 6 Nude ak 
Nc 


Fig. 3. Upper + and lower — experimental estimates of the mean T and mean square deviation 
o°(T) of the kinetic energy of protons evaporated from silver and bromine stars with Ne 
evaporation prongs. The curves are theoretical estimates for the parameters d, e with 
U« 7? and f, g with U« 7%. 


Theoretical values of T and o?(7) for the three sets of parameters, d, e (square 
law) and f, g (cube law), considered in this paper are plotted in fig. 3. The 


* The total excitation energies there quoted are too low: revised values are given in figs. 19, 21 
of the review of Camerini et al. (1952). 
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potential barrier is the same as in d, e of I, §4, but as explained in §4.1 above in 
the present analysis it is more easily represented by V’ = V;’ — $7 and so the mean 
proton kinetic energy was taken as T=V’+2r=4+57/3mev. The slight 
reduction (about 6% for §=2, 4°%, for 8=3) of o°(T) given by (4.15) for k=4 
was ignored; it would be partly compensated by the increase given by (8.2). 

Because of the scatter of the experimental points a very large range of star 
sizes must be considered if any conclusion is to be drawn. It must be remembered 
that the curves of mean energy but not of variance can be slightly displaced by 
making different assumptions about the potential barrier. All the curves lie 
close together and to the experimental points for stars of less than 12 prongs. 
The square law predicts a somewhat greater variance o7(7) for large stars than 
is found experimentally for stars of 15-22 prongs and, if taken seriously, this 
would be evidence in favour of the cube law or a similar law. However, the 
discrepancy, which is hardly significant, may arise because in such large stars 
thermal equilibrium is not set up until the energy has been reduced by emission 
of many particles. 


§ 12. DISCUSSION 

The main sources of fluctuation in nuclear evaporation have been analysed 
and evaluated numerically; the results for o(,, given N,), o(N,, given U) and — 
(aU, given N.) should be of some practical use. 

In all the problems considered the fluctuations have little effect on the mean 
number of particles emitted or on their energy spectra. ‘This justifies the 
procedure of I in which the very complicated equations of the actual physical 
problem were integrated for mean values only, ignoring fluctuations. It is out 
of the question to solve the cascade equations of § 2 exactly for the actual physical 
problem. 

It is known that mean value integrations are not adequate to solve the problem 
of the electron—photon cascade in cosmic-ray theory for which a refined treatment 
of fluctuations (Janossy 1950) is required. The reason is that the process is 
multiplicative so that any excess of electrons or photons produced by fluctuation 
at the start of the cascade gives birth to further particles and the fluctuations tend 
to grow as in other population problems (Kendall 1949). ‘There is no such effect 
in the evaporation problem, where the equations of motion are such as to damp 
out fluctuations, as has been discussed in detail for the neutron—proton 
fluctuation. A similar, but very much smaller, effect exists for the fluctuation 
in the total number of particles emitted. 

It has never been established how far temperature laws other than the 
Fermi gas, Ux7r*, can explain cosmic-ray stars. The results of §11 suggest that 
a law U=Ar’c! with 8 between 2 and 3, and « chosen as in f of $10 to fit the 
known temperature at 200Mev, would do, although only the square law has 
actually been shown to be consistent with the observed emission frequencies of 
a-particles, deuterons, fragments, etc. The different possibilities lead to very 
different average spacings D between nuclear levels at low energies. For 4 =93, 
U=8mey the results are 


d,e mee « = 12:4 Mev Dis: 6 6y, 
b= 7/3 « = 10-7 Mev D= 240 ev 
6= 5/2 e= 9-4Mev D=1300ev 


6=3 e= 6:8Mev D=8700ev 


i WS ek es bY 
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and according to measurements of Selove (1951) the first is very reasonable. 
On this evidence the Fermi gas law is the best adapted to span large ranges of 
energy. Its failure to follow the variation of D with A at very low excitation 
energies is probably due to the complexities of nuclear shell structures. 


Note added in Proof. Some predictions of this paper can now be checked 
against experimental results. 

(1) V. C. Tongiorgi and D. Edwards (private communication), using the 
Cornell Synchrotron as a meson source, studied the neutrons emitted when 
negative 7-mesons are captured in Pb, Sn, Al, C. If these results are inter- 
polated to mass 100, which is a mean for AgBr, one finds n,=5:5+0°5, 
o(m,)=1-9 for the mean and standard deviation of the number of neutrons 
emitted with energy below 10 Mev. ‘The considerations of §9 indicated a 
distribution w= 80 + 30 Mev of the thermal excitation energy following z~ capture, 
which leads to a neutron evaporation with average 4:5 and standard deviation 


o(My) = {(28 x 4:5)? +0-4}1/2 = 1-9, 


which agree surprisingly well with the observations. 

(2) Figure 1 of Bernardini, Booth and Lindenbaum (1952, Phys. Rev., 
85, 828) shows the distribution of numbers of black prongs in stars induced by 
400 Mev protons. ‘The number of events decreases steadily from 3 to 6 prong 
stars, but has a tail extending to 8 prongs. It seems likely that the steep part 
of the curve represents a distribution of excitation energy resulting from the 
primary collisions, and that the tail represents the fluctuation o(N,) of the 


charged particle evaporation. Table 2 gives o(N,)=1:2 for N,=7, which is 
consistent with this interpretation of the distribution. 
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ABSTRACT. The production of deuterons in high-energy collisions of neutrons with 
nuclei has been explained in terms of a direct ‘ pick-up ’ process by Chew and Goldberger 
and by Heidmann, the reaction being of the form nucleus A+ neutron-> nucleus B + deuteron. 
Here it is shown that an alternative reaction, in which the incident neutron interacts with 
a nucleon in the target nucleus, which in turn ‘picks up’ a second nucleon to form a 
deuteron, is more important at incident energies greater than approximately 200 Mev 
(nucleus A-+ neutron neutron-+-nucleus B+ deuteron). 

The total cross section for deuteron production is calculated numerically for energies 
of the incident neutron in the range 100-300 Mev and comparison is made with the results 
of the direct ‘ pick-up ’ theory. 


§1. INTRODUCTION 
MONG the products of high-energy collisions between nucleons and nuclei 
a large percentage of deuterons (of energies of the order of 50-100 Mey) 
have been observed. Detailed experiments on the deuterons produced 
by 90 Mev neutrons incident of carbon, copper and lead (first reported by 
Brueckner and Powell 1949) have been made by York (1949). Serber and his 
collaborators (Serber 1947, Serber, Fernback and Taylor 1949, see also 
Goldberger 1948) have shown that high-energy Collisions take place in two 
stages. During the first stage the incident nucleon traverses the nucleus, which 
at these energies is partially ‘transparent’, in a time of the order of 10° sec, 
interacting with and possibly ejecting a small number of the target nucleons; 
the nucleus is left in an excited state and its subsequent break-up, occurring a 
relatively long time after the primary events, may be described by an 
‘evaporation’ or ‘compound nucleus’ model. The deuterons observed in these 
experiments are produced at high energies mainly in the direction of incidence, 
during the first stage of the reaction. ‘The experimental results have been 
explained, in the main, in terms of a direct ‘pick-up’ process (Chew and 
Goldberger 1950, Heidmann 1950), whereby an incident neutron associates with 
a proton of a suitable momentum in the target nucleus to form a deuteron, the 
interaction of the incident neutron with other nucleons in the nucleus being 
neglected (the incident energy being sufficiently high). The process is thus, in 
some sense, the converse of the ‘stripping’ of deuterons to form neutrons. 

At even higher energies (from 500-1000 Mev) the primary products of 
cosmic-ray ‘stars’ (that is before ‘evaporation’ takes place) also include a large 
number of deuterons. At first sight it is difficult to account for these deuterons, 
as the direct ‘ pick-up’ cross section although large at 100 Mev falls off with energy 
as E~® (Heidmann 1950). However, at least two alternatives exist, either or 
both of which may be responsible for the production of these deuterons: 
(a) incident neutrons may react with neutrons in the nucleus producing deuterons 
and m-mesons (the rcaction p+ p—+d ++ has been observed); (b) deuterons may 
be formed as the result of a second-order process in which a nucleon of relatively 
small momentum (produced by the collision of the incident neutron with a 
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nucleon in the target nucleus) ‘ picks up’ a second nucleon in the target nucleus, 
to form a deuteron. This is the process described below, and will be referred 
to as the indirect ‘ pick-up’ process. 

This indirect process (neutron + nucleus Adeuteron+ neutron + nucleus B) 
will favour the production of deuterons of lower energies than those produced 
by direct ‘pick-up’, and it is possible that this accounts for the result of Chew 
and Goldberger (1950) that at 100 Mev, rather more low-energy deuterons are 
produced than is consistent with the direct pick-up energy. 

To facilitate comparison, the calculation has been made using, as far as 
possible, the same constants and methods similar to those of Heidmann, the main 
points of interest being the relative magnitude of the cross sections at 100 Mev 
and their variation with the energy of the incident neutron. 


§2. GENERAL THEORY 

It is supposed that the energy of the incident neutron 0 is such that it 
interacts, effectively, with only one nucleus1 in the target nucleus, the 
interaction of the nucleon 1 with another 2 in the target nucleus leading to the 
formation of a deuteron in the final state. To calculate the cross section of this 
reaction Born’s approximation is used, with V,, the interaction between 
nucleons 0 and 1 as perturbation energy. 

The target and residual nuclei are described by the Fermi model, but 
(following Heidmann 1950) the nucleons 1 and 2 are described by a wave 
function, generated by their mutual interaction V,,, which differs from a plane 
wave. The interaction Vj, is not introduced explicitly. 

The cross section for ejection of deuterons, wave vector K,, from nucleus A 
when nucleons 1 and 2 are in a state described by wave vectors k,, k, is 

2a 


a(Ky, ky) dK, = 7 — | dk,9(E— B)e(ba» Ky) dK, 
* wt I(ko, k,, k, | | K,, k,) |? ‘a I(Ko, k,, k, | se K,, k,) [?} 


where vy is the velocity of the incoming neutron, k, the wave vector of the 
scattered neutron, F;, EF; are the energies of final and initial states, p(k,, K,) the 
density of final states, and the wave functions are normalized to a box of volume V. 

R® is the matrix element for the transition when nucleon 1 is a proton and 
R"™ is the matrix element when 1 isaneutron. wis the number of neutron—proton 
pairs in the nucleus which are in even triplet states; then w=3(A—Z)Z. (It is 
assumed that the probability of the reaction occurring is small if 1 and 2 are in 
a singlet state and further that the exchange properties of the nuclear interaction 
are of the Serber type.) In eqn. (1) the use of the incoherent superposition 
{|R"?|?+ |R°™|?} rather than the coherent expression {|R"? + R™ |?} follows from 
the initial assumption that the energy of the incident neutron 0 is such that it 
interacts, effectively, with only one nucleon 1 in the nucleus—either a neutron 
or a proton— the subsequent interaction of nucleon 1 with some other nucleon 2 
in the target nucleus forming the deuteron. ‘lhe coherent expression 
{|R™? + R™|?} would occur only for a true three-body collision. 

To obtain the cross section irrespective of the initial states of 1 and 2 in the 
nucleus A, the expression (1) is multiplied by the probability of finding 1 and 2 
in a state described by k,, k, (denoted by P(k,, k,)dk,dk,) and integrated over 
all k,, k, such that |k,|, |k,|<L, L being the radius of the Fermi momentum 
sphere. 

48—2, 
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To evaluate the matrix elements A”? (or the similar matrix elements R™) — 
an approximation suggested by Chew (1948) is used. The matrix elements are — 


computed, using Born’s approximation, for an imaginary situation, in which 


particles 0, 1 and 2 are unlike and have no spin. ‘The matrix elements for the 


reaction then prove to be the product of two terms, one of which, G, is a function 
of the momentum lost by the incident neutron (k,)—k,) and the difference 
(k,—k,), only. The second term F (k)—k,) turns out to be identical with the 
matrix element for a free collision (with the same assumptions) of particles 0 and 1 
having the relative momentum (k,—k,). #(k )—k,) contains the potential Vo,, 
while G does not. The function F is expressed in terms of the differential cross 
section for free scattering of particle 0 by particle 1, for which ae 
experimental n—p scattering data are substituted. 

This procedure avoids the explicit introduction of the imperfectly known 
potential V,,. Assumptions about the form of potential have to be made when 
choosing the deuteron ground state wave function and for the wave function 
describing nucleons 1 and 2, but the exact form of these wave functions does not 
influence the result critically. 

If the momentum of the scattered or ejected particles lies within the occupied 


Fermi sphere of the residual nucleus, the cross section will be much reduced by ~ 


the exclusion principle. As the method used does not allow for exchange effects 
between nucleons (except that the scattering of 1 by 0 is described by 
experimental data and hence includes exchange effects correctly), the final 
integrations over the momenta of the scattered and ejected particles have been 
carried out excluding the occupied Fermi sphere. 


To this approximation R™= f ‘f*V),?V;,d7 where ‘’; the wave function — 


for the unperturbed system is 
Vi = Vb exp ki. ree = ee eee (2) 

#%, being the wave function of the target nucleus A, ‘Yy the wave function 
describing the final state of the system is 

Y= Vo Pdbgypp(r — re) exp (tk, .1r9) x V-¥? exp (cK, .(r,+1,)/2} ...... (3) 
by being the wave function of the residual nucleus B and ¢p that of the ejected 
deuteron. 

Following Heidmann (1950), the nuclei A, B will be described by a Fermi 
plane wave model, with the exception of nucleons 1 and 2. Then if » is the 
volume of nucleus A 


a1... . A)=y?7A3.... Ady (tyre) exp {27K,.(r,+1r,)/2} ...... (4) 
where +,,(;— Wg) is the wave function describing the relative motion of 1 and 2 
in the nucleus and differs from a plane wave by the action of the potential 
between 1 and 2, k,, K; being defined by k;=4(k,—k,), K;=3(k,+k,). If 
nucleus B and the remainder of nucleus A are described by a plane wave model 
then the integration over the coordinates of nucleons 3....A gives unity. 


[ $03... Als) ge nea ale Seen (5) 


With x=ryo—ry, y=r,—r, and z=4(r,+r,) one obtains for R™? 


RP = V8 H | exp (IK, )V™x) dx x | a (ibo"(y) exp (IK y/2) dy 


nucleus 


x | exp {-1(K, -2K,—K,).z} dz 
/ nucleus 


where K, =k, —k,. 
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The integration over z, the position vector of the centre of mass of 1 and 2, 
is limited by the nuclear boundary and gives therefore a factor 7 and 
2K,=K, — K, hence 

RP? = V-82.G(Kg, k,) F™(Ke) 


where G(Kg, Ks) =772 | “can, HAY DM(Y) exp (Ky. y/2) dy 
F™(K,) = | V(x) exp (iKy.x)dxX, eae (7) 


Since p(A,,k,)=(27)*V? and P(ks, K)dk,dK, =(9/27?L) dk, dK, (changing 
from k,,k, to k3,K; and normalizing so that fP(k;,K;)dk,dK,=1 the 
integrations being over all k;, K; such that |k,|, |k,|<Z, the cross section 
o(k,, K,) may be written finally as 


o( ky, K,) =(9/6477) (zo/fivyL®) | {|F™ (Ke) |? + [F™(Ke) [?} 
x G(k3, K,)5(E, — E,)8(2K;—K,—K,)dK,dk,dKy. —...... (8) 


(i) The Function G (ks, Ky) 
If a Yukawa potential is assumed, a very good approximation to the deuteron 
wave function y%p, which occurs in G is the Hulthén expression 


bp =Apy {exp (—a,y)—exp(—agy)} sees (9) 

Following Heidmann (1950), the wave function y,(y) is taken to be a Hulthén 
S wave function, for a Yukawa potential. The normalizing constant is adjusted 
so that in limit 50, ¥,,(y) is asymptotic to the S part of a plane wave, normalized 
to the nuclear dimensions. 

sin (k3y +6) 
ee 7 o .2)—1/2 3 2 FR] 
baa ty Moet kyty 19 | PCIE) _ gw] 
ov Hsin (Ra V)ikaywonen OU, 9 we nes (10) 
where cot 6=—a/ks, «1 being the triplet scattering length and p+ the range 
of the nuclear force for a Yukawa shape. 

The integration over the nucleus may be extended over all space without 
loss of accuracy as p(y) decreases sufficiently rapidly with y, and hence 
G(k,, K,) may be computed by standard methods: 

* + q1")(%a" + 92") 
G(k,, K,) = 87A)K,- a+ hey | } cot 6 lo der iredy Na + ai) 
(Kay Ma) = Sarin Ra =) i 2 (as? + 4a? )( 2? + 91°) 
(G1 + Jo)(% = 1) (1% + 9192) 
(041? — G1a)(%2” — 9192) + (Gi + Go)? %9 
= At G)(%1— %s) 
+4 tance ee eee 11) 
Pa HF aU 0) + (FP ae 
with q)=}Ky+hy qo=3Ko—hs. 


(ii) The Functions F??(K,) and F"*(K.) 
If o(K,) is the differential cross section in the centre-of-mass system for a free 
collision (in Born’s approximation) between unlike particles of relative momentum 
K,, then 


o(K,) dQ =(m/47h?)?| | V(x) exp {—iK, . x} dx |? = (m/4rh?)? | F(K,) |?. 
In accordance with the remarks above | F”?(K,) |? is replaced by 167?(h?/m)? x 
(the experimental cross section for n—p scattering at the appropriate energy) and 


+4tan 
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| F™(K,) |? is replaced by 167?(h?/m)? x (the experimental ‘p—p scattering cross j | 
section), in default of direct information about the n—n collision. 


(iii) The Energy Relations 


The energy of the nucleus A, ¢, differs from that of nucleus B, ¢g by the 
potential and kinetic energies of nucleons 1 and 2 in A, so 


€, —€p =(22/2m)(R, 24 Ry?) — 58 Mev 
where nucleons 1 and 2 are assumed to be eee in a uniform potential well 
of depth 29mev. ‘Then, as 


E,— Ep=e,— ep + (h?/2m)(ko? — k,.”) — (h*/4m) Ky? — eq 
where ¢, is the binding energy of the deuteron, the energy conservation relation 
becomes 

- W/2m(—$K,? + 2k,. K, — K,? + 2k,?4+ 2K? — B) =0 
where the variables have been changed from k,, k,, k, to Ky, ks, K; and where 


B=58Mev+e,=56Mev if the binding energy of the deuteron is taken as 
approximately 2 Mev. 


(iv) Integrations over K,, k;, K, and K, 


To obtain the total cross section, expression (8) has to be integrated over 
all Ks, k, such that |k,|, |k,| <Z and all k,, K, such that |k,|,|4K,] >Z. The 
integrations over K, and the angular variables of k,;, K, and K, may be performed 
analytically by standard methods, with the help of the conservation of momentum 
(2K, =K,—K,) and energy relationships. 

The total cross section op(ky) is then obtained by numerical integration 
over k,, K, and K,, the limits being 


(a) 0<k;<L—}|K,—-K,? 
(0) Kee Ue koakere 
(c) K,22L, K,2<2{k)2+ L? —(2m/h?)B}. 


If o®? and o?? are the experimental n—p and p-p scattering cross sections then 
on(hq) =(7200/n*L°) | fo" (Ky) + o°*(Ky)} 
G®(K5, kg)b( Ky, Ko, 2p) Ky, Ko, Rg? dG aRedk, =) eee (12) 


where G is given by eqn. (11) and ¢(K,, Ky,k3) is a complicated factor arising 
from the angular integrations. 


§3. NUMERICAL CONSTANTS 

As mentioned above the same numerical values were used in computing 
from eqn. (12) as were used by Heidmann (1950). 

(i) For the deuteron ground state wave function: Ay*=0-0613 x 10% cm}, 
a, =0-237 x 10% cm, a, =1:55x10'%cm. For the initial wave function of 
nucleus 1 and 2 #,,(y), «1 =5-39 x 107% cm,  =0-847 x 1013 cm-1, 

(11) For the Fermi nucleus: The nucleons were considered to move in a uniform 
potential well of depth 29 Mev and the maximum of the Fermi sphere was taken 
as L=1:0 x10" cm-* 

(iii) For the functions F®°(K,), F®?(K,): As explained above, these functions 
were replaced by the differential cross sections o”?, o?? for neutron—proton and 
proton—proton scattering. 
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(a) o°? proton—proton scattering. ‘The scattering was taken as isotropic in 
the centre-of-mass system the differential cross sections being of magnitude 
3-9 x 10-*cm? at 340 Mev (Chamberlain et al. 1950, 1951), 5-0 x 10-%cm?6 at 
146 Mev (Cassels et al. 1951), 5-5 x 10-28 cm? at 105 Mev (Birge et al. 1951). 

(5) o®? neutron—proton scattering. The data of Hadley et al. (1950) at 
260 Mev and Hadley et al. (1949) at 90 Mev was used. (The total n—p cross section 
of Hadley et al. is in good agreement with more recent measurements at 95 Mev 
by De Juren and Knable (1950) and at 97 Mev by Taylor et al. (1951)). The 
total cross sections are 0-076 x10-?4cm? at 90mev and 0-035 x10-%cm?2 at 
260 Mev. 

The differential cross sections are roughly symmetrical about 90°. The 
forward direction only is of importance and the cross sections in this region 
are approximately of the form a/(b+ K,?). 

(iv) The Number w=}(A—Z)A: If Z is taken as A/2 then w=3A?. As 
explained below, the number of nucleons effective in deuteron production is 
less than A and a correction has to be made. 


§4. NUMERICAL RESULTS AND CONCLUSION 
The total cross sections computed from eqn. (12) are shown in the table, 
compared with those of Heidmann, for energies of 100, 200 and 300mMev. As 
eqn. (12) describes a kind of double process, the cross section depends on A?, 
in contrast with the direct pick-up cross section which is proportional to A. 


Energy (Mev) 100 200 300 
Indirect cross section (from eqn. (12)) (10-28 cm?) 3:6 A? 5:5 A? 4-5 A? 
Direct cross section (Heidmann 1950) (10-?8 cm?) 194 A 37 A 10:5 A 


'' The asymptotic variation of the cross section with energy may be deduced’ 
from eqn. (12) and the expressions for Gand ¢. If it is assumed that the variation 
with energy of the function F(K,) is small above 300 Mev, then the total cross 
section will decrease as E,!. ‘This is in marked contrast to the variation of 
the direct ‘pick-up’ cross section which Heidmann found to decrease with 
energy at FE, °. 

In comparing theory with experiment certain difficulties arise. The Fermi 
model used cannot be expected to give good agreement for nuclei such as carbon 
(for which the most complete experiments have been made), and on the other 
hand, for heavier nuclei the incident neutrons may make several collisions in 
passing through the nucleus. Another difficulty is that if deuteron production 
occurs inside the nucleus the deuteron will have a large probability of disintegrating 
before it can escape. Heidmann supposed that deuteron production occurred 
in a layer of thickness ry (the nuclear volume being 4779?A/3) on the escape side 
of the nucleus, thus reducing the effective number of nucleons from A to }A4?". 
The theoretical total cross section for carbon given by Heidmann was 
8 x 10-26cm?, about three times the experimental value of 2-6 10-*°; the 
discrepancy is probably due in part to overestimating the thickness of the 
deuteron producing layer. However, these considerations do not affect the ratio 
of the cross sections for the direct and indirect modes of deuteron production. 

Figure 1 indicates the cross section as a function of energy for a carbon 
nucleus, where the deuteron producing layer is of thickness 75. At 100 Mev the 
cross section for the indirect process is 11°% of the direct pick-up cross section. 
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This is satisfactory, as at this energy the experimental data may be explained in 
the main by the direct ‘pick-up’ process. At 300 Mev, however, the indirect 
process is already twice as important as direct pick-up and at higher energies 
predominates. oq 

It is characteristic of the indirect process that the production of deuterons of 
low energies is favoured, in marked contrast to the direct pick-up distribution 
where at high energies the deuterons are produced in a narrow band at 8/9. 
The energy distribution for E,=300 Mev is shown in fig. 2. The production 
of deuterons below 41 Mev is not considered, the cross section probably 
decreasing rapidly below this energy, as then the individual particles of the 
deuteron would have momenta less than the Fermi maximum. 


I, 
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Total Cross Section (10”cm*) 
uw 
A*x Cross Section per Mev (10™°cm?) 
a 
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Fig. 1. Total cross section for the production of Fig. 2. The energy spectrum of deuterons 
deuterons by neutron bombardment of a ejected by neutrons of energy 300 Mev. 
carbon nucleus (indirect process). 


To conclude, it may be said that while at the lower energies of the order of 
100 Mev, deuteron production can be understood in terms of a direct pick-up 
process, at higher energies, indirect processes play an increasingly important 
part. Although comparison of the magnitude of theoretical and experimental 
cross section is difficult, the mode of deuteron production may be distinguished 
easily by means of the energy distribution of the emerging deuterons. 
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ABSTRACT. ‘The angular distribution of alpha-particles of long range from the reaction 
UB(p, «)*Be has been studied over the range 130 to 280 kev proton energy. A strong 
anisotropy varying with energy is observed and can be explained in terms of an interference 
between states of opposite parity in %C. The theoretical analysis gives the result that the 
well-known 16-11 Mev level in "°C formed by resonant capture of 162 kev protons in 1B 
has spin 2 and even parity. 


§1. INTRODUCTION 
HE reaction of “B with protons gives a high yield of alpha-particles at 
comparatively low bombarding energies and was one of the first nuclear 
reactions to be studied. Cockcroft and Lewis (1936) bombarding a thick 
target at 300kev found a continuous distribution of alpha-particles with a 
maximum range in air of about 4cm, and a much less intense homogeneous 
group of about 4-4cm range. 

Dee and Gilbert (1936) showed that the observed spectrum could be 
produced by the competing processes 

BRE pee Bede, Po oe a raecn (1) 
and hip be ts ee meee, (2) 
followed in each case by break-up of the *Be into two alpha-particles. 

The present work is concerned only with the homogeneous group which 
leads to the ground state of *Be. No measurements have been made on the 
continuum. 

A resonance at about 160 kev bombarding energy in the yield of the long range 
group was discovered by Williams et al. (1937). The resonance has since been 
shown to be identical with that found for the 1'B(p, y)!2C reaction and has been 
studied repeatedly (Jacobs and Whitson 1941, Marvin 1945, Morrish 1949, 
Neuert 1937, 1939). Recent work by ‘angen (1946) and Morrish (1949) has 
established the natural width of the level and the absolute voltage at which the 
resonance maximum occurs. 

The angular distribution of the long range group has also been studied, but 
it proved difficult to make sure that none of the particles from the much more 
intense continuous distribution were being included. ‘The present work has 
shown that the continuum overlaps the long range group rather more in the 
forward than in the backward direction. Mica absorbers—as employed in 
almost all the older measurements—may not be very reliable for separating the 
two components, in view of straggling, etc. In the work to be described here 
a fast pulse-height analyser (Hutchinson and Scarrott 1951) and proportional 
counter was used. ‘This made it possible to observe the spectrum of 
alpha-particles while counting was in progress and to estimate the contribution 
of the continuum to the recorded count at each experimental point. 
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§2, EXPERIMENTAL TECHNIQUE 
2.1. General Arrangement 


The 1 mv High Tension apparatus of the Cavendish Laboratory was used to 
supply protons with energies between 100 and 300kev. The targets were thin 
films of ordinary borax (of about 4kev thickness for protons) evaporated on to 
thin aluminium foil and supported on a wire frame at the centre of an evacuated 
cylindrical tank. The beam of protons was limited to a diameter of $ in. by a 
molybdenum stop at the entrance to the tank. Two proportional counters, to 
detect the alpha-particles, were mounted on arms which could be rotated (using 
standard oil seals and Gaco rings) about an axis through the centre of the tank. 
Kovar-glass seals were used to bring in the counter leads, and the coupling 
condensers and pre-amplifiers were supported immediately outside. ‘The 
counter cases were maintained at earth potential. Figure 1 shows the general 
arrangement. 


| |] | | 
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Fig. 1. Diagrammatic view of apparatus. The cylindrical block A carries the counter shown. 
A’ carries the second counter, which is omitted for clarity. 


2.2. Counters and Target 


The counters were made from brass tubing of 2 in. internal diameter and 
their overall length was about 3}in. The wire, 0-002 in. diameter tungsten, 
was supported along the counter axis by two Kovar seals. A mica window of 
about 2 cm air equivalent with a circular aperture of about ? in. diameter was 
mounted on the side of each counter by sealing the mica between two brass 
rings with Araldite; no supporting grid was then needed. 

The counters were filled with slightly purified tank argon at 60 cm pressure, 
sufficient to stop all the alpha-particles from the "B(p, «)’Be reaction, and were 
normally used at 800 volts, giving a gas amplification of about 10. With no 
special precautions it was found that the counters gave peaks of about 200 kev 
width at half height for 5 Mev alpha-particles. 

Monitoring was done by keeping one of the counters at a fixed angle, usually 
150° to the beam. ‘The other counter was moved round to plot out the angular 
distributions. 

When thin foil targets were used there was no trouble from the deposition 
of carbon layers on the target surface, though solid backed targets did collect 
such layers. 


2.3. Counting Procedure 


The output of the amplifier associated with the moving counter was connected 
in parallel to the inputs of an ordinary discriminator and of the pulse-height 
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analyser. The visual display of this instrument showed the peak due to the 
long-range group at the end of a continuous distribution from the other 
processes. ‘The bias of the discriminator was related to the channel number of 
the pulse height analyser by injecting a standard test pulse into the pre-amplifier, 
and in an actual measurement the bias was set at a voltage corresponding 
approximately to the valley between the continuum and the group in the 
pulse-height distribution. Pulses of a height greater than this were counted on 
a scale of 100. It was felt necessary to adopt this procedure rather than to rely 
on the readings of the pulse height analyser directly, because in some of the 
measurements the average interval between pulses was shorter than the maximum 
dead time (1-2 msec) of the analyser. 

It will be seen from the plot of a typical spectrum (fig. 2) that it was not 
possible in practice to resolve the long-range group completely from the 
continuum. This was probably due to a combination of poor resolution in the 
counter and a genuine overlap in energy (caused in part by straggling in the mica 
window). In order to estimate the number of particles in the long-range group 
the spectrum was plotted on semi-log paper and the upper end of the continuum 
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extrapolated. The background thus obtained was subtracted from the total 
count recorded by the pulse-height analyser in the channels above the voltage 
level defined by the discriminator. ‘The result was used to correct the count 
registered by the scaler. (Plotting on semi-log paper allowed an almost linear 
extrapolation to be made.) ‘The background to be subtracted in this way 
amounted usually to 20° or less of the total count, and no substantial error is 
thought to have arisen from the process. A test showed that it produced a 
symmetrical peak for the long-range group even under quite unfavourable 
conditions, namely for a spectrum observed at a forward angle (see fig. 2). 
Readings were taken at angles between 30° and 150° to the proton beam. 
(The limits were set by the size of the counters.) The forward and backward 


748 D. M. Thomson et at. 


readings were taken with the target at 90° to the beam. In order to obtain points 
between 120° and 70° the target was rotated to 45° and the two distributions were — 
normalized to fit at 150° and 120°. 


The angular symmetry of the tank was checked by mounting on the target 


holder a ThC+C’ alpha-particle source of the same size as the target spot. 
No departure from isotropy was found greater than 1%, and no correction has — 
been made. 


§3. RESULTS 

Preliminary runs with a solid target backing confirmed the general result of 
Neuert (1937, 1939) that the angular distribution is markedly more anisotropic 
near the resonance than at higher energies. It was also found that the 
distribution appeared to contain a term in cos 0, and this was confirmed when 
foil targets were used, allowing measurements as far forward as 30°. Detailed 
angular distributions were taken at the resonance and at about 10 kev higher 
bombarding energy. The results, corrected for centre-of-mass motion, are 
shown in fig. 3. The solid curves have been fitted using the method of least 
squares, and the numerical values of the constants in expressions of the form 
I(@)=1+acos 0+6 cos? 6 are shown below each curve. Corrections to take ~ 
account of the finite angle subtended by the counter window were found to be 
always less than 1°% and have therefore been ignored. 

Following the discovery of this definite evidence for the existence of a cos @ 
term in the distribution, the magnitude of which varies with the bombarding 
energy, the technique of Devons and Hine (1949) was then adopted and the 
ratio of the yields at 150° and 30° measured over a range of voltage from 130 kv 
to 300kv. The results corrected for centre-of-gravity motion are shown in 
fig. 4 together with an excitation curve. It was not possible to obtain reliable 
results below 130kv because of the instability of the h.t. set and the extremely 
low yield; above 300 kv the long-range group becomes almost swamped by the 
continuum at the forward angle and the subtraction process becomes unreliable. 
It will be seen that the curve has a form very similar to that obtained by Devons 
and Hine near the 440 kev resonance in the 7Li(p, y)’Be reaction, though the 
voltage scale is less extensive. 

An uncertainty in all the readings has been in the absolute value of the 
bombarding voltage. ‘The h.t. set uses a high voltage (30kv) discharge tube 
ion source for which the initial energy of the beam depends on the rate of gas 
inflow. It was found that this caused the apparent resonance voltage, as read on 
the high voltage meter, to vary by as much as 10kv from day to day. It does 
appear, however, to have remained constant over the period of any one run, and 
each day’s readings have been adjusted to give the resonance voltage the same 
nominal value. 


§4. THEORY 
4.1. Preliminary Discussion 
‘The anisotropy of the alpha-particles at resonance implies that protons of 
/=1 or higher must be effective in forming this sharp state. Considerations of 
barrier penetrability at such low energies make it likely that 7=1 is in fact the 
correct choice, as the yield observed by Bowersox (1939) appears too large to 
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admit of J/=2 or more. The variation of angular distribution with energy 
(fig. 4) points to the interference of two states of opposite parity. We therefore 
suppose the broad state to be formed by protons with /=0. 
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Fig. 4. (a),Observed thin target excitation curve at 150°. (b) Variation of I (150°)/J (30°) with 
energy (same energy scale as (a)). The full curve is obtained from the theoretical analysis. 


We shall assume that "B has spin 3/2 and odd parity. This is in accord 
with experiment (Gordy et al. 1948) and with the predictions of the shell model. 
With these assignments of spin and orbital momentum, we arrive at the 
following possible designations of the interfering states in #2C*: 


Resonant state (J=1): (0, +), (1, +), (2, +), (3, +) 
Broad state (Z=0): (1, —), (2, —). 


Both the resonant and the broad states undergo a branched «-decay to the 
ground state or to the first excited state of SBe. Both of the *Be states break up 
into two alpha-particles, and therefore have even spin and even parity. Moreover, 
it seems probable that one or other of these *Be states has zero spin, and the 
formation of a (0, +) state of 8Be by alpha-particle emission from !?C* can occur 
only if the latter has spin and parity denoted by (0, +), (1, —), (2, +), (3, —), ete. 
We can thus infer that the broad compound state of #C is (1, —) and that the 
resonant state is (2, +). Zero spin is of course excluded, since it would lead to 
an isotropic angular distribution at resonance. 

If the spin of the *C resonant state is indeed 2, it seems likely that the ground 
state of 8Be has zero spin, and that the excited state has spin 2. ‘The evidence 
pointing to this in the present reaction is of two kinds. 

(a) The transition from the *C resonant state to the ground state of *Be is 
much less probable (~;,) than the competing transition to the first excited level. 
If we assign zero spin to the ground state, and spin 2 to the excited state, 
transitions to the former could occur only by d-wave alpha-particle emission, 
whereas s-wave as well as d-wave alpha-particles would be effective in the 


latter case. 
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(6) The strong anisotropy of the long-range alpha-particles at resonance 
would be surprising if the ground state of *Be has spin 2, since in this case one 
would expect that emission of s-wave alpha-particles would render the angular 
distribution predominantly isotropic. 

It must be admitted that neither of these arguments is very substantial, for 
the coulomb and angular momentum barriers do little to discourage d-waves as 
against s-waves in the alpha-particle emission. Our assignment is, however, 
favoured on theoretical grounds and also finds support in some recent work on 
the reaction *C(y, «) (Telegdi 1951). 


4.2. Detailed Analysis 


We now make use of the above conclusions to obtain the exact form of the 
angular distributions. The spins of 4B and a proton can combine to give a 
resultant j equal to 1 or 2. Only 7 =1 can be effective in forming the broad state, 
but both spin combinations can take part in forming the resonant state. The 
differential cross section at a given energy E can thus be written in the form: 


A & oe 2 
aie. > 4 m iB nit m 
o(8,B)~ S|, ms 1,012, m) cee Vo", #) +PB ee Yi, $) 
A et 
a5 sy ae £02 ———_—— Y,"(6 
Ep pe Ape ms | »m) OE Bad 2? a, ye 


Here p ={P,(E)/P,(£)}” where P)(£) and P,(£) are the barrier penetrabilities 
at energy # for s and p proton waves respectively (Christy and Latter 1948). 
(A, «) and (B, f) relate the amplitude and phase for the resonant state to that for 
the broad state after allowance has been made for resonance factors and the size 
of the barrier penetration factor. ‘The additional phase z that we have introduced 
for the broad state comes from our assumption that it has a resonance at a much 
higher energy. / is the relative probability that the resonant state should be 
formed from the initial spin configuration j=2. Inserting the appropriate 
values of transformation coefficients <¢j,m;/,0|J,m) (Condon and Shortley 
1935; ve aa spherical harmonics, we find 


o(8, B)~ ai +f)+(1—f) cos? A} + 3 a2p*+ ¥/10ap = cos («—B—8—7) cos 0. 


In ie expression we have used the following abbreviations: 
R(E) ={E-—E,)? +402)", a=B/A, 6(£)=tan {1 /2(E—E£,)}. 

The values of £, and [ are known from the work of Morrish (1949) and 
Tangen (1946). ‘The other parameters at our disposal in fitting the experimental 
results are f, («—f) anda. All these may be energy dependent, although we have 
fitted the results assuming them constant. It is reasonable to suppose that their 
variations are small, provided that we are far from any resonance for the broad 
state. The change with energy of long-range alpha-particle yield between 
200 and 500 kev (von Ubisch 1942) can in fact be explained in terms of barrier 
penetration alone. 


§5. COMPARISON WITH EXPERIMENT 
The value of f in the expression for the differential cross section was found 
from the coefficient of cos? @ at resonance, where the influence of the broad 
state is very small, and the term in a? can be ignored. The coefficient of cos? 0 
at resonance is then (1—/f)/(4+/). The mean of several determinations at 
resonance gave f=0-42 + 0-02. 
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The ratio of the differential cross sections at 150° and 30° is then easily found 
to be 
ZB) = 1 + a*p?k? — 1-84 apR cos ¢ 

1 + a*p?R? + 1-84 apR cos d 
where ¢(£) =a —B—8d—7. 

At a fair distance from resonance, such that (H—E,)>5I', we can assume 
that to the order of accuracy of the experiment R can be written as (E—E,) 
and ¢ remains constant. In fact, the error in R thus introduced, even at 
(E—E,)=5I, is only 1%, and ¢ changes by as little as 6° in going from 
(E—E,) =S5IT to (E—E,) = o. 

It is convenient to consider the function (Z+1)/(Z—1) in this region; 
using this function, a method of least squares may be devised to find the values 
of a and cos ¢ which agree best with the experimental points. 

By taking a suitable mean value of 6 in this region, we find a=0-022 + 0-002, 
a—B=60°+7°. 

With these values for the parameters the exact expression for Z(H) was 
calculated for the whole range of energies, from 120 kv to 300kv, and the result 
is shown as the full curve in fig. 4(). Its close agreement with the observed 
values may be taken as to some extent confirming the constancy of a. 

It is interesting to note that modern dispersion theory (Wigner and Eisenbud 
1947) predicts a definite value for «—f. It is determined by the phase angles 
of the complex penetration factors for incoming protons and outgoing alpha- 
particles. We have calculated this angle, using the tabulation of coulomb functions 
by Bloch et al. (1951), and find it to be about — 60° (with little dependence on the 
assumed value of nuclear radius). The curve obtained with the theoretical phase 
would resemble that of fig. 4(b), except that the cos @ term in the angular 
distribution would vanish a few kilovolts above resonance, instead of below as 
the experimental results suggest. [he discrepancy appears to be genuine, but 
the following points should be considered: (i) The definition of proton energy 
in the experiment was not very good, and an observed angular distribution 
represents an average over several times the width of the 162 kev resonance. 
Computation shows that this does not markedly affect the observed value of the 
phase; it could shift the curve of J(150°)/J(30°) in the desired direction, but 
only by about 2kv. (ii) The outgoing alpha-particle energy is near the top of 


_the potential barrier, and the theoretical phase shift may be very sensitive to the 


assumed shape of potential near the nuclear boundary. Bloch et al. (1951) give 
the phase retardations for a pure coulomb potential since it is usual, though 
perhaps not very realistic, to assume such a shape right down to the nuclear 
radius. One could thus argue that a comparison with theory is not profitable in 
experiments of this kind unless both ingoing and outgoing particles have energies 
well below their respective barrier heights. 

The shape of the curve of fig. 4(b) is determined in the main by a and «—f, 
and the value of I’ only affects the curve in the region of the resonance peak. ‘The 
experimental results cannot, therefore, be used to obtain an accurate value for I’, 
the resonance width. 


§6. CONCLUSION 
The main feature of the experiments described here is undoubtedly the 
demonstration of interference effects between states of opposite parity. It is 
difficult to draw any quantitative conclusions from the fitting of a theoretical 
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curve to the experimental points. It does, however, seem clearly established 
that the resonant state in 1*C formed at 162 kev proton energy has spin 2 and even 
parity. Although our analysis gives a fairly good fit to the measurements, its 
assumption that the interfering background is due to a single state of spin one and 
odd parity is not demanded by the results, and further measurements of the 
excitation function above 300 kv would be desirable. 
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Angular Distributions of Particle Groups from Deuteron 
Bombardment of Neon 


By R. MIDDLETON anp C. T. TAI 


Nuclear Physics Laboratory, University of Liverpool 
Communicated by H. W. B. Skinner; MS. received 8th February 1952 


ABSTRACT. The angular distributions of nine proton groups, the elastically scattered 
deuteron group and an inelastically scattered deuteron group from the 8 Mev deuteron 
bombardment of neon have been investigated using the photographic plate technique. 
‘The absolute differential cross sections for the production of the various groups have been 
calculated and the spins and parities of some of the states of the residual nucleus 24Ne 
determined. 
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Si INTRODUCTION 

N a previous communication (Middleton and ‘Tai 1951) the technique used 

to investigate the charged particles emitted when neon of normal isotopic 

constitution was bombarded by 8Mev deuterons was described. Also 
presented were the measured Q-values of the reactions ?°Ne(d, p)*4Ne, 
?°Ne(d, d’)?°Ne and *°Ne(d, «)!8F. In this communication the angular 
distributions of the protons, elastically scattered deuterons and inelastically 
scattered deuterons are presented. 


§2. PROCEDURE 

From the measurement of the total number of tracks belonging to a particular 
particle group in a well-defined area of emulsion and the known geometry of the 
scattering chamber (Chadwick et al. 1944, Rotblat et al. 1951) the relative cross 
sections of the various particle groups at different angles can be calculated. 
Since the intensities of most of the particle groups, particularly the elastically 
scattered deuterons, show a pronounced maximum at small angles it is necessary 
to make more than one exposure of different intensities to obtain a convenient 
track density at all angles. For our purpose two exposures were made, one with 
neon at a pressure of 18-6 cm Hg and a total beam intensity of 0-62 microcoulomb, 
the other, intended for measurement at small angles, with neon at a pressure 
of 10-0 cm Hg and a beam intensity of 0-25 microcoulomb. Comparison between 
the two exposures was made by measuring the numbers of elastically scattered 
deuterons within the same area at angle 30°. ‘The ratio of the exposures thus 
obtained differs by 7+3°% from the value computed from the above exposure 
data. ‘The weaker exposure was used for measurement of angles 20° and 25° 
and the mean value of the ratios of the two exposures was used to correlate the data 
with that of the other angles. 

Measurements have been made at fourteen angles, for each of which some 
3000 proton and deuteron tracks were measured within an area of about 10 mm’. 
Due to the excessive intensity of the elastically scattered deuterons at angles less 
than 30°, the major portion of the measurements were made neglecting tracks 
within the range limits of the elastically scattered deuterons. ‘The number of 
tracks contained in groups which were not entirely resolved was obtained by fitting 
each group with a gaussian distribution. 

Since the emulsions used were of thickness 100, there was a tendency for 
some of the long range tracks to escape either into the glass or out of the surface. 
To minimize this effect the positions of measurement were chosen so that the 
angle of incidence of the particles with respect to the emulsion lay between 
4° and 8°. The number and range of all tracks escaping were also recorded. 
These tracks were then assigned to the longer range groups in proportion to the 
number of tracks already contained in these groups. ‘This simple method of 
correction was adopted since corrections based on statistical considerations of 
scattering phenomena are both complicated and not justifiable for the small 
number of observed escape tracks. ‘The uncertainty in this correction for the 
four long range proton groups is believed to be less than 6°% of the total number 
of particles in the groups and for the remaining groups to be negligible. 


$3, RESULTS 
The results of the measurements of the angular distributions have been 
converted into the centre-of-mass system and in figs. 1 and 2 are shown those 
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obtained for the elastically and inelastically scattered deuteron groups. In 


figs. 3 and 4 are shown the distributions of the proton groups P, and P, 
corresponding to the formation of *4Ne in the ground and first excited states 
respectively. The unit of relative intensity is the same for all groups and 
corresponds to a differential cross section of 0-74 millibarn/steradian (see §4). 
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The probable errors shown are those of a statistical nature only. Other possible 
errors may arise from (i) the uncertainty in the ratio of the intensities of the two 
exposures which may affect the 20° and 25° data by about 3% and (ii) the 
uncertainty in the escape correction which is about 5%, for the proton groups 
Eytan and Pe 
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. §4. DISCUSSION 
(i) Elastically Scattered Deuteron Group 


The general behaviour of the angular distribution of the elastically scattered 
deuterons (fig. 1) is similar to those observed for 6-3 Mev deuterons scattered by 
oxygen and nitrogen (Guggenheimer et al. 1947) except for a slight shift in the 
position of the secondary maximum. It therefore appears likely that our angular 
distribution can also be explained by using the one-body theory of elastic: 
scattering (Bethe 1937). The shift in the position of the secondary maximum 
to a smaller angle and the low intensity at large angles is probably due to the 
wavelength associated with the collision being less than those associated with the 
collisions of 6-3 Mev deuterons with nitrogen and oxygen, and also the larger 
nuclear radius of neon. 

If we assume that the intensity of the scattered deuterons below 30° is: 
due to coulomb scattering alone then, by comparing our arbitrary unit 
of intensity with the calculated coulomb scattering cross section, we can 
obtain a conversion factor to correlate our unit with the absolute unit. In 
this manner we have determined one unit of relative intensity to be equivalent: 
to 0-68 millibarn/steradian. For the particle groups, arising from the 
bombardment of ?°Ne, the above conversion factor has to be corrected for the 
contribution of about ten per cent of ?*Ne to the elastically scattered deuteron 
group and is 0-74 millibarn/steradian per unit of relative intensity. 


(ii) Proton Groups 


The experimental angular distributions of nine proton groups from the 
reaction *°Ne(d, p)?4Ne have been determined and those corresponding to the 
formation of 71Ne in the ground and first excited states are shown in figs. 3 and 4. 
All groups show a marked increase in intensity at small angles. From the 
measured angular distributions and the conversion factor obtained we have 
calculated the approximate total cross sections for each group and the results. 
are given in table 1. ‘The cross section of the remaining proton groups not 
listed in this table is estimated to be about one quarter of the sum of the listed 
total cross sections. This corresponds to a total cross section for all the proton 
groups of about 0-5 barn. 

Table 1 
Particle group ee eee eee ee eee bo one meine lactic; 
Total cross section (millibarn) ie PAR 450) 21-297 138 41 54 98 


The theoretical angular distribution curves of the proton groups from the 
reaction ?°Ne(d, p)*4Ne have been calculated by Huby and Newns using a 
stripping theory which gives results nearly identical with those of the published 
theory of Butler (1950 a, b). Comparison at small angles of the observed and 
calculated distributions gives for the states corresponding to proton groups 
P,, P, and P; a unique assignment of the angular momentum /, transferred by 
the neutron to the corresponding state of ?4Ne. ‘The angular distributions of 
the remaining six proton groups can each equally well be fitted by two theoretical 
curves corresponding to different values of /,. ‘This ambiguity arises as a result 
of the poor statistics of the observations at small angles coupled with the fact 
that it was impossible to extend measurements below 20° due to the high 
intensity of the elastically scattered deuterons. Since °°Ne has zero spin and even 
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parity, the corresponding state of #4Ne has spin /,+4 and even or odd parity 
according to whether J, is even or odd. The results thus obtained are listed 


‘in table 2. 


‘Table 2. Spins and Parities of ?4Ne Energy Levels. (Where two sets of values 
are given, the one without the parenthesis fits experimental data better than 


the other.) 
Particle group Po Ea Ps P; P, P; Bs Re BS 
Energy (Mev) — 0:33.) 91:68,.5 92°79) S738 4-71 5-44 5-74 7-30 
In Kage 2 Oor1 0 2or1, Jor(2) don)” ZorG) Toa 
Spin of 4Ne: J¢ 3 or $ 3 or 3 4 4 Sor$ or} Zor4 § or 8 ord 
gor 3 Sort (Zorg) (Zor) ($ord) (for %) 
Parity =e +*  -+or— + +or— =(or+) —(or) -klor—)) —Cr 


The salient feature of these results is that for all the measured distributions 
the value of J, is less than or equal to 2. This indicates the absence of states of 
21Ne of spin larger than 5/2 below an excitation energy of 7-3 Mev (corresponding 
to group P,) or if such states do exist they are not significantly formed by 
7:8 Mev deuteron bombardment. 


It will be observed from figs. 3 and 4 that at large angles the theoretical — 


differential cross sections fall off to much lower values than are measured. 
Similar behaviour has been observed in the cases of (d, n) reactions at 8 Mev 
bombarding energy (El-Bedewi, Middleton and Tai 1951a,b). This is probably 
due either to an inadequacy in the stripping theory or to the reaction proceeding 
predominantly via a compound nucleus interaction at large angles of particle 
emission. The former possibility arises from the fact that the present stripping 
theory cannot be expected to predict the correct differential cross section for the 
particles emitted as a result of deuterons of low incident angular momentum. 
Such nearly head-on encounters would be expected to cause the emergence of 
particles at large angles, but would not contribute significantly to small-angle 
particle emission. Similar conclusions seem indicated by the work of Cohen 
and Falk (1951). 

It is interesting to note that both the ground and the first excited state of 24Ne 
have the same value of /,(=2) for the absorbed neutron and may therefore form 
a spin doublet of ?D5). 3). states. It follows from the stripping theory of Bhatia 
et al. (1952) that for two states for which the wave functions are believed to be 
similar the total cross sections are roughly in the ratio of the statistical weights. 
Applying to the present case, we find, from table 1, that the ratio of cross sections 
of Py to P, is 13/21 which is very near the ratio 4/6 of the statistical weights of 
Ds), and D5. states. ‘Therefore we are justified in concluding that the ground 
state of #4Ne has a spin 3/2 and the first excited state a spin 5/2. This is in 
agreement with the value 3/2 given for the spin of the ground state by observations 
of hyperfine structure (Koch and Rasmussen 1949). It disagrees with the 
prediction of a D5). state for the ground state of 7!Ne from the Mayer’s shell model. 
A similar disagreement has been reported in the case of *8Na (Mayer 1950) 
which contains an equal number of particles in the outer shell. This may be 
explained by assuming that three particles in a state D5). can couple to give a 
spin 3/2, rather than 5/2, which is the value on the principle, suggested by Mayer, 
that the spin of even numbers of particles in a shell balance out, leaving wate! 
the spin of the odd particle. 


Deuteron Bombardment of Neon GWE 


(iii) Inelastically Scattered Deuteron Group 


The angular distribution of the inelastically scattered deuteron group from 
°Ne (fig. 2) shows a resemblance to that from ®4Mg as observed by Holt and 
Young (1949). The total cross section for the process is calculated to be about 
98 millibarns. The similarity between the appearances of the angular 
distributions and the order of magnitude of cross sections for the inelastically 
scattered deuteron group and for the proton groups suggests the interaction 
concerned to be of a similar nature. Theoretical calculations similar to those for 
the (d,p) reactions have been carried out by Huby and Newns (1951). 
Comparison shows the measured distribution to be consistent with the calculated 
one when the incident deuteron is assumed to impart one unit of angular 
momentum to the ?°Ne nucleus. Since the ground state of ?°Ne has spin zero 
and even parity, the excited state of 1-66 Mev would therefore be of spin 0, 1, or 2 
and of odd parity. Since there is a y-ray transition from that state to the ground 
state (Jelley 1950) the possibility of spin 0 is eliminated. 

It is also interesting to note that the appearance of the angular distribution 
for inelastically scattered protons from ?°Ne at an energy of 4-2 Mev (Heitler et al. 
1947) shows similar features to that of the inelastically scattered deuterons 


observed in the present work. 


ACKNOWLEDGMENTS 
The authors wish to acknowledge their indebtedness to Professor H. W. B. 
Skinner for his continuous interest and helpful advice and to Professor J. Rotblat 
and Dr. J. R. Holt for some useful discussions. We also thank Dr. R. Huby 
and Mr. H. Newns for making the results of their calculations available to us 
and Miss B. Anderson for assistance in observations. 


REFERENCES 

BeTHE, H. A., 1937, Rev. Mod. Phys., 9, 179. 

Buatia, A. B., Huanc, K., Husy, R., and Newns, H. C., 1952, Phil. Mag., 43, 485. 

But er, S. T., 1950 a, Harwell Conference Report; 1950 b, Phys. Rev., 80, 1095. 

CHADWICK, J., May, A. N., PrcKavaNcE, T. G., and PowELL, C. F., 1944, Proc. Roy. Soc. A, 
183, 1. 

CoHEN, B. L., and Fak, C. E., 1951, Phys. Rev., 84, 173. 

E.-Bepewl, F. A., MIDDLETON, R., and Tat, C. T., 1951 a, Proc. Phys. Soc. A, 64, 756 ; 
1951 b, Ibid., 64, 1055. 

GUGGENHEIMER, K. M., HEITLER, H., and PowELt, C. F., 1947, Proc. Roy. Soc. A, 190, 196.. 

HeitTier, H., May, A. N., and PowE Lt, C. F., 1947, Proc. Roy. Soc. A, 190, 180. 

Ho tt, J. R., and Younc, C. T., 1949, Nature, Lond., 164, 1000. 

Husy, R., and Newns, H. C., 1951, Phil. Mag., 42, 1442. 

JELLEY, J. V., 1950, Proc. Phys. Soc. A, 63, 538. 

Kocu, J., and Rasmussen, E., 1949, Phys. Rev., 76, 1417. 

Mayer, M. G., 1950, Phys. Rev., 78, 16. 

MIDDLETON, R., and Tat, C. T., 1951, Proc. Phys. Soc. A, 64, 801. 

Rorsrat, J., Burrows, H. B., and Powe i, C. F., 1951, Proc. Roy. Soc. A, 209, 46r. 


758 


LETTERS TO; THE EDITOR 


Local Penetrating Showers from Water and Carbon 


A cloud chamber of diameter 30 cm and illuminated depth 15 cm with a 2:5 cm Pb plate 
-across the middle was operated for 44 months under a layer of water 90 cm thick. The 
‘bottom of the water was 20 cm above the top of the cloud chamber. The cloud chamber 
-was operated by two counter trays (each of area 750 cm?), one immediately beneath the. 
‘water and the other ina lead pile beneath the cloud chamber. At least a two-fold coincidence 
from each of the trays was required. Extensive showers were detected by a tray of area 
750 cm? placed 1:5 m from the cloud chamber. At the end of the run a second run was 
made for an equal time with an amount of carbon equivalent in mass and geometry to the 
oxygen in the water substituted for the water. 

An event was classified as a local penetrating shower if it had either (a) at least two 
penetrating particles visible in the cloud chamber, or (b) at least one penetrating particle 
plus at least one nuclear reaction visible in the cloud chamber. Events in which the 
extended tray was discharged were disregarded. 

The number of events in which various multiplicities of penetrating particles occurred 
are as follows: 


No. of penetrating particles 1 2 Br) A +t, © Og EO een Total 
‘No. of showers from water 42° 256. 923" 2 OG ee —_ oie - 137 
No. of showers from carbon 52) 45, 21 Sn, a Siar tet ea 132 


The average multiplicity of penetrating particles from carbon is 2:25+ 0-13 and that for 
showers from water is 2:16+0:13. If we assume geometric cross section for interaction 
for carbon, oxygen and hydrogen and that the multiplicity from the oxygen of the water is 
equal to that from the carbon, this gives an average multiplicity for showers from hydrogen 
of 1:88+ 0-68. 

The average multiplicity of the showers from carbon is in good agreement with the 
results of Salant et al. (1950) working with photographic plates and primaries of similar 
energy; it also agrees well with the results of Chang and Del Castillo (1951) working with a 
cloud chamber; it is however only about half the value of the average multiplicity for events 
from carbon obtained by Walker, Duller and Sorrels (1952). It seems very likely that this 
is due to differences in geometry and selection. On the other hand the low multiplicity for 
the events from hydrogen is in good agreement with the conclusion of these workers that 
showers from hydrogen are of low multiplicity. This disagrees with the conclusion of 
Bertolino et al. (1952) that the average multiplicity in nucleon—proton collisions is of the 
order of 10. However it seems likely that this last result was due either to the simultaneous 
arrival of several nucleons as part of an extensive air shower or to the cascading of secondary 
particles in the lead shielding. 


Dublin Institute for Advanced Studies. C. B. A. McCusker. 
23rd June 1952. D. D. Miiar.* 


BERTOLINO, G., Cin1, M., COLOMBINO, P., WaTAGHIN, G., 1952, Nuovo Cim., 9, No. 5, 407. 
Cuanc, W. Y., and Det CastiILLo, G., 1951, Phys. Rev., 84, 504. 

SALANT, E. Oe ” HoRNBOSTEL, Nes Fisk, (on B., and SMITH, bi E, 1950, Phys. Rev., 79, 184. 
WALKER, W. D., DuLter, N. M., and SorrELs, J. D., 1952, Phys. Rev., in the press. 


* Now at University of eater 


Resonance Effects in the Photo-Disintegration of ‘*O Nuclei 
into four «-Particles 


The disintegration of '*O nuclei into four «-particles by high-energy y- rays has been 
studied by several workers. The reaction is 


"’O-+hv + 4 *He—14:5 mev (1) 


0 fe esl tele 


and the data relate to y-ray energies up to 23 Mev. In this energy region, stars were found 
with an average total w-particle energy of about 7:5 Mev, corresponding to y-rays of 22 Mev 
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energy. It was also reported that about half of these stars were apparently produced by the 
disintegration in flight of *Be nuclei in the ground state. Goward and Wilkins (1950) 
suggested that these *Be nuclei were formed from an excited level in the 12C nucleus at 
9-7 Mev, so that the entire reaction may be represented as three processes in rapid succession 
6’O-+hv + *He+2C* 
12C* _, ‘He+®Be 
®Be + 2 4He. 

We have undertaken further experiments to investigate first, the nature of the excitation 
function of reaction (1), and second, the detailed mechanism of the disintegration and the 
possible excited levels of intermediate *Be and C nuclei. 

Ilford C2 nuclear plates of 100 emulsion thickness were irradiated in. a synchrotron 
beam of maximum energy c. 32 Mev (University Department of Radiotherapeutics. 
Cambridge). 

The results obtained refer to 83 stars and the histogram (fig. 1) shows the number 01 
stars plotted against the y-ray energy. Although the number of stars is not high and 
although the histogram is not corrected for the spectral distribution of y-rays from the 
target, the results strongly suggest that the stars fall into various energy groups. The first 
of these corresponds to the distribution reported by Goward and Wilkins. The shaded 
histogram, showing stars apparently involving *Be nuclei in the ground state, agrees with 
previous results, in that 21 out of 54 stars in this low-energy region fall into this category. 
A more sharply defined group of stars occurs close to a y-ray energy of 25 Mev and there is 
some evidence of a further group at 29 Mev but this is near the maximum energy of the 
y-rays used. In this respect the reaction shows significantly different properties from those 
other photo-reactions of the (y, p) and (y, n) types which have already been investigated. 
Up to the present, the results already published on oxygen «-stars show no group structure 
such as we have found. In some cases this certainly is due to the fact that the maximum 
energy of the photons was too low (e.g. Goward and Wilkins 1950, Millar and Cameron 
1950). 

The occurrence of ®Be in the ground state seems to be confined to the low-energy group 
near 22 Mev y-ray energy. In a study of stars of this kind it is easy to recognize the two 
a-particles coming from the ground state *Be. ‘The other two «-particle tracks may then 
be used to calculate possible values for the excitation energy of a !2C nucleus as parent to 
the *Be. Only one half of these results can be significant; therefore fig. 2, which shows the 


10 10 
5 5 
ui 
» Jt OD 25 3, %s i i5 
7- ray Energy (Mev) "C Excitation (Mev) 
Fig. 1. Yield curve. Fig. 2. Possible excitation energies of !2C 


producing *Be in ground state. 
numbers plotted against excitation energy of #C nucleus, should exhibit a background 
continuum with significant peaks superposed, if the reaction proceeds via an intermediate 
12C nucleus. Interpretation of these results is rendered difficult by the experimental errors 
in determining ranges and angles of very short «-particle tracks. Nevertheless it is clear 
that, although 10 out of the 21 stars can be interpreted in terms of a level at 9-7 Mev, the 
remaining stars cannot be explained in this way. An alternative level near 11-3 Mev would 
account for the other stars. ‘This deduction is confirmed by theoretical calculations 
predicting the appearance of stars of this kind. If we take a star produced via a 9-7 Mev 
level in °C by a 22 Mev y-ray it will in general have one long a-particle track of energy 
3-9 mev, and three shorter tracks with energies totalling 3:6 Mev. Such stars are in fact 
observed; but in others the two a-particle tracks, not derived from the *Be in the ground 
state, are of comparable length, and these may be explained in terms of a C level near 
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11-3mev. The other stars in the low-energy group, not involving *Be in the ground state, 
are apparently produced from the 3 Mev level in *Be, 

An analysis of that group of stars produced by y-rays of energy about 25 Mev shows that 
it is improbable that the disintegration goes via an excited state of *C but there is evidence 
that an excited state in ®Be at 4-3 Mev is involved. In the high-energy group the number of 
stars is too small for any detailed analysis to be significant. 

If these deductions are confirmed, it would seem that the photo-disintegration of *O 
into four a-particles is, at least partly, a resonance process for y-rays below 30 Mev, and that 
in the lowest energy group there are alternative modes of disintegration via *Be. 


Oundle School, D. L. Livesey. 
Northants. 
Department of Radiotherapeutics, C. L. SmirH. 


University of Cambridge. 
7th May 1952. 


Gowarb, F. K., and WILKINS, J. J., 1950, Proc. Phys. Soc. A, 63, 1171. 
Mitiar, C. H., and Camsron, A. G. W., 1950, Phys. Rev., 78, 78. 


The Nuclear Electric Quadrupole Moment of Erbium 167 


Recently Bleaney and Scovil (1951) reported paramagnetic resonance measurements on 
erbium ethyl sulphate diluted with lanthanum ethyl sulphate in the ratio Er : Lax1 : 200. 
Their observations were made at 1:2 cm wavelength and at temperatures down to 13° k; 
at this temperature, where the line width resulted predominantly from the short spin-lattice 
relaxation time, the hyperfine structure was just resolved perpendicular to the hexagonal 
axis of the crystal and only partially resolved parallel to the axis. The spectrum was 
analysed in terms of the usual Hamiltonian 

H =B{g HS +e (H,S,+H Spits 1+ BS letS,1)+PU2—WI+ 1} 
and values were given for the first four parameters and an upper limit assigned to P. This 
upper limit was based on measurements of the spacings between transitions of the type 
(M, m) > (—M, m). 

We have extended the experiments to 4° K at 3 cm wavelength. Because of the longer 
spin-lattice relaxation time at this temperature we have been able to resolve the spectrum 
completely in all directions. As a consequence it is possible to give more accurate values for 
A and B as follows: |A|=0-0052+0-0001 cm; |B|=0-0314+0-0001 cm-. 

A nuclear electric quadrupole term can be found from anomalies in the spacings of the 
Am=0 transitions (Bleaney 1951). This effect not only requires very precise measurements 
of the hyperfine structure, which is not always possible because of the line widths, but can 
be masked by higher order terms resulting from the nuclear magnetic moment. A second 
method of finding a quadrupole term is to observe the otherwise forbidden transitions 
corresponding to the selection rules Am=+1 and +2 for the nuclear quantum number. 
According to Bleaney these transitions are split into doublets by a nuclear quadrupole term; 
he gives formulae for this splitting and for the intensity of the lines. 

We have observed the lines of both types. |The great anisotropy in this crystal means 
that the accuracy of observations based on the Am= + 1 transitions is entirely limited by the 
precision with which the crystal can be aligned with its axis parallel to the magnetic field. 
There is no such limitation for the Am= +2 lines, which have greatest intensity perpen- 
dicular to the axis, because irrespective of the mounting it is always possible to rotate the 
crystal through the perpendicular plane and this direction can easily be found from the 
position of the maximum g value. 

From measurements of the doublet splitting of the Am=+2 lines we find that 
|P| =(30+ 3)x10-4cm~!. Insertion of this value in the intensity formulae gives a value 
which agrees with the measured intensity within the experimental error of 20%, thereby 
providing a rough check on the result. 

Elliott and Stevens (1951) have estimated the gradient of the electric field and 1/r*. 
They gave two possible ground states for the ion; but since then the uncertainty 
has been removed (Elliott and Stevens, in course of publication) in favour of 
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cos 6| +7/2)+sin 6| £5/2) with @~45°. Using these results we can give the nuclear 
electric quadrupole moment Q as defined (Mack 1950) by O=J p7(32*—1’) dr where 
Pr denotes the charge density for m;=J. We find |Q|=(10-:2+ 3) x 10-*4 cm®. 

We wish to make grateful acknowledgment to Dr. A. H. Cooke for his co-operation 
during these experiments. One of us (H.J.D.) is indebted to the National Research 
Council of Canada for financial support. ; 


Clarendon Laboratory, G. S. BOGLe. 
Oxford. H. J. Durrus. 
24th June 1952. Ee EXD scovir 


BLEANEY, B., 1951, Pail. Mag., 42, 441. 

BueaneY, B., and Scovit, H. E. D., 1951, Proc. Phys. Soc. A, 64, 204. 
EuuiotTt, R. J., and STEVENS, K. W. H., 1951, Proc. Phys. Soc. A, 64, 205. 
Mack, J. E., 1950, Rev. Mod. Phys., 22, 607. 


Pile Neutron Absorption Cross Sections of Lead 206 and 207 


At the request of B. B. Kinsey (Atomic Energy of Canada Ltd., Chalk River, Ontario, 
Canada) we have measured the pile neutron absorption cross sections of two lead samples 
which he supplied. One was a sample of natural lead and the other a sample of radio lead; 
both were measured on the G.L.E.E.P. oscillator by comparing them with boron, the thermal 
neutron absorption cross section of which was assumed to be 710 barns. The results 
obtained were o(natural lead) =162+5 mbarns, o(radio lead) =84+ 3 mbarns. 

Neither of the samples was analysed for chemical impurities, and therefore these cross 
sections come into the category B of Colmer and Littler (1950). 

The value obtained for the cross section of natural lead is very different from the one 
obtained by Colmer and Littler, namely 280+10mbarns. The sample used by them 
was remeasured, and the value of 280 mbarns was again obtained. A spectroscopic analysis 
was therefore made on the sample, and the presence of cadmium was detected. 
Mr. A. A. Smales then determined the cadmium content by a radio-activation method to 
be 30+10 parts per million. This involves a correction of 134+45 mbarns to the cross 
section, thus bringing it into agreement with the present measurement on natural lead. 

Mass spectrometer analyses have been made of the isotopic abundances in the two kinds 
of lead. ‘The analyses were made using a solid-source method with lead iodide as the source 
material. The results obtained were 


Isotope 204. 206 207 208 
°% abundance natural lead 1:-46+0:04 23:944+0:06 22:324+0:06 52-27+0-06 
°% abundance radio lead — 88:19+ 0-02 8-86+ 0-02 2:95+0-01 


We have not been able to find any published value for the activation cross section of 
lead 204, and have assumed that it has not been measured to date because it has a very low 
value. The activation cross section of lead 208 to produce 3:3 hour lead 209 has been 
measured by J. S. Levinger and is quoted by Way and Haines (1948); the value obtained 
was 0-45 mbarns. 

Therefore, taking o(Pb 204) =0, o(Pb 208) =0-45 mbarns, we derive 

o(Pb 206) =25+ 5 mbarns, o(Pb 207) =698 + 25 mbarns. 
The values obtained for the cross sections of lead 206 and 207 are not very sensitive to the 
- yalues assumed for lead 204 and 208; in fact, the activation cross section of lead 204 could 
be as high as 300 mbarns without changing the values for lead 206 and 207 by more than 
one standard error. 

The mass spectrometer analyses were made by Aitken and Palmer, and the cross section 
measurements by Littler and Lockett. 


Atomic Energy Research Establishment, K. L. AITKEN. 
Harwell, Berks. Dy a lInrieR: 
30th June 1952. Ee RLOCKETT: 
G. H. PaLmer. 
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Artificial Activity Induced in Carbon by Fast *He Ions 


Natural helium 3 was discovered by Alvarez and Cornog (1939) by using the Berkeley 
60 in. cyclotron as a mass spectrograph. Using a beam of *He ions, they were able to 


produce some 2} minute #°P from silicon. Since these experiments some lower energy — 


work has been done with *He, but the only use of *He as a bombarding particle in the 
cyclotron which appears to have been reported is a study of reactions of the (*He, w) type 
which were observed in nuclear emulsions by Lukirsky, Mescheryakov and Khrenina (1947) 
using ®He ions of 5:7 Mev energy. 

In similar experiments made recently with the 60 in. Nuffield cyclotron of the University 
of Birmingham (1952) ordinary atmospheric helium containing about 1-3 parts per million 
of ?He (Aldrich and Nier 1948) was fed into the ion source. The cyclotron was operating 
at its normal frequency (10-2 Mc/s) but at three-quarters of the normal magnetic field. 
The exact magnet-current required for resonance for >He was found by means of a tube 
carrying a fluorescent screen and terminated by a Perspex window at the outer end. 
Using this it was found that at 20 mev the width of the resonance peak was about 25 gauss 
(in 10.000). 

The identity of the particles accelerated was confirmed by the use of Ilford C2 plates 
placed successively at various radii inside the cyclotron tank, plenty of tracks being obtained 
in one-second’s exposure with an arc current of about one-tenth of normal. ‘The grain- 
density of the tracks showed that the particles were doubly charged and the ranges had 
the expected average values. 

Carbon in the form of pure graphite was bombarded for 10 minutes with the *He ions 
at a radius in the cyclotron corresponding to 20 Mev. The decay curve of the resulting 
activity is shown in curve I, and can be interpreted as that due to a mixture of activities 
of 20 minute and 10 minute half-lives in nearly equal initial strength. ‘These are the 
half-lives of 4C and !N respectively, which could be produced by the reactions 
12C(®He, «)4C+ 1-86 mev and #C(?He,-d)#N—3-55 or ®C@He, *H)8@N—2:-29. 
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Decay of activities produced in bombardments of carbon targets. 
Curve I. 10 minutes. Nominal *He energy 20 mey. 
Curve II. 50 seconds. Enriched helium. Nominal *He energy 18-5 mev. 


‘The experiment was repeated, using in the ion source helium enriched to 1/1000 *He. 
The decay of the activity resulting from a bombardment of 50 seconds with 18-5 mev °He 
ions is shown in curve II. Comparison of these results shows that an increase in amount 
of bombarding *He by a factor of about 800 increased the 'N activity by a factor of 200 
and the "'C by 600. ‘The latter agrees with the enrichment factor within experimental 
error and indicates that the whole of the !C activity is due to 7He. On the other hand, 
some '*N activity in the first experiment must arise from other particles accelerated in the 
cyclotron under *He resonance conditions; these have been shown, by experiments to be 
described elsewhere, to be deuterons. 

Bombardment of carbon with ordinary and enriched *He at energies of about 4:2 Mev 
showed that for this energy the only activity attributable to the He ions was that of UC, 
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Roughly 10 uc of 1'C was produced by 1paAsec of 18:5 Mev *He ions. The current 
could only be estimated from the a-currents with similar source conditions, and might 
therefore easily be wrong by a factor of two or three. 

If the **N was produced from "C, the integrated cross section for its production was 
about one-third of that for “C, while if it was produced from !°C it must have been thirty 
times as great. 

The author is grateful to Professor W. E. Burcham for many valuable suggestions and 
to Mr. K. E. A. Effat and Mr. G. L. Munday for their help with the counting. He is also 
grateful to Mr. W. Hardy for his help in finding the exact conditions for acceleration and 
for making the necessary bombardments. 


Department of Physics, J. H. FREMLIN. 
University of Birmingham. 
7th December 1951; in final form 13th June 1952. 
AvpricH, L. T., and Nikr, A. O., 1948, Phys. Rev., 74, 1590. 
ALVAREZ, L. W., and Cornoe, R., 1939, Phys. Rev., 56, 379 and 613. 
Luxirsky, P. I., Mescueryakov, M. E., and Kurenina, T’. I., 1947, C.R. Acad. Sci., U.R.S.S., 
yap Oar A 
University of Birmingham Cyclotron, 1952, Nature, Lond., 169, 476. 


The Deuteron Stripping Reaction with Aluminium 


'Fheories giving the angular distribution of the particles emitted in the (d, p) and 
(d, n) reactions based on the capture of one of the particles in the deuteron in a single-stage 
or ‘stripping’ process have been given by Butler (1951) and by Bhatia et al. (1952). With 
a beam of deuterons of approximately 8 Mev energy and using a triple proportional counter 
as the detector of emitted protons, we have carried out angular distribution measurements 
extending to the forward direction for proton groups from a variety of elements. The 
isotopes which have so far been investigated are the following: %’Li, *Be, 7%?**Mg, 
27 Al, 78Si, 92S, 4°Ca and 88Sr. We report here briefly the results obtained with the reaction 
27 Al(d, p)?8Al. 

We have measured the angular distributions of the proton groups (pp and p,) relating 
to the ground state and the first excited state of **Al. The measurements for the group pe 
are shown in fig. 1. ‘The curve is a theoretical one which could be derived from either 
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«of the two theories mentioned above with suitable choice in each case of the parameter 
representing the radius at which the neutron is captured, and with the unambiguous value 
.of zero for the orbital angular momentum /7 of the captured neutron. ‘Thus neutrons 
zare captured with zero orbital angular momentum to form one or both members of the 
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unresolved ground-state doublet of ?8Al (Enge 1951). Since the spin of *’Al is 5/2 the 
spin of one or both members of the doublet must have either the value 2 or 3. According 
to the nuclear shell model, the lowest unfilled neutron orbit in ?’Al is an s-orbit. Our 
results are consistent with a neutron being captured directly into an s-orbit to form the 
ground state of 78Al. 

The experimental angular distribution for the proton group p, (fig. 2) is unusual since 
it could be fitted only by adding together in suitable proportions two theoretical curves, 
one having 7=0 and the other having /=2. The curve drawn in fig. 2 was calculated on 
Butler’s theory with the ratio of the peak values of the two component curves with /=0 
and /=2 adjusted to be 1°85. 

Such a double curve could, of course, be due to the energy level in question being an 
unresolved doublet. In the present case such an explanation is unlikely in view of the 
magnetic analysis of Enge. However, since the target nucleus has a spin greater than 
zero, the vector addition rules for angular momenta allow transitions to a final state of 
definite spin to take place with more than one value of /. In the present case, since the 
angular distribution shows that / can have the value zero, the first excited state of "SAI 
must have a spin of either 2 or 3. The transition from 2’Al having spin 5/2 to a state of 
28Al having spin 2 or 3 can take place with values of / of 0, 2 or 4. The p, angular 
distribution shows the presence of the first two of these values. ‘The component curve 
of this distribution having /=0 has a peak value about 10% of that of the p» distribution. 
The ratio of the peak value of the component curve of the p, distribution having /=2 to 
the peak value of the py distribution is in agreement with the ratio given by Butler assuming 
the same probability of capture of the neutron by the nucleus in both cases. Thus the first 
excited state of ?8Al is formed by the capture of neutrons with zero orbital momentum 
with a probability about 10% of that for the capture of neutrons with two units of orbital 
momentum. If then we assume that the neutron is finally bound in the nucleus with the 
same orbital angular momentum with which it is captured, the results show that its final 
state in *8Al is not a pure d-state, but has a 10%, admixture of s-state. Bethe and Butler 
(1952) have recently suggested a search for transitions with double /-values in order to 
examine the purity of nuclear states. Bohr and Mottelson (1952) on the basis of a modified 
form of the independent particle model of the nucleus, have predicted the existence of such 
mixed states. It would seem possible that our case of ?8Al is an example. 


Nuclear Physics Research Laboratory, JseRe orn 
University of Liverpool. T. N. MarsHam. 
30th June 1952. 
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Note added in proof to paper entitled ‘‘ Calculation of Scattering Amplitudes ”’,. 
by G. J. Kyncu. 


Equations (2.9) and (2.10), using an integral representation for the wave function, are 
only true when f(k, : p) is the same for forward and backward scattering. In general, to. 


correct these equations f(k, : p) should be replaced by f,(k, : p)-+(R/p) fo(k, : p) when fo. 


and f, contain even and odd angular momenta respectively. 


Two other mistakes are that the factor (k?/2m%) should be replaced by (—1/472) in 


these equations and that df/R?dR in eqns. (2.7) to (2.10) should be multiplied by 4 and 
not divided by it. 
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Physical Properties and Analysis of Heavy Water (National Nuclear Energy 
Series, Vol. 4A, Division III), by I. KirsHensaum. Edited by G. M. 
Mourpuy and H. C. Urey. Pp. xv+438. (New York and London: 
McGraw-Hill, 1951.) $5.25; 45s. 


This volume is devoted to the experience obtained in the Manhattan Project of the 
complete isotopic analysis of water. 

The earlier investigations of Joliot, Halban and Kowarski established the suitability of 
heavy water as a moderator for an atomic pile. One facet of the multilateral attack on the 
technical problems involved in the release of atomic energy, organized under the Manhattan 
District of the U.S. Corps of Engineers, was concerned with the industrial production of 
large quantities of deuterium oxide. The design, development and control of this plant 
required reliable, complete isotopic analyses of water. The project eventually accumulated 
a unique body of experience in this field, and, incidentally, critically appraised the existing 
data on the isotopically homogeneous varieties of water, redetermining some of the less 
satisfactorily established constants. 

The first chapter deals with the physical data and includes a description of the new 
determinations of the specific gravity of deuterium oxide of normal oxygen content made by 
Kirshenbaum, Graff and Forstat and the confirmatory independent determinations of 
Voskuyl and Barach. ‘The second chapter compares the experimental and calculated values 
of the equilibrium constants for the principal protium—deuterium exchange reactions. 

The next two chapters, comprising nearly half the book, are devoted to the mass spectro- 
metric analysis of water. After a brief description of the constructional details, including 
circuit diagrams, of the Nier types I and II 3:4 mass spectrometers and of the similar 
instrument designed for the carbon dioxide estimation of the oxygen isotopes, follows an 
excellent account of the difficulties that beset these determinations. Sufficient data are 
tabulated to substantiate estimates of the errors arising from a variety of causes. The 
‘precautions necessary when analysing samples of very low protium content are emphasized. 
These include the memory effect in the mass spectrometer, and the effect of handling this 
hygroscopic material in the ordinary atmosphere. Several techniques of preparation of the 
sample, both for hydrogen and oxygen estimations, are compared. Particular attention is 
paid to the method of equilibrating gaseous hydrogen with the liquid sample. Besides the 
mass spectrometric method, several densitometric procedures are discussed. ‘These include 
pyconometric and the various float and falling drop methods, the latter being extended to 
cover low protium content waters. A very interesting continuous reading float gauge, 
suitable for recording the output from a deuterium separation plant, is described. Attention 
is drawn to the importance of careful pretreatment of the sample. ‘The significance of some 
of the earlier densitometric estimations of deuterium is lost because of failure to normalize the 
18O content of the samples. Normalization can be effected by equilibration of the water with 
excess sulphur dioxide. 

Refractometric analysis is also described in detail and a method of complete isotopic 
analysis by the combination of refractive index and density determinations is mentioned. 
The data presented show that the mass spectrometric and densitometric procedures give 
equally consistent results for relative isotopic analysis of the hydrogen in water but that 
absolute determinations involve much greater errors. Although this book may confidently 
be expected to become the vade mecum of those making mass spectrometric or densitometric 
isotopic analyses of water there would still appear scope for the development of new 
procedures. It is surprising that the measurement of the freezing point is barely mentioned 
although certainly convenient for rough determinations of the isotopic analysis of waters 
containing comparable amounts of protium and deuterium. ‘The large difference in the 
neutron scattering cross sections of protium and deuterium might form the basis of an 
analytical procedure for the protium content of the essentially pure deuterium oxide. 
Spectroscopic methods do not appear to have been exhausted. 

The concluding chapter deals with new and existing data on the natural abundance of the 
oxygen and hydrogen isotopes in various materials. ‘The distribution of the oxygen isotopes 
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is especially interesting since the nature of the principal fractionation mechanisms is not yet — 
clearly established. 

The book contains a mass of previously uagublished dua culled from reports of the 
Manhattan Project. It has been most carefully edited, contains less typographical errors 
and is more homogeneous in treatment than previous volumes in this series. A. G. MADDOCK. 


X-ray Analysis of Crystals, by J. M. Bijvort, (the late) N. H. KoLKMEYER and 
C. H. MacGttiavrey. Pp. xii-+304. (London: Butterworths Scientific 
Publications, 1951.) 50s. 


This book falls into three parts. Five chapters deal with methods, three with results, 
and these are followed by a set of appendices. Chapter 1 is introductory, dealing with 
diffraction by a single molecule and by a crystal, the plane groups and space groups. In 
Chapter 2 the derivation of the Laue and Bragg equations is followed by a brief description 
of the more important types of x-ray photographs. Here one reads: “‘ The absence of a 
symmetry centre in a crystal can readily be established from the unequal development of 
plane and complementary plane, or from the presence of the piezoelectric effect etc.” The 
reader is then referred to Chapter 4, where he reads: ‘“‘ Frequently, however, the macro- 
scopic symmetry is not completely known. The decision whether or not a direction is 
polar is then often uncertain from crystallographic considerations, i.e. from the development 
of the faces, piezoelectric behaviour, etching patterns, etc.” In Chapter 4 the space group 
Pmnb is assigned to mercuric chloride without further comment because it forms bipyramidal 
crystals. There is no mention in the book of Wilson’s work on this important subject. 

Chapter 3 deals with the factors affecting the intensities of diffracted beams, the whole 
treatment being based on the powder photograph. Chapter 4 is devoted to the determin- 
ation of structures by the ‘ trial’ method. No less than ten pages are given to the deter- 
mination of the structure of HgCl, (carried out in 1934) from rotation photographs. Struc- 
ture determination by Fourier methods is discussed in Chapter 5. The impression is 
given in these chapters that the older method (comparison of intensities) is a trial method 
whereas the Fourier method is in some ways more direct, an impression heightened by the 
emphasis placed on the use of pairs of isomorphous crystals. Not everyone would agree that 
‘“* It is this method which now makes it possible to determine the structures of complicated 
organic compounds.” Patterson techniques are only mentioned in Chapter 5, the reader 
being referred to Appendix 6. 

Chapters 6, 7, and 8 contain a well-illustrated survey of the crystal structures of some 
inorganic and organic compounds, a survey useful for the general reader but obviously 
inadequate for the serious student of chemistry. Chapter 8 contains a number of minor 
blemishes such as the peculiar reference to the “‘ spidery model of the molecule of cyanuric 
triazide (N;CN) ;.... which strikingly demonstrates the fragile structure of this explosive 
compound ”’, the omission of the molecules of SO, from Figure 146 c, and a thoroughly 
confusing discussion of bond lengths in certain aromatic molecules on page 198. 

Nine Appendices account for about one-quarter of the book. They deal with the point 
groups, space groups, indexing of photographs, the formulae of Fourier and Patterson 
syntheses, the reciprocal lattice, and miscellaneous topics such as disorder in layer structures, 
electron and neutron diffraction, and finally Grimm’s classification (1934) of structures in: 
terms of bond type. ‘The reciprocal lattice is first mentioned on page 248 in connection 
with the indexing of Weissenberg photographs (which is carried out with a special type of 
set-square, though it is mentioned that it is possible to use curves drawn on tracing-paper), 
but it is not discussed in detail until page 255. "Two pages (244-5) are devoted to indexing 
a rotation photograph without use of a chart (none are illustrated in the book), and though 
the reference to the classical paper by Bernal is given no fewer than three times (twice 
incorrectly) there is no clear statement of the standard way of indexing oscillation photo- 
graphs. 

In the case of a book which is translated from the Dutch it is justifiable to ask whether it 
contains much information not readily available in English texts or a presentation of out- 
standing value for teaching purposes. This question is particularly relevant at the present 
time, for no expense has been spared in the production of this handsome book, a fact which 
is reflected in the price. It is claimed in the Preface (by W. H. Taylor) that this book will 
be valuable to students and stimulating to their teachers. Without doubting the truth of 
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the second statement the reviewer feels that a student will not find a clear statement of how 
crystal structures are determined at the present time, witness the large amount of space 
devoted to powder and rotation photographs as compared with Weissenberg photographs, 
and the archaic methods of indexing. 'There is also an undesirable lack of precision in many 
statements. For example, if the number of molecules in a unit cell does not calculate to 
exactly an integer “‘ it is mostly due to errors in the gravimetric determination of density ’’. 
Or, with reference to space group determination, ‘‘ The further choice of the space group 
conforming to the crystal symmetry and the absences is most/y made with the help of auxiliary 
tables compiled for the use of x-ray crystallographers ”’ (reviewer’s italics). The reviewer 
suspects that this is unfortunately too often true and may account for many of the erroneous 
space groups to be found in the literature. A deeper understanding not only of space 
groups but also of the fundamental ‘ classical’ crystallography would be preferable ; this 
is only too obviously lacking in many of the ‘ x-ray crystallographers ’ produced today. 
Although parts of this book might be suitable as general reading for certain classes of 
science students, it cannot be recommended to the serious student of x-ray crystallography ; 
more detailed and more up-to-date treatments are already available in this country. 
A. F. WELLS. 


Perturbation Methods in the Quantum Mechanics of n-electron Systems, by 
E. M. Corson. Pp. xii+308. (London: Blackie, 1951.) 65s. 


No one can get very far in the theory or application of quantum theory without meeting 
many-particle systems, and being at that moment introduced to plenty of difficulties and new 
ideas not encountered in the more familiar but simpler single-particle problems. For 
example, we can no longer hope for exact solutions of wave equations, but must seek methods 
of approximation which, while not too complex, yet do justice to the physical characteristics 
of our system. And we have to take into account the symmetry properties of wave functions, 
which result from the indistinguishability of our particles and the operation of the Pauli 
Exclusion Principle. Dr. Corson has set out to provide a theoretical account of some of the 
ways which are open to us. ‘This is particularly necessary since no other such account is 
available in English. He has adopted the Dirac formalism rather than the Schrédinger one, 
with the result that matrices occur far more frequently than differential equations. There 
are certain advantages in this, for the vector model of Dirac and Van Vleck provides elegant 
answers to many problems which are dealt with much more clumsily by more conventional 
methods; and the density matrix is recognized at once to occupy the central role in all 
discussions of many-electron systems. The chief disadvantages are that, with our present 
knowledge, a good deal of preliminary analysis (about one third of the book) involving 
representation theory and group theory is necessary before any physical problems can be 
properly dealt with. But when the stage has been prepared, the play proceeds at a great 
pace, and perturbation theory, variation theory, Thomas—Fermi statistical theory, molecular 
valence theory, second quantization and S-matrix theory all fall into place. 

As may be anticipated, the book does not make easy reading. ‘The Dirac bra and ket 
notation does not lend itself to a neat appearance, and at times the argument seems a little 
abbreviated. But it is a scholarly book, worth serious reading and study by postgraduate 
workers. The absence of all but a few illustrative examples will appear to some as a great 
gain: to others, interested more in the relation between experimental results and their 
interpretation, it will appear as a minor deficiency. ‘The book is extremely free from errors, 
and its approach is one which recognizes difficulties without glossing them over. Parti- 
cularly in view of its unusual contents, it is a most welcome addition to the literature of 
quantum mechanics. C. A. COULSON. 


Diélectriques solides, by R. Jouaust. Pp. 84. (Paris: Editions de la Revue 
d’Optique, 1949.) 500 fr. 

This little book gives an elementary introduction to the physics of dielectrics. It is 
primarily written for engineers and does not presuppose any extensive knowledge of atomic 
physics. It starts with an introduction into the atomic and electronic structure of solids, 
which is followed by simple facts on polarization, conductivity and dielectric loss. Further 
chapters deal with the influence of surface layers, and with various types of dielectric 
breakdown. Naturally a book of this length can deal only superficially with such a large 
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subject. Thus, for instance, only seven pages could be allowed for dielectric loss, and the — 
author has not found it possible to go much beyond the phenomenological theories. In other 
chapters, too, the explanation of the various properties in terms of atomic structure is only 
indicated. 

As a whole, I think, the book will be valuable by showing to electrical engineers the 
importance of modern physics for the understanding of dielectric phenomena. H. FROHLICH. 


Isaac Newton, by S. 1. WawiLow. Pp. viii+214. 1st German Edition. (Berlin: 
Akademie-Verlag, 1951.) 8.90 DM. 


This book is the authorized translation from Russian into German of a biography of 
Isaac Newton written by a Russian scientist who, at the time of his death early in 1951, 
was President of the Academy of Sciences of the U.S.S.R. His own interests, which lay 
predominantly in the optical realm, had led him to analyse and translate into Russian the 
Opticks of Newton, the tercentenary of whose birth subsequently called forth this wider 
survey of the philosopher’s life and achievements. It was first published by the Russian 
Academy in 1943 as one of a series of scientific works apparently intended for the serious 
general reader; and this translation is based upon the second edition of 1945. 

Following a brief sketch of Newton’s boyhood and student years, six chapters (com- 
prising nearly two-fifths of the book) are devoted to his researches on light, in which the 
author recognizes the germ of all Newton’s subsequent discoveries and the mature 
expression of his characteristic philosophy of physics. It was from optics that Newton 
passed to gravitational theory by way of speculations on the role of the aether. Wawilow 
recognizes the contributions which lesser men such as Borelli and Hooke made to the 
Newtonian synthesis; and he explains the criticisms to which the underlying assumptions 
of the Principia have been subjected by Einstein. The later chapters deal more briefly 
with Newton’s work in mathematics and his controversy with Leibniz, his researches 
in chemistry, his theological writings and public life. There is a short classified bibliography 
but no index. The plates include eight likenesses of Newton (several of admittedly doubtful 
authenticity). 

Wawilow obviously took pains to acquaint himself with the unfamiliar historical setting 
of Newton’s career. His book is soberly written, free from bias, and even from the expected 
‘social-economic’ emphasis. The author has made considerable use of L. T. More’s 
biography of 1934, but he has taken account of more recent amendments to the traditional 
story. In the course of translation and re-translation the versions of passages quoted 
here have become noticeably free; also some errors have crept into the text. On page 5 
Bath should read Bate; the date of Rémer’s determination of the velocity of light (page 24) 
was 1676; the Royal Society list reproduced in fig. 13 relates to 1675, not 1671; something 
has gone wrong with the formulation of Kepler’s third law on page 99, and Meshemami 
on page 152 should presumably read Masham. a A. ARMITAGE. 


Ultraviolet Spectra of Aromatic Compounds, by R. A. FRIEDEL and M. ORcHIN. 
Pp. vi+52+579 diags. (New York: Wiley; London: Chapman and 
Hall,.1951.) 80s: 


This book consists mainly of a collection of solution spectra of nearly 600 benzenoid 
and heterocyclic compounds, partly taken from the literature and partly determined in 
the authors’ laboratory. It is primarily intended for organic chemists, but will also 
constitute a valuable reference work for the growing number of mathematical physicists 
interested in the theory of the spectra of complex molecules. A diverse range of structures 
extending alphabetically from acenaphthene to xylenol are covered, and while many of the 
compounds included are probably too complicated chemically to be of much interest to 
physicists, the collection provides a convenient source of information on the ultra-violet 
and visible light absorption properties of most basic types of benzene derivatives which 
will be widely appreciated. ‘The paucity of up-to-date reviews of reliable experimental 
data in this field has frequently been deplored in the past, and Friedel and Orchin’s book 
represents a most welcome and timely addition to the literature of spectroscopy. 

E. A. BRAUDE. 
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ABSTRACTS FOR SECTION B 


The Load Dependence of Rubber Friction, by A. SCHALLAMACH. 


ABSTRACT. It is shown experimentally that the load dependence of rubber friction can 
be explained as being due to the load dependence of the true area of contact between rubber 
and track if the surface asperities of the rubber are assumed to be hemispherical. 


The Frictional Properties of Plastics, by K. V. SHOOTER and D. 'Tasor. 


ABSTRACT. The frictional properties of a group of linear polymers. have been investi- 
gated at loads ranging from a few milligrams to several kilograms at slow speeds of sliding. 
At light loads the coefficient of friction tends to increase but at loads above 100 g the 
coefficient of friction is generally constant and is almost independent of the size and shape of 
the surfaces. From the results obtained at loads above 100 g a tentative theory to explain the 
origin and magnitude of the friction has been developed. The theory is similar to that of 
Bowden and Tabor for the friction of metals. When a plastic slides on a harder metal 
strong adhesion occurs between the surfaces and shearing takes place within the bulk of the 
plastic rather than at the interface. The frictional force is essentially equal to the product 
of the area of contact and the bulk shear strength of the plastic. The experiments also show 
that the area of contact is proportional to the applied load so that the plastic behaves as if it 
possesses an effective yield pressure which is constant. The coefficient of friction is thus 
primarily determined by the bulk properties of the plastic and is equal to the ratio of the 
shear strength to the effective yield pressure of the plastic. When a plastic slides on a softer 
metal shearing occurs within the metal and the friction is primarily determined by the bulk 
properties of the metal itself rather than those of the plastic. Examination of the surface 
damage produced during sliding shows marked transfer of the softer to the harder material. 

The frictional behaviour of polytetrafluoroethylene differs from that of the other plastics 
in that the adhesion during sliding is small and shearing apparently takes place preferentially 
at the interface. However, when this plastic slides on a sufficiently soft metal marked 
adhesion and transfer of the metal to the plastic occurs. 


The Transfer of Metal to Plastics during Sliding, by E. RaBiNowicz and 
K. V. SHOOTER. 


ABSTRACT. Investigations have been carried out on the amount of metal transferred 
when a radioactive metal is slid over the surface of a plastic. It is found that in every case 
metal fragments are transferred to the plastic in amounts that are of the same order of 
magnitude for the various metals and plastics examined. The close similarity of the results 
with those obtained with sliding metals suggests that, by a process analogous to the welding 
that occurs between metal surfaces, strong local adhesion occurs when metal and plastic 
are pressed and slid together. 


On the Dielectric Constant of the Water in Wet Clay, by L. S. PALMER. 


ABSTRACT. Some recent experimental values of the variation of the dielectric constant 
of wet clay with moisture content cannot be explained on the assumption that the clay 
particles and water molecules form a simple mixture. The results can, however, be theoreti- 
cally deduced by considering relatively dry clay to consist of closely packed water-coated 
clay particles in an air matrix and relatively wet clay to consist of clay particles uniformly 
distributed in a water matrix, the water associated with the clay having an effective dielectric 
constant which varies exponentially from that of ‘ bound ’ water (about 3) to that of ‘ free ’ 
water (about 80) as the percentage of water is gradually increased. The conclusion follows 
that water films have a very low dielectric constant compared with that of water in bulk, 
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Magnetic Viscosity under Discontinuously and Continuously Variable Field 
Conditions, by R. STREET, J. C. Woo.iey and P. B. Situ. 


ABSTRACT. For a ferromagnetic specimen exhibiting magnetic viscosity, the intensity of 
magnetization increases continuously with time under steady external field conditions. An 
account is given here of investigations to determine the influence on magnetic viscosity of 
changes in the magnetic field applied to the specimen, two cases of discontinuous and con- 
tinuous changes in field being considered. The latter case is of practical importance when 
direct measurements of magnetic viscosity are made using specimens with demagnetization 
coefficient other than zero, since the effective field acting in the specimen is then continuously 
variable in time. It is found that both cases can be adequately described by extending the 
previously proposed formal theory of magnetic viscosity in which it was supposed that the 
domain processes responsible for magnetization can be activated by thermal agitation. 

‘The required extension of the theory consists in assuming that the activation energy of the 
domains is a function of the field. The experimental results show that the changes in 
‘domain activation energy are linearly proportional to the changes in field strength, if the 
latter are small. Denoting the constant of proportionality by g, and considering a possible 
domain model, it is shown that q is related to the elementary volume of the material in which 
thermal agitation leads to activation. Measurements on alnico show that the linear 
dimensions of this elementary volume are of the same magnitude as the width of a domain 
boundary wall. Results of the variation of g and other parameters characteristic of magnetic 
viscosity as functions of the temperature of the specimen are also presented and it is shown 


that the observed behaviour is consistent with the present views on the metallurgical structure 
of alnico. 


The Threshold Gas Pressure Required to Sustain a Stable Arc in a Magnetic Field, 
by N. L. ALLEN. 


ABSTRACT. Measurements of the minimum pressure required for the operation of 
stable, pulsed, high-current arcs in hydrogen and air, in strong longitudinal magnetic fields 
have been compared with earlier results obtained with helium. The values obtained are 
shown to be approximately those predicted by theoretical reasoning from the earlier work. 
Some relevant processes occurring in the discharge are discussed. 


The Thermal Conductivity of Rare Gases, by W. G. KANNULUIK and 
E. H. Carman. 


ABSTRACT. ‘The thermal conductivity of the rare gases is determined in absolute 
measure by a ‘ hot wire’ method at selected temperatures in the range —183° c to 306° c. 
The experimental values of the thermal conductivity are shown to be in good agreement 
with the values calculated according to the recent theory of Hirschfelder, Bird and Spotz in 
the above range of temperature. The quantity e=K/ync, is shown to be independent of 
the temperature, but the value of € increases slowly with the molecular weight from a value 
2:43 for helium to a value 2:58 for xenon. 


On the Variations of Lattice Parameters of some Semiconducting Oxides, by 
L. D, BROWNLEE and E. W. J. MITCHELL. 


ABSTRACT. Measurements have been made of the lattice parameters of the sub- 
stitutional semiconductors Ni(Li)O, Fe,(Ti)O,; and also of reduced Mg,TiO,. The 
contraction of the nickel oxide lattice, which was briefly reported by Verwey et al. (1950), 
has been observed for a series of compositions. We have also observed the expansion 
of the Fe,O, lattice associated with the production of Fe® ions. Similarly, the expansion 
which accompanies the reduction of the spinel Mg,TiO, is attributed to the formation 
of Ti? ions. The lattice parameter of stoichiometric nickel oxide is believed to be 
4:1726+ 0-0002 kx (4:1811 A) while that of unreduced magnesium titanate is found to be 
8-425+0-0005 kx (8-442 A), 
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The Formation of Order in the Alloy AuCus, by 1. G. EomMunps and R. M. HINpeE, ‘i 


ABSTRACT. 'The approach to order in AuCu, has been studied by x-ray examination 
of single crystals. The distribution of intensity in the diffuse super-lattice reflections was — 
measured from a series of photographs of a stationary crystal, and again from a moving-film ~ 
photograph, the effect of instrumental broadening being eliminated by Stokes’ method. 
The intensity in the Aki reflection with (k+/) even is given as the Fourier summation 


I(uvw)= % X & K exp {— (Bx? + ay?+ a2") } exp 27i(xu+yu+ 2); 
UY z@ 


the constants « and f tend to zero as complete order is approached. The experimental 
results of Strijk and MacGillavry conform to this expression, and a measure of agreement 
is obtained with Cowley’s observations at a very early stage of ordering. 

The discrepancy between the present observations and crystal models proposed by 
A. J. C. Wilson is discussed, and expressions are derived for order parameters at successive © 
short ranges. It is suggested that the ordering process is similar at all stages, and depends 
on a tendency for gold atoms to avoid each other; ordered antiphase domains are formed as" 
a consequence of this process. ; 


Iterative Ray-tracing, by W. WEINSTEIN. 


ABSTRACT. An iterative ray-tracing method described by T. Smith is shown to give 
a divergent sequence of approximations in certain cases; when it converges it has only 
first order convergence according to Hartree’s classification of iterative processes. An 
iterative method which always gives second order convergence is described and the 
significance of the successive approximations in terms of higher order aberrations is 
discussed. 


